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SECTION 1 


INTRODUCTION 


A study at Williamstown , Massachusetts, in August 1969, jointly 
sponsored by the National Aeronautics and Space Administration 
and by the Massachusetts Institute of Technology* recommended; 

Develop navigation systems capable of economically 
determining ship velocities to +5 cm/sec over several- 
minute time averages and ship positions to +100 m at 
the same time, plus free instrument float locations with 
+2-km accuracy at intervals of about 5 days: to obtain 

better knowledge of current patterns and diffusion rates. 

In response to this recommendation and requests from other 
Governmental agencies for precision position fixing equipment, 
NASA Wallops Station by direction of the NASA Office of Space 
Sciences and Applications defined a program for the study and 
experiment described in this final report. This work was per- 
formed under NASA Contract No. NAS6-1962 by Applied Information 
Industries, using unique proprietary equipment. During the 
course of this contract effort the applications of precision 
position fixing became increasingly apparent to the individuals 
involved in the program. Requirements for position fixing with 
accuracies of less than 50 feet have been expressed by all three 
Armed Services, the Coast Guard, the Maritime Administration, 
the Federal Aviation Administration, the National Oceanic and 
Atmo‘spheric Agency, and the Department of the Interior. Many 
areas of private industry are also most concerned; these include 
geophysical exploration teams, ship owners and operators and air 
lines . 

The test results described here verify the ability of a satellite 
based system to satisfy the requirements of precision position 
fixing. The potential for serving the user needs by using 
communications satellites in synchronous orbit and low-cost 
ground terminal receiving equipment is most attractive. 


* The Terrestrial Environment, Solid-Earth and Ocean Physics , 
Application of Space and Astronomic Techniques , August 1969, 
NASA ERC, Cambridge, Massachusetts 
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SECTION 2 


SUMMARY 


2.1 PROGRAM SYNOPSIS 

The L-*Band Precision Positioning Experiment, sponsored by the 
National Aeronautics and Space Administration, Wallops Station, 
encompassed seven months of experiment design, experimentation, 
and data reduction and analysis. In the experiment the ATS-5 
synchronous satellite L-band transponder was used in conjunction 
v/ith the modified All ALPHA II navigation receivers to demonstrate 
the technical capability of precision position fixing for oceano- 
graphic purposes. The test program proved -the feasibility of 
using relative ranging techniques implemented by two identical 
receiving systems, properly calibrated, to determine a line of 
position accurately on the surface of the Earth. As this final 
report shows, the program achieved the objectives established in 
the work statement; i.e., to demonstrate 

• Level of resolution 
o Repeatability 
9 Precision, and 
9 Accuracy 

of existing modest-cost effective navigation equipment. 

Figure 2-1 summarizes the major tasks and events that comprised the 
program. The experiment configuration and data reduction techniques 
were developed in parallel with the hardware modification tasks. 

In the hardware modification phase, two existing All L-band receiving 
systems were specifically revised for the test and then calibrated 
by operating with the All system modulator; also, a tape control 
unit was developed for interfacing the receivers with a magnetic 
tape drive for data recording. 

A PN code modulator was installed at the NASA STADAN Station at 
Rosman, North Carolina, to provide the ranging signal for trans- 
ponding by the ATS-5 satellite. 

After the installation of the receiving systems at the NASA Wallops 
facilities, three weeks of testing were conducted to evaluate the 
equipment precision and determine the accuracy of relative position 
fixing. The initial data reduction was accomplished in near-real- 
time using the NASA computer facilities at Wallops Station. 

The location of the average LOP for each test period was offset 
from the remote site, the location of which was to be determined 
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Figure 2-1 . Program Plan and Schedule 


relative to the position of the fixed site, by not more than 8.5 
feet. This offset is believed due to the inaccuracies in the 
satellite ephemerides. The rms noise was on the order of 30 feet. 
Lines of position (LOP's) resulting . from individual pairs of 
range measurements, one measurement resulting from each receiver, 
indicate that for oceanographic and geophysical purposes 20 
successive bursts from the ATS-5 satellite will yield an LOP which 
is offset from the true location by not more than 50 feet. These 
data were collected in 16 seconds. Since the duty cycle of the 
ATS-5 satellite is such that data can be collected for 50 milli- 
seconds every 780 milliseconds, an improvement in duty cycle, as 
would result from a nonspinning satellite, will permit the same 
accuracy to be realized in approximately one second of ranging on 
the satellite or an equivalent improvement in accuracy if data 
were smoothed over 16 seconds. 

2.2 INSTRUMENTATION CONFIGURATION 

The synchronous ATS-5 satellite is the only satellite of the ATS 
series with L-band transmission capability and thus was selected 
as the platform for the relative ranging test. Because the ATS- 
5 satellite is rotating, its directional antennas scan the Earth 
approximately once every 780 millise'conds . This results in 
common illumination of the transmitting and receiving terminals 
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for approximately 50 milliseconds every revolution. During the 
experiment period, the ATS- 5 L-band receiver was temporarily 
inoperative, thus necessitating C~band transmission from the 
Rosman site. Operating the satellite transponder in the C-band 
to L-band cross-strap mode, permitted L-band transmission from 
the satellite to the receiving systems located at the NASA 
Wallops site. However, the transponder bandwidth was limited to 
approximately 6 MHz when operating in this mode, thus degrading 
the 10-megabit code-modulated signal by approximately 6 dB, 

The All ALPHA II receiving system employed during the experiment 
was an L-band ranging system operating with a 10-megabit PN code. 
The PN code is transmitted as coherent phase-shift-key (CPSK) 
modulation from the NASA Rosman STADAN station, using an ALPHA II 
modulator. The 10-megabit code provides the capability for deter- 
mining the coarse range granularity to within 100 nanoseconds. 
Additional circuitry within the ALPHA II code correlator provides 
quantization to 1-nsec intervals. The ALPHA II correlation 
receiver employs a demodulation technique which does not require 
carrier acquisition prior to code lock-on. This accelerates the 
lock-on time and enables operation in the presence of the inter- 
mittent signal from the ATS-5 satellite. A digital memory exists 
in the range tracking loop of the code correlator to provide 
continuity between the satellite bursts. 

Figure 2-2 shows a simplified version of the test configuration 
for the relative ranging experiment. An All ALPHA II modulator 
was used to drive the NASA C-band transmitter, while the ATS-5 
transponder transmits at L-band to the two independent ALPHA II 
receiving systems, located at the NASA Wallops Station. A common 
magnetic-tape data recording subsystem provided a permanent record 
of the range measurement data. 

The data recording subsystem consisted of a tape control unit 
developed by All to interface the ALPHA II code correlator with 
a 7-track incremental tape recorder. The receiver at the fixed 
site was connected directly to the tape control unit, while the 
receiver located at the remote site interfaced with the tape 
control unit via a data modem and a conventional telephone line. 
Range readings were taken for each satellite burst and if the two 
receivers measured the range reading within 50 milliseconds of 
each other, the tape control unit would accept the readings as 
being valid and record them, with time information. The resulting 
data tapes were then reduced using the computer facility at 
the NASA Wallops Station. 

2.3 TEST CONFIGURATION 

The test configuration illustrated in Figure 2-2 was implemented 
at the Wallops facility as shown on the simplified site map of 
Figure 2-3. The fixed receiving . site was located according to a 
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Rosman, N.C. 


Wallops 
Main Station 


Wallops 

Island 


Figure 2-2, Test Configuration for the Relative 
Ranging Experiment 


first order survey and the range to the satellite measured by the 
fixed receiver and the remote receiver. The geometry is such 
that by knowing the location of the fixed receiver and the pre- 
dicted satellite position as a function of time and by measuring 
the ranges R;^ and Rg from the satellite with the two receivers as 
shown in Figure 2-4, the line of position (LOP) containing the 
remote receiver can be calculated. This line of position can then 
be compared to the surveyed position of the remote receiver to 
determine the precision and repeatability of the system. 

The range measurements include several biases and error sources 
which were evaluated prior to computing the line of position. 

A fixed bias representing the delay difference between the antenna 
feeds and the correlation point within the code correlators was 
calibrated. Rubidium standards were used to drive the PN code 
generators in each code correlator and these contributed what 
appeared to be the most significant noise and bias errors en- 
countered in the experiment. Calibration of the slowly varying 
bias term was accomplished by comparing the phase difference 
between the two rubidium standards before and after each satellite 
ranging test period and fitting the resultant points to a straight 
line. The magnitude of the frequency standard noise errors was 
determined by performing a series of special tests designed to 
measure the stability characteristics between the two rubidium 
standards . 
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Figure 2-3. Site Map Showing Relative Locations 
of Fixed and Remote Receivers 
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Errors contributing to inaccuracies in the ranging measurements 
arise from either equipment or environmental noise sources and 
include receiver noise, equipment non-linearity, quantization 
frequency-standard errors, multipath, ionospheric and tropospheric 
delays and geometric effects. Errors associated with satellite- 
position uncertainty represent the present tracking and prediction 
capability and have been reduced to second order effects by virtue 
of the method of data reduction used. The equipment noise error 
(which includes both the receiving system noise and the frequency- 
standard noise) and the satellite position uncertainty are the two 
largest error sources contributing to the relative range uncertainty. 
Table 2-1 summarizes the error sources and the total relative range 
uncertainty. 

The link budget for the experiment is summarized in Table 2-2. 

This reflects the use of one TWT in the satellite L-band trans- 
ponder and a 3-foot receiving antenna at the Wallops Station. 


TABLE 2-1. ERROR SOURCE SUMMARY 


Error Source 

Relative Range 
Uncertainty 

Receiver Noise 

5.1 feet 

Correlator Linearity 

2.8 feet 

Phase Alignment of Frequency 
Standards 

14.3 feet 

Satellite Position 
Uncertainty 

7.0 feet 

Other Error Sources 

0.8 feet 

Range Accuracy (RSS) 

. 17.0 feet 

Position Accuracy 

21.3 feet 


2.4 TEST RESULTS 

The average line of position (LOP) for each test period was 
computed from the individual LOP's. The standard deviation of 
the individual LOP's about the mean was then determined for each 
test period. Figure 2-5 summarizes the average LOP’s for each 
test period and Table 2-3 presents the offsets and rms noise for 
each test. 

A test period dedicated to the calibration of the receiver front- 
end delay yielded a bias of 45 feet between the two receiving 
systems which is equivalent to an LOP offset of 56 feet when pro- 
jected through the elevation angle. Calibration of the phase 


2-7 




TABLE 2-2 


SUMMARY LINK BUDGET OF ATS“5 SATELLITE 
OPERATING IN C-TO-L CROSS-STRAP MODE 


Parama te 


Net Antenna Gain (dB) * 

ERP (dBW) 

Net Propagation Loss (dB) 

Received Carrier Power (dBW) 

Received System Noise Temp 
(dB “ °K) 

Noise Power Density (dBW/Hz) 

Carrier/Noise Power Density 
(dB-Hz) 

Receiver Bandwidth (dB) 
Carrier/Noise Ratio (dB) 


*85-foot Transmit Antenna; 3 and 4-foot receive antennas 


^ Projection of Line of 
/ Sight from Receiver to 
Satellite 


Rosman tc 

3 Wallops 

6214 MHz 

1551 MHz 

Up-Link 

Down- Link 

+76.4 

+31.5 

+96.5 

+24.0 

-203.8 

-189.5 to -191.5 

• 

1 — 1 
c\ 
\ 

-147 to -149 

+30,7 

+28.2 

-197.9 

-200.6 

+106.5 

51.6 to 53.6 

+67.8 

+32 

+ 38.7 

19,6 to 21.6 


AVERAGE LOP’S FOR TEST 
PERIODS 2,5,8,12 and 13 


Remote Site 



Test # 


Offset 

6.9 feet 

8.5 feet 
1.2 feet 

4.6 feet 

4.6 feet 


Test Number 2 13 8 12 5 















TABLE 2-3. SUMMARY OF THE TEST DATA 


Test 

Date 

Offset of 
Average LOP (feet) 

Rms (feet) 

2 

July . 12 

6.9 

25.88 

5 

July 15 

8.5 

31.63 

8 

July 21 

1.2 

16.77 

12 

July 27 

4.6 

33.26 

13 

July 28 

4.6 

30.52 


drift in the rubidium frequency standards by fitting the pre- 
and post-test calibration data to a straight line resulted in a 
correction term of approximately 600 feet per hour. 

The noise on the LOP 's. was analyzed by comparing the data from 
several test configurations. This enabled a noise value of 12 
to 15 feet to be determined for both frequency standards. The 
noise error of an individual receiver was determined to be less 
than 4 feet by removing the above frequency standard noise from 
the results, 

2.5 CONCLUSIONS 

The following conclusions are warranted by the results of the 
study and experiment: 

o Conmiunications satellites in ' synchronous orbit can be 

utilized as transponders for relay of precision positioning 
signals, 

o Using those signals at L-band, positioning accuracies 
relative to another receiver of less than 30 feet are 
achievable even without accurate ephemeris data and on 
a moving ship. With an integration period of an hour 
or more, relative positioning accuracies of 8 feet are 
achievable. 

• The major source of errors encountered appears to be 
the rubidium frequency standard. 



\ 
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SECTION 3 


EXPERIMENT DESIGN 


3.1 EXPERIMENT TEST CONFIGURATION 

The objective of this study was to demonstrate the technical 
feasibility and level of accuracy of precise range measurements 
at L-band via synchronous satellite for application to oceano- 
graphic and geophysical position-fixing. The demonstration test 
was required to operate with the spinning ATS-5 satellite and 
interface with NASA-supplied equipment and facilities. The 
experiment design was, of course, dictated by the test objectives, 
the equipment capabilities, and the constraints, of the proposed 
test site. 

To accomplish the experiment objectives, a test configuration 
employing a fixed and a transportable receiver, and operating in 
conjunction with a NASA transmitter and the ATS-5 satellite, was 
designed {see Figure 3-1) . The demonstration test profile con- 
sisted of calibration of both receivers in close proximity, 
followed by ranging measurements through the ATS-5 satellite, and 
generation of the line of position (LOP) on the surface of the 
Earth upon which the transportable receiver was located. By 
knowing the precise locations of the two receivers through first- 
order survey markers, the repeatability and precision of the 
ranging measurements could be established. The use of two 
proximate stations enabled most of the errors caused by propaga- 
tion effects and lack of knowledge of the satellite's position 
to be cancelled. 

As shown in the test configuration of Figure 3-1, an All ALPHA II 
modulator was used to drive the C-band up-converter and trans- 
mitter system located at the NASA Rosman, North Carolina, station. 
(C-band transmission to the ATS-5 satellite was necessitated by 
the temporary inoperative condition of the ATS-5 L-band receiver 
during the test period.) The 10-megabit PN-code-modulated signal 
was transponded by the ATS-5 satellite while operating in the 
C-band to L-band cross-strap mode, which has a 1-dB bandwidth of 
6 MHz. The identical and independent L-band down-converters 
located at the Wallops facility translated the received signal 
to a 70-MHz IF frequency. The IF signals were then processed by 
identical code correlators and decoded by PN code generators 
operating from separate rubidium frequency standards . To accom- 
modate the data collection function, the range readings from the 
code correlator were transmitted to a tape control unit, formatted 
and recorded on 7-track magnetic tape. 
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3.2 EQUATIONS FOR RELATIVE RANGING 

The relative ranging technique used in this test is a method which 
yields a line of position of one receiver (B) relative to another, 
fixed or reference, receiver (A) . The information used consists 
of the ranges measured by the two receivers, and Rg, at the same 
instance of time, and a calibration curve describing the behavior 
of the B frequency standard relative to the A frequency standard, 
the satellite ephemeris, and the location of receiver A. 

The complete derivation of the equations for position determina- 
tion in a relative ranging environment can be found in Appendix 
B, and are summarized here. 

The position of a point on the surface of the Earth can be defined 
by latitude, longitude and distance to the center of the Earth. 
(Refer to Figure B-1 in Appendix B.) These coordinates are de- 
noted by 1^, X and R. The terms A, B and S (used as subscripts) 
refer to the fixed and movable receivers and the satellite, 
respectively. In order for a system using range measurements for 
positioning purposes to determine these variables, an equal number 
of satellites is required. Therefore, two of these variables 
(latitude ij^, and geocentric range R) must be assumed when only 
one satellite is available. 

The measurements and Mg are ambiguous measurements of the sum 
of up-link distance D and down-link distances and Dg. The 
sums D + D;!^ and D + Dg are on the order of 2.46 x 10° nsec. The 
code D is only 204,700 nsec long. Therefore, Ma and Mg represent 
the excess above the largest integer multiple of the code length 
not exceeding the distances Da and Dg, respectively. The geo- 
metric range difference Dg - Da and the measured range difference 
Mb - Ma are related by 

Dg - D^ = Mg Ma + At + B + 45 ft (3-1) 

where A and B are calibration constants for the rubidium frequency 
standards and the 45 ft term is the measured fixed bias due to a 
difference in front end delay of the receivers. 

From spherical trigonometry, the range from each receiver to the 
satellite is given by 

+ Rg^ - 2R Rg cosy (3-2) 

where 

cosy = cos4> cosing cos (Xg-X) + sini|^ sin4^ (3-3) 
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The procedure is now as follows , The quantity cos Ya 
calculated for the fixed site and predicted satellite position 
from equation (3-3) . is obtained from equation (3-2) while 

D 3 results from equation (3-1) . Next cos yb calculated from 
equation (3-2) using the known value of R 3 . This value, together 
with assumed permits the determination of Xg from equation 
(3-3)* A slightly different value of leads to a different Xg. 
Two such pairs of latitude and corresponding longitudes define 
the line of position. 

The accuracy with which Xg is determined and in particular the 
effect of the uncertainty associated with satellite position are 
considered in Section 4. It is shown that the uncertainty in 
satellite position, i.e., a 10 x 10 x 10 km box centered at the 
predicted satellite position, will introduce a maximum error of 
12 feet. If the deivation of predicted position from true 
satellite position is consistently in one direction, the differ- 
ence will introduce a bias factor. An alternate method of LOP 
computation is discussed in Appendix F. 

t 

3*3 DATA REDUCTION METHODS 

^ typical test period consisted of three distinct time periods : 

(a) Pre-calibration, 

(b) Ranging on the satellite, and 

(c) Post-calibration. 

Receiver outputs were compared in the satellite mode and in the 
local calibration modes. During the calibration modes, when the 
two code correlators were co-located, a laboratory modulator 
simulated the link from the Rosman transmitter site to the two 
code correlators. The modulator was driven by the frequency 
standard of receiver A and its code output fed both to correlators 
A and B. The initial range values measured by the correlators 
reflected the code positions that the correlators were required to 
search through, starting from an aribtrary value in the range 
counters. This starting value was a random nvimber within the 
length of the code, zero to 204,700 nsec (201,337.02 feet). Range 
measurements made during the calibration mode were used to model 
the phase and frequency drifts between the two rubidium standards. 
Interpolation of this curve at the times for which measured 
satellite ranges were available yielded a correction factor to 
those measurements describing the long-term drift of the standards 
The correlators also acquired their references during this pre- 
calibration mode ('.'reference" being the initial reading of the 
code generator at lock-on) . 

During the periods that ranging data with the satellite were 
collected, receiver B was located at the test site near Wallops 
Island, while receiver A remained in the laboratory as the fixed 
site. 


3-4 



3.3.1 System Calibration 


Site Layout 


The fixed receiver was located in a building designated F-160 on 
the main Wallops base. The transportable receiver was located in 
Building U-15 near the gate to Wallops Island. Table 3-1 summa- 
rizes the geodetic position coordinates of these two sites based 
upon the Fischer Earth Model. 

The indicated-altitudes are the sum of local altitudes of the 
sites above mean sea level and the antenna heights. 

Rubidium Standard Calibration 


The rubidium standard drift which occurred during any one of the 
test periods was determined using range data collected during the 
pre- and post-calibration periods. These data and the times at 
which the individual data points were collected, were fit to a 
parabola and to a straight line using a least squares fit. 
Evaluation of two special ‘3-hour continuous calibration runs and 
the individual pairs of calibration sets showed that the straight 
line fit was preferred. (See Section 4.) 

Meaningful calibration of the rubidium standards and evaluation 
of the ranging measurements on the satellite resulted only when 
both receivers maintained their code reference during the test 
period. Loss of reference occurred several times during trans- 
portation of the rubidium standards to or from the remote site. 

It was caused by an interruption in the power or frequency 
standard feeding the correlator and resulted in a different range 
number appearing in the range counters. Since the new range 
difference was not related to the previous range differences, the 
drift between the standards could not be properly modeled from 
the data taken. 


TABLE 3-1. GEODETIC POSITION COORDINATES OF 

RECEIVER SITES (FISCHER EARTH MODEL) 


Parameter 

Building 

F-160 

Buidling 

U-15 

North Latitude (deg) 

37.936409 

37.864341 

West Longitude (deg) 

75.473952 

75.507593 

Altitude Above Mean 



Sea Level (ft) 

43.52 

45.44 
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Front End Calibration 


When each system was switched from the local pre-*test calibration 
mode to ranging on the satellite, bias errors from antennas, 
receivers and cables connecting each pair of subsystems were 
introduced. 

Tests were conducted on July 16 to assess the additional delay 
difference due to these error sources. Thus, no absolute measure- 
ment was made on the individual antenna-to-receiver links; rather, 
the difference in delays through these two links was determined. 

During those tests, the antennas were co-located on the rcpof of 
Building F-160. Care was taken to insure that the elevation and 
azimuth headings of the antennas were identical and that the 
antennas were aligned along the same perpendicular (LOP) to the 
line-of -sight to the satellite. A regular test period was then 
executed. The data collected was analyzed in the same manner as 
data resulting from other test periods. When the straight-line 
fit to the calibration data was removed from the range differences 
observed on the satellite, a residual of 45 feet remained which 
represents the delay difference between the two receiver front 
ends. This fixed bias was removed from all satellite ranging 
data so that the same internal points on the code correlator were 
referenced throughout the analysis. 

3.3.2 Computer Programs 

The range data generated by the two receivers during the entire 
test period was recorded by an incremental tape recorder (Kennedy 
Model 1510) on magnetic tape. The data was then analyzed by computer 
using two programs specifically written for the relative naviga- 
tion experiment. 

The first computer program, the Raw Tape Decoder Program, was used 
to read and decode the records , one at a time , and convert the 
binary data to decimal form suitable for further analysis. A 
standard format data tape was generated at this time containing 
a listing of the data for preliminary analysis. In particular, 
the following items were printed out: time, coarse and fine 

range for both receivers, the range difference, and status data. 
Measured range is the sum of coarse range, in hundreds of nano- 
seconds, and fine range, in nanoseconds. Status data recorded 
consisted of search/ track mode, threshold level detected, station 
and channel number, and calibration or track data. 

The second computer program, the Line of Position (LOP) Program, 
operates on the computer tape generated by the Raw Tape Decoder 
program. Initially, the LOP program translates the indicated 
time of each ranging data record, which represents the elapsed 
time since turn-on of the tape recording system, to Greenwich 
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Mean Time. All pre-test and post-test calibration range 
differences are then read and fit to a straight line. Next, 
sets of approximately 20 consecutive range differences of 
satellite data are fit to parabolas. Each such parabola is 
evaluated at its center point in time. The smoothed range differ- 
ence value is reduced by the calibration line interpolated at 
this point in time. The remainder is then used to calculate a 
line of position using the equations described in Section 4. A 
flow diagram of the Line of Position Program is shown in Figure 
3-2. 

3.4 SATELLITE CHARACTERISTICS 

The spinning ATS-5 satellite is the only satellite of the ATS 
series containing an L-band transponder and, thus, provided the 
only vehicle for performing testing in the allocated band for 
marine/aeronautical service. The ATS-5 satellite is rotating 
with its directional antennas scanning the Earth approximately 
once every 780 milliseconds. The 3-dB beamwidth of the satellite 
antenna provides joint illumination of the transmitting and 
receiving sites for approximately 52 milliseconds every revolution. 
This scanning illumination limits the period during which the 
satellite can relay the up-link signal to approximately 52 milli” 
seconds out of every 780-millisecond rotational period. This is 
true for both the L-band and C-band transponders. 

During the test, a malfunction in the reciever of the L-band 
transponder necessitated operating the satellite system in a C- 
band to L-band cross-strap configuration. Normally, the satellite 
is capable of operating in either a C-band to C-band or L-band to 
L-band frequency translation mode with a bandwidth of 25 MHz. The 
cross-strap mode provides reception on either frequency and trans- 
mission on the opposite frequency. However, in the cross-strap 
mode, the bandwidth is limited to approximately 6 MHz. Figure 3-3 
illustrates the actual passband characteristic of the satellite 
when operating in the C-band to L-band cross-strap mode. The 
measurements were conducted by NASA at the Mojave ground station. 
The horizontal axis is the transmitter IF frequency and the 
vertical axis is the output power level normalized to full output 
power in the L-band, wideband frequency-translation mode. 

The ALPHA II down-converter system was designed to receive exactly 
1550.0 MHz, The equivalent Rosman transmitter IF frequency for an 
ATS-5 down-link transmission of 1550 MHz is 69.0495 MHz and 
examination of the passband characteristics of Figure 3-3 indicate 
that operation at this frequency would result in a output power 
approximately 12 dB below normal. Therefore, the Rosman trans- 
mitter was required to operate at an IF frequency of 70.1 MHz and 
the ALPHA II down-converters were offset by slightly more than 1 
MHz using a government-supplied frequency synthesizer in the local 
oscillator injection chains. 
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Figure 3-2. Flow Chart for Line of Position Calculations 
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Figure 3-3. ATS-5 Bandwidth Characteristic in C- to L-band 
Cross -strap Mode 

The ATS-5 satellite is in a near-synchronous Earth orbit with a 
daily in-orbit excursion between +2® in latitude, coupled with a 
slow drift westwards and a geocentric range variation on the 
order of 30 nautical miles. The satellite thus oscillates about 
a point near 106 ®W longitude. The errors in the quoted satellite 
position are governed- by the accuracy to which its range can be 
measured from the NASA Mojave and Rosman tracking stations. 

3.5 EXPERIMENTAL PROCEDURE 

The experiment hardware configuration at NASA Wallops Station 
consisted of two identical receiving systems, a data recording 
system and a duplicate of the modulator installed at the NASA 
Rosman STADAN transmitting site. A simplified functional diagram 
of a single receiving system is shown in Figure 3-4. 

The satellite RF signal received by the antenna is applied to 
the L-band receiver where the signal is down-converted“’^from 1551 
MHz +12.5 MHz to an IF frequency of 70 MHz +12.5 MHz and amplified 
•to a nominal level of 0 dBm. The code correlator demodulates the 
signal, extracting the range information. This range data is 
processed by the tape control unit of the data recording system 
and recorded on magnetic tape. Three identifiable equipment 
biases exist which had to be calibrated in order to assure mean- 
ingful range readings. These are:. (1) the RF front-end bias, 

(2) the PN code generator offset, and (3) the frequency standard 
drifts. 
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Data to Recorder System 


Figure 3-4. Simplified Functional Diagram of 
ALPHA II Receiving System 


The front-end bias encompasses the difference in signal path 
delay from the antenna feed to the input of the code correlator . 
This bias is fixed throughout the experiment provided the physical 
integrity of the equipment is maintained. The RF front-end^ bias 
can be calibrated, as shown in Figure 3-5, by operating both 
receiving systems at the same site with both antennas physically 
located with their feeds adjacent on the same line of position. 

The resulting data is then processed to determine the magnitude 
of the bias. 

The bias produced by the PN code generator offset is a random 
offset due to the positions of the code-generating shift registers 
when power is applied during each test period. This bias is fixed 
and will not change if the integrity of the povrer source and the 
frequency source is not violated. However, the two independent^ 
frequency standards drift with respect to each other as a function 
of time and, therefore, produce a time dependent bias which must 
also be calibrated. The calibration configuration used to deter- 
mine both of these biases is shown in Figure 3-6. 

The internally generated PN codes are derived from the 14-bit 
shift registers in each correlator and the modulator. The cali- 
bration mode enables a determination of the daily phase discrep- 
ancy between each receiving system. As can be seen in Figure 3-5, 
the ALPHA II modulator is used to simulate the transmitted RF 
signal. Since the System ”A" frequency standard is used to pro- 
vide the reference clock to the modulator, the System "A" calibra- 
tion reading drifts in direct proportion to the drift between the 
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Figure 3-5. Calibration Configuration for Frequency 
Standard Drift and Code Offset 


two frequency standards, provided power and clock integrity are 
not violated. To insure this integrity, the System "B" standard 
is operated from its own battery supply and the System "B" 
correlator is powered from batteries and DC to AC solid-state 
inverters . 

At the beginning of each test day the operating condition of 
each receiving system was verified. Once normal operation was 
demonstrated, both systems were configured in the calibration 
mode. After recording the intial system offsets and drift for 
approximately fifteen minutes. System "B" was transported to the 
remote site just prior to satellite transmission. The first few 
minutes of each satellite test period consisted of CW trans- 
mission, which was used to maximize signal strength by manually 
"tuning" the frequency synthesizer to compensate for the ATS-5 
master oscillator offset and drift. The code was then transmitted 
to obtain the precision ranging data. Following the transmission 
period, System "B" was returned to the fixed site and a second 
calibration was performed to determine the total drift of the 
frequency standards. Table 3-2 summarizes the events of each 
scheduled test day. Tests were performed for 14 out of 16 
scheduled test periods for a total satellite use time of 33 hours 
10 minutes. The first six completed test periods were conducted 
in mid-afternoon, the next four just prior to sunrise and the 
remaining tests near 2400 hours local time. 
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Figure 3-6. RF Front End Bias Calibration Configuration 
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TABLE 3-2. SUMMARY OF WALLOPS PRECISION POSITIONING TEST PHASE 


TEST 

NO. 

r— 

OATE/DAY 

SCHEDULED 
TIME (GMT) 

ACTUAL 

TIME 

I^TIEES 

P RE- CAL 

READINGS 

A/B 

5NANOSECONDS) 

POST-CAL 

READINGS 

A/B 

(NANOSECONDS) 

TIME BTWN 
CAL . RDNGS 
(HRS) 

OFFSET 
DF AVG 
LOP 
(FEET) 

♦ OF 
TWTs 

CW 

S/N dB 
U.6KH2) 

ATS-5 

MAST.OSC, 

FREQ.DRFT 

(HZ/HR) 

SPECIAL NOTES 

- 

7-7-71_Wed. 

2200-0100 

CANCELL] 

SD -- NAS 

A PRIORITY - 








• 

1 

7-8-71-Thur 

1630-1930 

1750-2030 

Clear 

8290/202047 

N.A. 



1 

NA 

750 

1. No remote site 

2. Circuit Breaker 
interruption 

2 

7-12-71-Mon 

1630-1930 

1645-1935 

Clear 

127061/ 

45115 

127069/47397 

3.1 

6.9 

1 

NA 

6500 

1. Heaters on for 
one hour before 
test in ATS-5 
2.3-hr rubidium std 
calibration run 

3 

7-13-71-Tues 

1630-1930 

1815-2130 

Clear 

127070/ 

129920 

127070/ 

116323 

B.O 

” 

1 

+18.4 

250 

1. Heaters on in ATS-5 
overnight 

4 

7-l4-71-Wed 

1630-1930 

1700-2030 

Clear 

5399/82943 

181092/85442 

5,0 


1 

+19.8 

<250 

2."B" lost reference 

1.5MH2 line to 
rt ' Correlator 
interrupted 

5 

7-l5-?i-Thur 

1630-1930 

1850-2030 

Clear 

39523/89701 

39522/91342 

2.8 

8.5 

1 

+16.4 

<250 


■ 

7-16-71-Fri 

1630-1930 

1700-1930 

Clear 

140574/ 

100779 

91597/102075 

3.0 

45* 

1 

+16.4 

<250 

1. Conducted front 

end bias calibration 

2. "A" lost reference 
during remote site 
test 

1 

7-20-7X-Tue 

0630-0830 

0730-0830 

Clear 

196369/ 

38330 

80498/128446 

N.A. 

1 

1 

+18.6 

<250 

1. Modulator reference 
lost during test 

2. Power down at 
remote site 

3. Phone coupler 
difficulty. 
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TABLE 3-2. SUMMARY OF WALLOPS PRECISION POSITIONING TEST PHASE (Continued) 


TEST 

NO. 

DATE/ DAY 

SCHEDULED 
TIME (GMT) 

ACTUAL 

TIME 

WEATHER 

PRE-CAL 

READINGS 

A/B 

NANOSECONlJ 

POST -CAL 
READINGS 
A/B 

(NAl^OSECONDS) 

TIME. BTWN 
CAL.RDNGS 
(HRS) 

)FFSET 
3F- AVG 
^OP 
(FT) 

# OF 
TWTs 

•— 

cw 

S/N dB 
(1.6KHZ) 

ATS-5 

MAST.OSC. 

FREQ.DRFT 

(HZ/HR) 

SPECIAL NOTES 

8 

7-21-71-Wed 

0630-0830 

0700-0845 

Clear 

196722/ 

131649 

196722/ 

134074 

3.5 

1.2 

1 

+ 17.6 

4750 

1. Heaters were NOT 
on in advance of 
test in ATS-5 

9 

7-22-71-Thur 

0630-0830 

0700-0845 

Clear 

101265/ 

78496 

101268/ 

23159 

3.25 


1&2 

+18.6 

300 

l."B” ref lost 

10 

7-23-71-Fri 

0630-0830 

0700-0845 

Clear 

149222/ 

203669 

149233/ 

77240 

2.33 

— 

1 

+18.8 

<100 

l."B'* ref lost 

11 

7-26-71-Mon 

1300-1600 

1300-1600 

Cloudy 

195624/ 

74941 

195624/9435 

3.85 


1&2 

N.A. 


L.Sys "A” would 
not maintain 
lock — no 
data 

12 

7-27-71-Tue 

1300-1600 

1300-1600 

Clear 

62029/ 

143253 

62033/145842 

4.1 

4.6 

1&2 

+20.6 

<250 

1.3-hr rubiditan 
std calibration 

13 

7-28-71-Wed 

1300-1600 

1300-1600 

Cloudy 

135300/ 

94875 

135299/ 

97645 

4.5 

4.6 

l£2 

+19.2 

<250 

run 

14 

7**29-71-Thur 

7-30-71-Fri 

1300-1600 

0330-0600 

Rain 

102965/ 

154360 

141760/ 

166589 

3.75 


1 

N.A, 


l.Sys "A" 
reference 
lost 














NOTE: The difference in antenna cable lengths and front end delays =45 feet. 

TOTAL ACTUAL ATS-5 TEST TIME * 33 Hours 16 Min 



The entries of Table 3-*2 are generally self-explanatory. The 
number in the OFFSET OF AVERAGE LOP column is calculated elsewhere 
in this report. The S/N measurement was recorded during the CW 
transmissions for the peak signal strength observed during the 
burst of signal. The SPECIAL NOTES COLUMN indicates that during 
six test periods the power and/or frequency integrity of either 
system A or B or the modulator, which was always the reference, 
was violated. This was due to either operator error or mechanical 
deficiency during the transportation of system B to the remote 
site. 

The range reading was measured by the digital subsystem of the 
receiver code correlator consists of a pulse count which is gener- 
ated by enabling a counter to accumulate standard clock pulses 
when the reference PN code is at a specified phase point and 
stopping the counter when the received PN code is at the same 
phase point. The resultant count accumulated in the counter is 
then representative of the time phase difference or the range 
between the transmitted signal and its point of origin and the 
received signal and its point of reception. This count can be 
con^^erted to any suitable unit of measurement by knov^ing the value 
contained in the counter and the clock rate which was being accumu- 
lated. The lane width of the 10-megabit ALPHA II system is 30 
meters (100 nanoseconds). As described in Section 5, a special 
four-phase clock allows this bit interval to be divided into quarters, 
giving a granularity of 7.5 meters (25 nanoseconds). Range differ- 
ences finer than 7.5 meters are resolved by digitizing the D-C 
voltage appearing at the phase detector output of the correlator 
with an A/D converter. The complete system is capable of quantizing 
range readings down to one nanosecond. A front panel display in- 
dicates the resultant range reading as it is being measured. The 
displayed number is indicative of the range in nanoseconds, which 
is approximately the range reading in feet. The range reading is 
then formatted and buffered to the magnetic tape control unit for 
recording. 

The tape control unit performs a parity check on the data received 
from either code correlator. The tape control unit also generates 
a time signal from the 5-MHz frequency standard which is recorded 
with each data record. The time signal is an elapsed time measure- 
ment from the beginning of the test. The time is expressed in 
minutes, seconds and tenths of seconds and the total accumulation 
overflows at 59 minutes, 59.9 seconds. The data record is not 
written if parity errors are detected in the words received from 
either correlator. 

The first 5 bytes of the 13-byte data record contain the range 
reading of the fixed receiving site; the second 5 bytes contain 
the range reading of the remote receiving site, and the last 3 
bytes are the time code. Following the 13 bytes, an inter-record 
gap is generated which will contain the standard LRC check 
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character. The magnetic tape format for the individual bytes 
described above is shown in Table 3-3. The range readings are 
recorded as 8-4-2-1 positive convention BCD digits. The composite 
"FINE" reading will assume values between 0 and 99 nanoseconds. 

The composite "COARSE" reading will assume values between 0 and 
2,046; this value must be multiplied by 100 to convert to nano- 
seconds. The total reading is then the decimal sum of the "FINE" 
and "COARSE" values. 

3.6 EXPERIMENTAL ERROR ANALYSIS 

The following paragraphs systematically identify and evaluate the 
various error sources which contributed to the ranging inaccuracy 
in the position determination experiment. The error calculations 
are based upon actual measurements or equipment characteristics 
whenever possible. Certain error sources, e.g., delays due to 
the troposphere and ionosphere, have effects which can only be 
estimated. Table 3-4 contains the ATS-5 C-band to L-band link 
budget calculations based upon the actual equipment used, with 
the exception that it does not consider the 6 dB loss of signal 
strength due to the satellite passband constraint. The equivalent 
S/N ratios observed each day in the CW mode at the fixed receiver 
site are tabulated in Table 3-2. 

The important error sources identified during the experiment 
include: (1) receiver noise, (2) equipment nonlinearity, 

(3) quantization, (4) phase alignment of the rubidium standards, 

(5) multipath, (6) tropospheric delay uncertainty, (7) ionospheric 
delay uncertainty, (8) atmospheric and galactic noise, (9) satellite 
position uncertainty, and (10) fixed receiver position uncertainty. 
These, of course, are further subject to Geometric Dilution of 
Precession (GDOP) . 

3.6.1 Signal-to-Noise Ratio 


Measurements on the received CW signal indicated that the signal- 
to-noise ratio (S/N) varied between 16.4 and 20.6 dB in the 1.6- 
kHz IF bandwidth, which is about 3 dB less than the calculations 
of the link budget. The tracking loop filter in which the range 
measurement is made is 50-Hz wide. The S/N of the modulated 
signal in this filter varied between 21,4 and 25.6 dB. The 10 dB 
discrepancy is due to a loss of roughly one dB in the transponder 
output power and about 3 dB due to the satellite antenna rotation. 
The additional 6 dB loss was due to the 6-MHz filter in the C-band 
to L-band cross-strap mode in the satellite which had to pass the 
10 megabit PN code. 
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TABLE 3-3. MAGNETIC TAPE DATA FORMAT 


Byte 

No. 

B 

A 

8 

4 

2 

1 

] 

Y Channel A 

Data 

< 

V Channel B 

1 Data 

J 

\ Time Data 

1 

F/40 

F/80 

Spare 

B3 

Spare 

B2 

Ch.l=l 

Ch.2=0 

1 

2 

F/1 

F/2 

F/4 

F/8 

F/10 

F/20 

3 

C/1000 

C/2000 

C/4000 

C/8000 

C/10000 

CAL = 1 
NORM = 0 

4 

C/40 

C/80 

C/100 

C/200 

c/400 

C/800 

5 

C/1 

C/2 

C/4 

C/8 

C/IO 

C/20 

6 

F/40 

F/80 

Spare 

B3 

Spare 

B2 

Ch.l=l 

Ch,2=0 

0 

7 

F/1 

F/2 

F/4 

F-8 

F/10 

F/20 

8 

C/1000 

C/2000 

C/4000 

C/8000 

C/10000 

CAL = 1 
NORM = 0 


C/40 

C/80 

C/100 

C/200 ■ 

C/400 

C/800 

10 

C/1 

C/2 

c/4 

C/8 

C/10 

C/20 

11 

M/2 

M/4 

M/8 

M/10 

M/20 

M/40 

12 

S/4 

S/8 

S/10 

S/20 

S/40 

M/1 

13 

S/.l 

S/.2 

S/.4 

S/. 8 

S/1 

S/2 


F/X = 

Fine Range/BCD X 

(X = 

= nanoseconds) 

C/Y = 

Coarse Range/BCD Y 

(Y = 

: 100 nanoseconds) 

M/A = 

Minutes/BCD A 

S/B 

= Seconds/BCD B 

S/.C = 

Seconds/BCD , C 

All 

BCD are positive 


convention, 8-4-2-1 code. 


3-17 




8T- 


TABLE 3-4. LINK BUDGET FOR THE ATS-5 SATELLITE, IN THE C-BAND 
TO L-BAND CROSS-STRAP MODE 




Parameter 


Transmit Antenna Gain (dB) 

Transmit Network Loss (dB) 

Receive Antenna Gain (dB) 

Receive Network Loss (dB) 

Net Antenna Gain (dB) 

Transmitter Power Out (dBW) 

ERP (dBW) 

Free Space Loss (dB) 

Atmospheric Loss (dB) 

Polarization Loss (dB) 

S/C Antenna Point Loss (dB) 

Ground Antenna Point Loss (dB) 

Net Propagation Loss (dB) 

Received Carrier Power (dBW) 
Receiver Noise Temp (^K) 

Receive Antenna Noise Temp ('^K) ^ 
Receive System Noise Temp (dB - K) 
Noise Power Density (dBW/Hz) 
Carrier/Noise Power Density (dB-Hz) 
Receiver Bandwidth (dB) 
Carrier/Noise Ratio (dB) 


Rosman to Wallops . 

6214-MHz 

1551- 

MHZ 

Uo-Link 

Down Link \ | 

(85 ft Antenna) 

3 ft Antenna 

4 ft Antenna 

+ 62.0 

+ 14.3 

+ 14.3 

- 1.5 

- 1.3 

- 1.3 

+16.3 

+ 20.0 

+ 22.6 

- 0.4 

- 1.5 

- 1.5 

+76.4 

+ 31.5 

+ 34.1 

+ 36.0 

+ 11.0 

+ 11.0 

+ 96.5 

+ 24.0 

+ 24.0 

-199,9 

-188.0 

-188.0 

- 0.2 

- 0.2 

- 0.2 

- 3.0 

- 0.5 

- 0.5 

- 0.5 

-0.3 to -2.3 

-0.3 to -2.3 

- 0.2 

- 0.5 

- 0.5 

-203.8 

-189.5 to -191.5 

-189.5 to -191,5 

-91.4 

-147 to -149 

-144.4 to -146.4 

890° 

600° 

600° 

290° 

55° 

55° 

30.7 

28.2 

28.2 

-197.9 

-200.6 

-200.6 

106.5 

51.6 to 53.6 

54.2 to 56.2 

67.8 

32 

32 

38.7 

19.6 to 21.6 

22.2 to 24.2 


NOTE: 


see paragraph 3.4 



3.6,2 Internal Errors 


Receiver Noise' r i 


The correlation function "Z” curve represents the error between 
received and locally generated codes. This curve is shown in 
Figure. 3-7. ‘ . ... 

As derived in Appendix D, 

50 

/ S/N 

The resulting range error corresponding to a signal-to-noise 
level of 23 dB is then 

3.5 nsec 

3.6 ft. 

Since two recei vers were used, the resulting noise error increases 
by a factor of , thus ^ 5.1ft. 

Correlation Curve Nonlinearity 


The correlation curve of Figure 3-7 in reality exhibited certain 
nonlinear characteristics, as opposed to the idealized character- 
istics shown. This resulted in a prevalence of certain fine range 
readings (between zero and 99 nanoseconds) and a relative shortage 
of other fine range readings. This characteristic was noticeable 
near the terminal points of the analog-to-digital converter inter- 
vals, i.e., (0 and 24), (25 and 49), (50 and 74) and (75 and 99). 

The values 0 and 24 were recorded much more often than the values 
21, 22 and 26. The same situation was observed at 49, 74 and 99 
nanoseconds. The error due to this behavior is difficult to 
determine theoretically. However, laboratory tests indicate that 
a 2% nonlinearity of the "Z" curve is reasonable. Therefore, 
0^-2 feet for one receiver or 2.8 feet for two receivers, as 
used in the test. 
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Quantization Error 


Quantization error refers to the inability to determine the exact 
location of the correlation maximum which occurs somewhere within 
the smallest resolvable interval. This interval is 1-nsec wide 
and represents the nearest value to which range is resolved. 

Since the correlation maximum can occur anywhere within the +.5 
nsec interval around the measured fine range value, then 


o = — • 

^ ,/T2 

= .29 (ft) 

where q = 1 ns. 

The / Vz' term results from the assumed uniform distribution over 
the 1-nsec interv al . Since two receivers are used, the error 
increases by the / 2 , thus 

a = .42 feet 

q 

Phase Alignment Error 


In the test configuration, the two rubidium standards drift apart 
in phase as well as in frequency. The noise on the rubidium 
standards was estimated on the basis of the typical stability of 
these standards, 5 parts in 10^0 per day. The resulting drift 
per day would then be 5 x 10“i0 x 8.64 x lO^-^ (nsec/day) or 
4.32 X 10"^ nsec/day. The corresponding drift per minute would 
be 30 nanoseconds, a portion of which is clearly correctable. 

Results of 3-hour noise-free calibration runs on July 12 and 
July 27, yield rms values, of 15.9 and 13.2 feet, respectively. 
These rms values result from the deviations observed about the 
least squares parabolic fit. However, the rms values are con- 
taminated by the error resulting from the nonlinear correlator. 
When 2.8 feet is removed from the observed rms values to account 
for this nonlinearity, the remaining rms due to the two rubidium 
standards is 15.7 and 12.9 feet for the July 12 and July 27 tests, 
respectively. The average of the two test periods is 14.3 feet, 
which is in reasonable agreement with published test data (see 
Figure 3-8 taken from Reference 3-5) . From Figure 3-8 it can be 
seen that a single frequency standard would yield an error of 
about 10 feet if the nearest calibration point were an ho ur and 
a half away. Thus, with two standards, an error of / 2 ' x 10 
feet, or about 14 feet, would be expected. Therefore, although 
the error due to the rubidium standards is clearly dependent on 
the time from the test to the closest calibration point, it was 
(jgcided to account for this with a non— time— varying term. 
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Figure 3-8. Barnes' Variate Difference Data and 
GTC 304B Rubidium Oscillator Test 


3-21 



3.6.3 Environmental Errors 


Multipath 

Error due to multipath is caused by the simultaneous reception 
of the signal over more than one path. One is the direct path, 
while the others result from reflections. The multipath at 
Wallops Station is stationary since the reflections are caused 
by land. 

The rms on the range measurement due to multipath is given by 
(Reference 3-1) 


_ p h sinE 


m 


/G 


SR 


where p is the surface relfectivity (30%) , h is the height of the 
receiving antenna (20 feet for the fixed receiver, 3 feet for the 
remote receiver), E is the elevation angle (36 degrees) and Gsr 
is the power ratio of main lobe response to side lobe response in 
the direction of the reflected ray (20 dB) . Thus, is .36 feet 
for the fixed site, and .05 feet for remote site. In addition, 
the characteristics of the ranging code itself have a significant 
effect on reduction of the multipath rejection. However, this 
effect will not be considered in this analysis. 


Trop ospheric Delay 

The principal measurement error introduced by the ionosphere and 
troposphere is an increase in the effective path length between 
the satellite and ground station over the corresponding path 
length in free space. Part of this refractive delay is predictable, 
assuming a priori knowledge of the geometry between the ground 
station and the satellite, time of day and year, sun-spot activity, 
and surface conditions. A portion of the effect will be unpre- 
dictable, varying with anomalies in the ionosphere and troposphere. 
The following estimates are a condensation of many references 
which have been reviewed to arrive at a reasonable estimate of 
both ionospheric and tropospheric delay uncertainties. No attempt 
has been made to model the tropsopheric delay and, hence, reduce 
the uncertainty introduced in the range measurements. 

The incremental increase in path length due to tropospheric 
effects is, in comparison with ionospheric effects, quite pre- 
dictable and independent of carrier frequency in the range of 
100 MHz to 10 GHz. The refractivity of the troposphere is given 

by 

N = 77.6 + 3.73 X 10^ 


3-22 



where P is total atmospheric pressure in millibars; e is water 
vapor pressure in millibars; and T is temperature in degrees 
Kelvin. 

The Central Radio Propagation Laboratory of the National Bureau 
of Standards developed an Exponential Reference Atmosphere 
model and adopted this model for prediction of refraction 
phenomena in the troposphere. The results of this model are 
plotted in Figure 3-9 for two extremes of surface refractivity 
(Ns) versus elevation angle. 

The incremental increase in path length due to downlink tropo- 
spheric delay, at an elevation angle (E) of 36®, is 

ALt = 14 feet 


The rms variation on this path length is typically taken as 10% 
of path delay. Hence 


a 


ALt 


1.4 feet 


Both ionospheric and tropospheric inhomogeneities tend to be 
correlated over small distances with the result that nearly 
parallel paths tend to suffer equal range-error fluctuations 
which tend, therefore, to cancel in terms of range difference. 



Figure 3-9. Tropospheric Effects; Incremental Path 
Length vs Elevation Angle 
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Figures 3-10 and 3-11 indicate the extent of this correlation or 
range-difference error reduction factor as a function of the 
range separation. These data were recorded as a part of the 
National Bureau of Standards series of Maui experiments (Reference 
3-7) . Figure 3-10 is a plot of interferometer angular error (Atj>) 
due to tropospheric fluctuations. The data is taken in a bandwidth 
of the order of 10 Hz, but tends to be dominated by very low fre- 
quencies below approximately .002 Hz for baselines of interest. 

For the normal' baseline , the angular fluctuations (o(|)) is directly 
related to range difference fluctuation (Ofy;)) through the baseline 
B. 


For baseline length of 8 kilometers, a<j> 1.5 pirad so 
Oj^D - .012 meters 

Figure 3-11 is the range difference directly taken from a similar 
experiment and indicates for a baseline length of 8 kilometers, 
^RD ~ meters. 

The path for these measurements is that from Haleabrela beach on 
Maui to the island of Lanai and is intentionally chosen as one 
highly subject to cloud and high wind turbulence, probably among 
the most inhomogeneous type of air that could be encountered. 
Similar measurements in Colorado have shown almost an order of 
magnitude smaller fluctuations. 

Thus we conclude that the errors due to the troposphere are 
negligible in this experiment. 

Ionospheric Delay 

The incremental path-length difference (ALj) between effective 
path length and free-space path length, due to ionospheric 
refraction is 


40.25 N Q(E) 

ALj = (feet) 

f2 

where Nrp is the total electron content (number/m^) ; Q(E) is the 
obliquity factor, a function of elevation angle (E) to the 
satellite; and f is the carrier frequency. 

The total electron content (Nij) varies about its approximate mean 
of 10l*^/m2 by at least 30% in the course of a day and also from 
day to day. Figure 3-12 presents a plot of the variation in 
incremental path length versus elevation angle for a nominal 
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value of Nt extreme Nij values of 0.7 x 10^^ and 1.3 x 10 

The principal variation of Nt in the ionosphere is due to tne 
sun spot number, time of day, season and latitude. In addition, 
there are day-to-day fluctuations, simultaneous spatial variations 
due to ionospheric irregularities, and short-term (order or 
minutes) variations due to motion of the irregularities. 

'The uncertainty in Nt at a particular time and placesis estimated 
to be approximately 30%. This uncertainty includes a +20% day- 
to-day variation and the uncertainty in smoothed sun spot nu^er. 
The minimum, nominal and maximum path lengths errors on the down- 
link can be obtained from Figure 3-12, 

ALj = .5 ft. minimum 

= 7 ft. nominal 


= 100 ft. maximum 

In addition, at any one time there will be a random spatial 
variation of Nm due to the presence of irregularities in the 
ionosphere, which is dependent on the distance between the points 
through which the signal passes, as indicated in Figure 
For this experiment, the fluctuation can be seen to be = 0.1%, 
a negligible value. 



Figure 3-12. Ionospheric Effects at 1600 MHz 
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Figure 3-13. Short Distance Fluctuations in Electron Content 


Atmospheric and Galactic Noise 


In addition to receiver noise, noise also originates in the 
atmosphere from lightning discharges, as well as from electrical 
machinery (see Reference 3-3) . However, since receiver noise is 
usually the dominant noise source at frequencies above 300 MHz, 
atmospheric and galactic noise are neglected in this analysis. 

3.6.4 Geometric Errors 

Satellite Position Uncertainty 

Errors in satellite position are governed by the accuracy to which 
the ranges are measured from the Mojave, Calibornia, and Rosman, 
North Carolina, ground stations. The present NASA tracking 
procedure can estimate the position of the satellite within a 
10 km X 10 km X 10 km box, centered at the predicted point. The 
motion of the satellite relative to the predicted orbit is not 
known . 

The 10 X 10 X 10 km box would produce an error vector of 8.7 km 
from the predicted position to the farthest corner of the error 
cube. The angle subtended at the satellite by the 4.7-nmi distance 



between the receiver sites is 2.4 x 10*''^ radians. Hence, the 
maximum error due to satellite position is 7.13 feet {8.7 km x 
3281 ft/km X 2.4 x 10“4 rad). 

Position Error of the Fixed Receiver 


The range measurements from the fixed and movable receiver to the 
ATS-5 satellite, denoted by and Mg, respectively, are ambiguous 
measurements of the true ranges and D^. In order to determine 
the location of the movable receiver with respect to the fixed 
one, M^, Mg and are used, with calculated from the geometry. 
However, the location of the fixed receiver must be accurately 
known, as well as the location of the satellite. Errors due to 
uncertainty in satellite position have already been discussed. 
Since the receivers are located in areas for which type I surveys 
are available, it is assumed that no new errors are introduced 
into the relative range measurements by the position of the 
fixed receiver. 

Geometric Dilution of Precision (GDOP) 


Geometric dilution of precision refers to the attenuation or 
amplification of range errors. These range errors contribute 
to horizontal errors in the computation of receiver position to 
an extent determined by the relative geometry of the receivers 
and satellite. GDOP thus is not an error source itself, but 
represents an error multiplication. 

The magntidue of the GDOP is given by 1/cosE, where E is the 
elevation angle of the satellite. For Wallops Station, the 
elevation of the ATS-5 satellite is 36 degrees. Hence, the 
GDOP is 1/cosE, or 1.25. 


If is the rms range error due to all contributing error 
sources, the rms position error is given by 


pos 


= O- X 1.25 
1 \ 


3.6.5 Position Error Summary 


The error sources and the resulting uncertainties are presented 
in Table 3-5, 


Range accuracy is determined as the square root of the sum of the 
squares of the individual range errors. The position accuracy is 
then found by multiplying the range accuracy by the GDOP factor. 
Thus, 


a - 21.3 feet 

pos 
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TABLE 3-"5. RELATIVE RANGING ERRORS AT WALLOPS STATION 


Error Source 

Relative Range 
Uncertainty 
(ft) 

Receiver Noise 

5.1 

Nonlinear Correlator 

2.8 

Quantization 

0.42 

Rubidium Frequency Standard 

14.3 

Multipath 

0.4 

Troposphere 

Negligible 

Ionosphere 

Negligible 

Satellite Position Uncertainty* 

7.0 

Survey Error 

Negligible 

Range Accuracy (RSS) 

17.0 

*This would be the result of a 10 x 10 x 10 km error which 
we believe to be high. 
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SECTION 4 


EXPERIMENTAL DATA AND ANALYSIS 


The experimental test data was measured by the ALPHA II receiving 
systems in several test configurations, as described in the 
previous section. , The calculation of the position of the remote 
receiver relative to the fixed known site necessitated calibration 
of several equipment characteristics. The results of the calibra- 
tion tests, the satellite ranging tests and the resulting calcu- 
lated lines of position describing the location of the remote 
receiver are presented in this section. An evaluation of the 
contribution of several external error sources is also included to 
explain the data reduction method. 

4.1 SUMMARY OF THE TEST RESULTS 

A summary of the complete test cycle is given in Table 4-1, showing 
valid LOP test data measured during six test periods. During five 
of these test periods the B receiving system was located at the 
remote site and data suitable for relative position determination 
was collected during those periods. 

The sixth fully successful test period was used to calibrate the 
receiver front-end delay bias by co-locating the two systems with 
the antennas mounted on the same line. of position. The average 
LOP calculated for this test was 56 feet removed from the fixed 
site. Thus, the delay bias between the two systems is 45 feet. 

The average calculated LOP’s for the five successful test periods 
with receiver B located at the remote site are summarized in 
Figure 4-1. A statistical summary of these tests is given in 
Table 4-2. 

The offsets of the average LOP's in Figure 4-1 were measured by 
using the scaling factor of one inch equal to 18.23 feet. The 
standard deviation (a) of the LOP's about the average LOP for 
each test is also given in Table 4-2. 

The estimate of (o) is obtained from the detailed LOP plots in 
Figures 4-6 to 4-10 as follows. The number of individual LOP's 
to one side of the average LOP are counted and the LOP corres- 
ponding to 2/3 of this total is singled out. The same procedure 
is repeated for the LOP's on the opposite side of the average 
and the distance between the two selected LOP's is measured. This 
distance is assumed to be 2a. 
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TABLE 4-'l. SUMMARY OF EXPERIMENTAL TEST CYCLE 


July 

1971 

7 

8 


Test 

Number 


1 


12 

13 

14 

15 

16 
20 


2 

3 

4 

5 

6 
7 


21 

8 

22 

9 

23 

10 

26 

11 

27 

12 

28 

13 

29 

14 

30 



Test Disposition 

Cancelled - NASA priority 

Circuit breaker interrupt at fixed 
site - co-located test 

Valid data 

Receiver B lost reference 
Receiver A lost reference 
Valid data 

Front end calibration, co-located test 

Battery of rubidium standard B low; 
Power failure at remote site, no data 
link. 

Valid data 

Receiver B lost reference 
Receiver A lost reference 
Equipment problems with receiver A 
Valid data 
Valid data 

Receiver B lost reference 
Cancelled 







Figure 4-1. Average 



's for Test Periods 2, 5, 8, 12 and 13 








TABLE 4“2. SUMMARY OF TEST DATA 


Test 

Date 

Offset of 
Average LOP 
(ft) 

Rras 

(ft) 

GMT 

2 

July 12 

6.9 

25.88 

1856-1930 

5 

July 15 

8.5 

31.63 

1702-1734 

8 

July 21 

1.2 

16.77 

0759-0845 

12 

July 27 

4.6 

33.26 

0131-0403 

13 

July 28 

4.6 

30.52 

0200-0402 


4.2 CALIBRATION OF RF FRONT END BIAS 

On July 16/ both complete receiver systems, including antennas, 
were co-located at the fixed site in order to determine the rela- 
tive delay due to the front ends. This test was important since 
any delay resulting from these systems, if uncorrected, would affect 
the determination of the lines of position, resulting in a bias on 
the LOP ' s . 

The antennas were placed on the roof of Building F-160 and physi- 
cally aligned along the same perpendicular (LOP) to the line of 
sight to the satellite and matched in azimuth and elevation. A 
normal satellite test period was then performed. 

The pre- and post-calibration periods were fit to a straight line 
in the least squares sense. The average difference in feet between 
the data taken on the satellite and the calibration line was then 
determined. This average was 45 feet, with a spread from 30 feet 
to 61 feet. The corresponding lines of position in Figure 4-2 were 
obtained by projecting the observed differences just defined through 
the elevation angle along the line of sight to the satellite. Since 
*:.the position predictions of the ATS- 5 satellite were not used in 
'■■'■determining the receiver front-end delay bias, inaccuracies associ- 
• '■■ated with the ephemeris are eliminated. 

4.3 CALIBRATION OF THE FREQUENCY STANDARD DRIFT 


Proper modeling of the relative drift in frequency and phase between 
the two rubidium standards is the most important step in the deter- 
mination of the lines of position. Prior to the tests, it was 
assumed that the relative drift in frequency was significant. Since 
measurement of distance was sought, it was decided to model the 
relative behavior of the rubidium standards by a least squares 


4-4 





s RF Front End Bias Calibration 














parabola (At2 + Bt + C) , where the quadratic term in time repre- 
sents the frequency drift, the linear term in time the phase 
drift, and the constant represents the initial offset between the 
range measurements of the two receivers. In order to test this 
hypothesis, two special 3-hour calibration test periods were 
conducted on July 12 and July 27, respectively. The parabolas 
were based on 4789 and 5000 contiguous points, respectively. The 
resulting coefficients calculated from the test results are pre- 
sented in Table 4-3. 

These test results show that the quadratic coefficients were very 
small compared to the linear coefficient, being 4.5 and -6.3 feet 
in three hours, versus -1644.9 and -1952.4 feet. It is, therefore, 
concluded that the rubidium standards did not drift relative to 
each other in frequency. 

The same model was also used to fit the calibration data of the 
various experimental test periods. A typical test resulted in 
the following parabola. 

P2 (t) = 10.3t2 ~ 732. 8t + 93217.0 

(Test 2 based on 712 points) 

On both July 12 and July 27, a normal satellite test period was 
executed in addition to the special 3-hour calibration runs. The 
quadratic coefficient observed in the satellite test was two 
orders of magnitude larger than for the 3-hour calibration runs. 

To understand the cause of this conflicting result and obtain a 
better model, the 3-hour calibration run of July 27 was subjected 
for further analysis. To simulate a satellite test sequence, the 
first 38 data points of the calibration run were entered as time 
and observed range difference, and 94 points, beginning almost 3— 
hours later, were entered into the curve- fitting routines, thus 
simulating the discontinuity that occurs in an actual satellite 


TABLE 4-3. COEFFICIENTS OF LEAST SQUARES FIT PARABOLA* 
FROM SPECIAL 3-HOUR CALIBRATION TESTS OF 
RUBIDIUM STANDARDS 


Date 

Number of 
Points 

A 

B 

C 

Rms 

(ft) 

July 12 
July 27 

4789 

5000 

.5 

-.7 

-548.3 

-650.8 

78114.0 

62877.0 

15.9 

13.2 

*At2 + Bt + C, where t is in hours 
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test period. The total number of points entered (132) was typical 
for an average satellite test period. The time interval between 
sets was smaller than for a typical satellite test period. The 
least squares parabola and least squares straight line were then 
calculated based upon this set of 132 points. The continuous data 
resulting from the three-hour test was also entered. All three 
functions, namely, 

a the least squares line - L(t) = 62,876.6 - 650. 92t 

o the least squares parabola - P(t) = 62,878.5 - 733. 31t + 

27.48t2 

9 and the 3-hour parabola - 0(t) = 62,877.0 - 650. 8t - 

. 7t2 

were evaluated for a continuous set of points for which difference 
readings were available. The selected set was near the center of 
the time interval. The results are given in Table 4-4. Time and 
difference readings are given in Columns 1 and 2. The parabolic 
fits, namely the ones based on 5000 and 132 points are given in 
Columns 3, 4, 5 and 6. The deltas from the observed range differ- 
ences are also given. The last two columns represent the straight 
line fit. 

The P(t) parabola differs from the observed range differences by 
40 to 70 feet but the straight line fit and the original parabola 
based upon the continuous data set do not differ from the observed 
range difference readings by more than about 15 feet. Thus, it 
was concluded that the parabola is not a good fit to drift char- 
acteristics of this standard. 

The data points collected in the pre- and post-calibration periods 
were then fit to a straight-line model of the frequency standard 
drift. The individual lines are summarized in Table 4-5. 

The rms values about the moving mean are also shown in Table 4-5. 
The rms noise for Test 8 is large. The pre- and post-calibration 
periods are shown in Figures 4-3 and 4-4, respectively. The pre- 
calibation is very noisy due to an improper threshold setting in 
code correlator B. The rms for Test 13 is also larger than 
expected; however, in this case the range difference did not 
decrease at a near-constant rate. The range difference remained 
nearly constant for periods as long as two minutes in some 
instances, while during other 2-minute time spans the range 
difference decreased at nearly 1000 nsec/hour. 

Figure 4-5 shows the pre-calibration data for Test 2. The cali- 
bration line is also shown. This figure is typical of the cali- 
bration data observed during the other satellite test periods. 
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TABLE 4-4. LINEAR AND PARABOLIC LEAST SQUARES FITS TO TEST DATA FROM THE 
THREE-HOUR JULY 27 CALIBRATION RUN 


Time 

hr min sec 



Parabolic 

Fits 

1 

Linear 

111 4 

Difference 

Rnnn points 1 

132 Points 1 

. 132 Points I 

Keauing 
8 ft 

Value 

Observed 

-0(t) 

Value 

Observed 
-P (t) 

1 III" ”1 

Value 1 

Observed 

-L(t) 

1 27 48.4 

61920,7 

61923.1 


61864.2 

56.54 

61924.0 

-3.26 

1 27 49.5 

61948.2 

61922.9 

25.33 

61864,0 

84 . 24 

61923.8 

24.44 

1 27 50.0 

61944.3 

61922.8 

21.52 

61863.9 

80.43 

61923,7 

20,63 

1 27 50.5 

61941.4 

61922.7 

,18.71 

61863.8 

77.62 

61923.6 

17.82 

1 27 51.1 

61919.7 

61922.6 

-2.88 

61863.7 

56.03 

61923.5 

-3.77 

1 27 52.6 

61924.6 

61922.3 

2.29 

61863.4 

61.20 

61923.2 

1,40 

1 27 53.2 

61931.5 

61922.2 

9.30 

61863.3 

68.21 

61923.1 

8.40 

1 27 54.2 

61986.6 

61922.0 

64.58 

61863.1 

123.47 

61922.9 

63.69 

1 27 55.3 

61929.6 

61921.8 

7.78 

61862.9 

66.69 

61922.7 

6.88 

1 27 56.3 

61935.5 

61921.7 

13.86 

61862.7 

72.77 

61922.6 

12.97 

1 27 57.4 

61934.5 

61921.5 

13.06 

61862.5 

71.97 

61922.4 

12.16 

1 27 59.5 

61931.5 

61921.1 

10.44 

61862.2 

69.25 

61922.0 

9.54 

1 28 00.0 

61934.5 

61921.0 

13.53 

61862.1 

72.44 

61921.7 

12.63 

1 28 01.6 

61937.4 

61920.7 

16.72 

61861.8 

75.63 

61921.6 

15.82 

1 28 02.1 

61927.5 

61920.6 

6.91 

61861.7 

65.92 

61921.5 

6.11 

1 28 03.7 

61918.7 

61920.3 

-1.60 

61861.4 

57.31 

61921.2 

-2.50 

1 28 04.2 

61921.7 

61920.2 

1.49 

61861.3 

60.40 

61921.1 

. 59 

1 28 04.7 

61905.0 

61920.1 

-15.12 

61861.2 

46.80 

61921.0 

-16.02 

1 28 06.8 

61904.0 

61919.8 

-15.73 

61860.8 

43.18 

61920.7 

-16.64 

1 28 08.9 

61898.1 

61919.4 

-21.25 

61860.5 

37.66 

61920,3 

-22.16 

1 28 09.4 

61919.7 

61919.3 

.44 

61860.4 

59.35 

61920.2 

- .47 

1 1 28 10.5 

61917.8 

61919.1 

-1,26 

61860.2 

57.65 

61920.0 

-2.17 
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TABLE 4-4. LINEAR AND PARABOLIC LEAST SQUARES FITS TO TEST DATA FROM THE 
THREE-HOUR JULY 27 CALIBRATION RUN (Continued) 


Parabolic Fits | Linear Fit 


hr 

Time 
min sec 

Difference 

Reading 

5000 Points 1 

132 Points 

134 Points 

.8 ft 

Value 

"Observed 

-0(t) 

Value 1 

Observed 

-P(t) 

Value 

Observed 

-L(t) 

1 


13,6 

61900.0 

61918.5 

^18.50 

61859.6 

40.41 

61919.4 

-19.41 

1 

28 

14.1 

61916.8 

61918.4 

-1.61 

61859.5 

57.30 

61919.3 

-2.52 

1 

28 

15.7 

61910.9 

1 61918.1 

-7.22 

61859.2 

51.69 

61919.0 

-8.13 

1 

28 

16.2 

61928.6 

61918.0 

10.57 

61859.1 

69.48 

61919.0 

9.66 

1 

28 

16.8 

61922.7 

61917.9 

4.78 

61859.0 

63.69 

61918.8 

3.87 

1 

28 

17.8 

61898.1 

61917.8 

-19.64 

61858.9 

39.27 

61918.7 

-20.55 

1 

28 

18.3 

61904.0 

61917.7 

-13,65 

61858.8 

45.26 

61918.6 

-14.56 

1 

28 

18.9 

61917.8 

61917.6 

.26 

61858.7 

59.17 

61918,5 

- .56 

1 

28 

19.9 

61917.8 

61917.4 

.44 

61858.5 

59.35 

61918.3 

- .47 

1 

28 

20.4 

61897.1 

61917.3 

-20,17 

61858.4 

38.74 

61918 . 2 

-21.08 

1 

23 

21.0 

61911'.8 

61917.2 

-5.36 

61858.3 

53.55 

61918.1 

-6.27 

1 

28 

21.5 

61915.8 

61917.1 

-1.27 

61858.2 

57.64 

61918,0 

-2.18 

1 

28 

23.1 

61910.9 

61916.8 

-5.88 

61857.9 

53.03 

61917.7 

-6.79 

1 

28 

23.6 

61915.8 

61916.7 

- .89 

61857.8 

58.02 

61917.6 

-1.80 

1 

28 

25.2 

61913.8 

61916.4 

-2.60 

61857.5 

56.31 

61917.3 

-3.51 

1 

28 

25.7 

61920.7 

61916.3 

4.39 

61857.4 

63,30 

61917.2 

3.48 


i 




TABLE 4-5. STRAIGHT LINE FIT* TO THE CALIBRATION DATA 


Test 

Date 

A 

B 

Rms (ft) 

2 

July 12 

-547.54 

89660.22 

14.02 

5 

July 15 

-618.84 

-38074.84 

18.43 

8 

July 21 

-673.28 

68018.83 

75.13 

12 

July 27 

-605.74 

-79019.48 

10.45 

13 

July 28 

-638.78 

39765.47 

10.25 

★At + 

B (t in hours) 




4.4 CALCULATED LINES OF POSITION 

The data for the satellite test periods are tabulated in Appendix 
A and the resultant lines of position are plotted in Figures 4-6 
through 4-10. These figures are Mercator projections; thus, any 
measurement of length made on these figures is directly related 
to change in latitude. The scale used is one inch in latitude 
(.000050 degree) equals 18.23 feet. 

In each figure, the remote site is indicated at latitude 37.864343 
degrees and longitude 75.507593 degrees. The two other latitudes 
shown are the assumed reference latitudes of 37.864300 and 37.864253 
degrees used in the calculations. 

The longitudes of the LOP*s were calculated for each of these 
latitudes. The average LOP for the complete test period is 
indicated as the dashed line. This average LOP results from 
averaging the sets of longitudes pertaining to each reference 
latitude. The azimuth heading from north for all LOP’s is 
approximately 132 degrees. The line perpendicular to the LOP*s 
is the line of sight to the satellite. 

The lop's presented in Figures 4-7 through 4-10 are based on 
approximately twenty successive range differences which were fit 
to a parabola. Thus, the noise on the range differences has been 
smoothed to some extent. Any reference to time has been lost in 
these figures in that they show latitude versus longitude. Since 
the headings of the LOP's are all approximately 132 degrees from 
true north, any LOP is equally well defined by the longitude at, 
for example, the reference latitude of 37.864253 degrees. The 
longitude at the other reference latitude of 37.864300 degrees 
exceeds the first longitude by .000063 degree. It is then possible 
to plot longitude versus time as is done in Figures 4-11 through 
4-15 for Tests 2, 5, 8, 12 and 13 for the first 30 minutes of each 
test. The indicated longitudes are still the result of smoothing 
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Figure 4-4. Test 8 Post-Calibration Data 

Reference Time : 0 ~ 9:25:54.3 GMT 














Figure 4~5. Test 2 Pre-Calibration Data 

Reference Time: 0 = 16: 31: 37 ,9 GMT 

a group of successive measured range differences (AR) . Figures 
-4-16 and 4-17 are also plots of longitude versus time. However/ 
each longitude indicated is the result from one range difference. 
Hence, no smoothing was done on this data, A set of successive 
AR's were analyzed for Test 2 (Figure 4-16) and two sets of 
successive AR’s for Test 12 (Figure 4-17). Twenty successive 
longitudes in these figures describe the structure within one LOP 
or longitude in Figures 4-11 and 4—15. The data for all these 
figures is presented in Appendix A. 

Consider now a configuration in which a user wishes to determine 
his position (LOP’s relative to one or two satellites) in a short 
period of time. In particular, we are addressing relative navigation 
for a ship. If Or represents the position error on the surface of 
the Earth, while 03 represents the permissible error in the LCP’s 
and if 

1 ^R 

Os < Or 


The reduction in error as a function of integration time or number 
of samples is schematically depicted in Figure 4-18. 


then one measurement suffices; or if 
o2 


then N = 


R 


measurements are required. 
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Thus, initially the residual error decreases as — when N samples 

are integrated. However, if integration is increased above 
certain limit, no further reduction in o can be obtained. This 
constant level is due to inherent noise sources in the system, 
e.g., short term noise due to the frequency standard. 

The results from the Wallops tests indicate that a can be reduced 
to about 30 feet if 20 data points are smoothed. The time frorn^ 
measurement to measurement is 0.8 second due to satellite rotation 
while each measurement lasts 20 microseconds. The 16 seconds could 
be considerably reduced without loss in accuracy if the navigation 
signals originated from non-spinning satellites. 

4.5 LOP ERRORS FROM EPHEMERIS PREDICTS 

In order to determine a line of position, the position of the 
ATS-5 satellite at that instance of time must be known. NASA 
Gphemeris predictions for tlie ATS-5 satellite yield the geocentric 
range, latitude and longitude at 30-minute intervals. Typically, 
five or six consecutive ephemeris points are required over a test 
period. The sets of latitude, longitude and range points are fit 
to least squares parabolas. These parabolas are then evaluated at 
the points in time corresponding to the ranging test data. 



Figure 4-11. West Longitude of Test 2 LOP ' s at 
Latitude 37,864253 Degrees 
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Figure 4-12. V7est Longitude of Test 5 LOP ' s 
at Latitude 37.864253 Degrees 



Figure 4-13. West Longitude of Test 8 LOP ' s 
at Latitude 37.864253 Degrees 
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Figure 4-14. West Longitude of Test 12 LOP's 
at Latitude 37.864253 Degrees 



Figure 4-15. West Longitude of Test 13 LOP's 
at Latitude 37.864253 Degrees 
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Figure 4-16. West Longitude Crossing Points at Latitude 
37 .864253 Degrees for Test 2 











Figure 4-18. Residual Error Following 
Smoothing of N Samples 


Two possible errors result from this approach. The first is due 
to the inaccuracies associated with the predicted satellite 
positions; the second is due to the selection of the parabolic 
model of the orbit, which is insignificant, as shown below. 

According to NASA (see reference 3-5) , the position of the satellite 
is true within a cube of 10 x 10 x 10 km centered at the predicted 
point. This inaccuracy in epheraeris prediction is sufficient to 
effect the relative position determination. Thus, it was neces- 
sary that the latitude and longitude of the fixed site be known 
to six decimal places with the Type I survey. 

To determine the effect of a satellite position error within the 
confines of the ephemeris uncertainty cube, the predicted satellite 
position on July 21 at 2100 GMT was selected and the position per- 
turbed in range, latitude and longitude individually. The range 
was perturbed by +5 km and the latitude and longitude were per- 
turbed by +.016 degrees (5 km at 22,000 nmi) . The results are 
summarized in Table 4-6. 

The degree to which LOP longitude changed is a nearly linear 
function of the variation in ephemeris. Therefore, individual 
component errors can be root sum squared. It should be noted 
that the effect of change in range is very small compared with the 
effects of change in angle. It is not known in which manner the 
true position of the satellite differs from the predicted position, 
but the deviation was probably periodic. If so, over a short 


4-23 



TABLE 4-6. LOP ERROR AS A FUNCTION OF SATELLITE ' 
POSITION ERRORS 


Parameter 

Change in Longitude 

Deg 

Feet 

Range +5 km 

+0.000002 

+ 0.4 

Latitude +0.016 deg 

+0.000023 

+ 4.6 

Longitude +0.016 deg 

+0.000059 

+11.4 


period of time, the deviation would be consistently in one direction, 
introducing a bias in the LOP's. This may be the reason the 
average LOP of any given data set is offset by several feet from 
the remote site. 

The second error source was due to the parabolic fit to the 
satellite's latitude, longitude and geocentric range values. A^ 
parabolic fit was selected to approximate the nearly closed ellip- 
tical path of the ATS-5. The rms values of the data points about 
the parabolas are given in Appendix A. These values are generally 
two orders of magnitude smaller than the uncertainty in predicted 
values, indicating very good fits. The rms on the longitude, the 
major sources of LOP error, was on the order of 0.00020 degree. 

The interpolated satellite longitude would, at worst, introduce 
a 3 a error (0.00060 degree), which corresponds to 1/400 of the 
errors exhibited in Table 4-6 of 0.029 foot. Therefore, errors 
resulting from the parabolic fits to the ephemeris data are 
negligible . 

4.6 RUBIDIUM STANDARD NOISE EVALUATION 

The 3-hour calibration runs of July 12 and July 27 indicate that 
the relative drift in frequency of the standards is negligible. 

These data points were fit to least squares parabolas and contin- 
uous data comprising roughly 5000 data points in each run were 
available. The relative drift in phase of the standards was 
significant. When the observed phase drift per hour (^600 nsec) 
was extrapolated to a per day drift and then compared with the 
total number of nanoseconds in a day (8.64 x 10^^), the relative 
phase stability was found to be (600 x 24) (8.64 x 10^^) or 1.5 

in 10^® per day. 

The pre- and post-calibration runs were always conducted in a 
noise-free environment. Thus, any noise on the data was due to 
noise generated by the frequency standards and internally by the 
correlator non-linearities and logic noise. The rms noise ob- 
served on the July 12 data was 15.9 feet versus 13.2 feet for 
July 27 data. The rms noise observed on the range readings of 
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code correlator A were typically 3.5 feet. If the range readings 
of the two correlators are assumed to be independent (which they 
may not be, see below), rms noise values of 15.5 and 12.5 ft. for 
July 12 and July 27, respectively, result from the range readings 
for code correlator B. Removing 3,5 feet of noise for this 
correlator leaves a noise-component magnitude of 15 feet {July 
12) and 12 feet (July 27). This noise then, is due to the 
instability of the rubidium standards. From the error budget 
given in Section 3, it is concluded that the noise due to the 
rubidium standards is a dominant error source. 

Figure 4-19 shows the relative, unsmoothed, range difference 
versus time for data collected during the 3-hour run of July 27. 

The data gaps are due to absence of readings during those partic- 
ular time intervals. Figure 4-20 shows the* effects of performing 
minimal smoothing. The arithmetic means of all readings observed 
in a 30-second interval, i.e., 15 to 25 data points, are given as 
a function of time. The calibration lines for these three-hour 
tests are also given in Figures 4-19 and 4-20. The calibration 
line shown in Figure 4-20 does not pass through all the smoothed 
range values? however, this^ is not unexpected considering the time 
varying nature of rubidium standard errors. 

The following implications result for the test data collected. 
Typically, 10 to 20 range differences (roughly 20 seconds) collected 
on the satellite were fit to parabolas and evaluated at the center 
point of the time span, yielding a value , say, AR^. This averaged 
range difference was reduced by L^, the calibration line evaluated 
at the same point in time. Thus AR^ - could be in error by 3 a. 

The corresponding LOP would be offset from the remote site by 3 a 
times the GDOP or typically 19 to 56 feet. 

4.7 RANGE DATA RMS NOISE 

All satellite test periods were analyzed in order to obtain estimates 
of the noise on the range measurements. The results of these analyses 
are presented in Table 4-7. The rms noise on the five drift cali- 
bration lines calculated from the pre- and post-calibration tests 
are shown, as are the rms noise for each receiver, the range 
differences between the receivers, and the rms noise on the pre- 
and post-calibration data. 

Two methods were used to determine the standard deviation of the 
ranging data; both differed from the previously described method 
of determining the a of the LOP's. In the first method, small 
continuous sections of data were fit to parabolas and the a 
computed. 

In the second method it was assumed that the ranging data was 
normally distributed. Thus, the formulation to determine a was 
that 99.2% of the time not more than 50% of the points deviated 
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TABLE 4-7. RMS NOISE 


Test 

No. 

Date 

Time 

GMT 

2 

7/12/71 

1715-1935 

3 

7/13/71 

1850-2130 

4 

7/14/71 

1700-1945 

5 

7/15/71 

1645-2030 

6 

7/16/71 

1640-1930 . 

7 

7/20/71 


8 

7/21/71 

0700-0845 

9 

7/22/71 

0700-0846 




(FEET) ON INDIVIDUAL RANGE READINGS 



3.1 12.5 

2.0 14.7 


3.0 


11.1 

11,0 

9.7 


2.1 10.0 

2.1 8..0 


0.2 7,8 

0.2 8.1 



19.2 
16 .4 
8.2 


6.7 12.2 

8.9 11.9 


Range 

•Difference 

Measurement Mode 

14 .02 

Drift Calibration Lina 

10.6 

Pre-Calibration 

7.2 

Post-Calibration 

5 . 2 

Post-Calibration 

6 . 1 

Post-Calibration 

3.9 

Post-Calibration 

13-23 

Satellite Ranging 

15.3 

Pre-Calibration 

14.9 

Post-Calibration 

14-22 

Satellite Ranging 

11.3 

Pre -Calibration 

11.3 

Pre-Calibration 

10.1 

Post-Calibration 

9-19 

Satellite Ranging 

18.43 

Drift Calibration Line 

10.3 

Pre -Cali brat ion 

8.3 

Post-Calibration 

10-20 

Satellite Ranging 

7 . 8 

Pre-Calibration 

8.2 

Post-Calibration 

8-14.8 

Satellite Ranging 


. No Data 

75.13 

Drift Calibration Line 

24 . 8 

Pre-Calibration 

17.2 

Pre-Calibration 

8.2 

Post-Calibration 

8-17 

Satellite Ranging 

13.3 

Pre-Calibration 

14.0 

Post-Calibration 

8-18 

Satellite Ranging 
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TABLE 4-7. BMS NOISE (FEET) ON INDIVIDUAL RANGE READINGS (Continued) 


Test 


Time 

Receiver 

Range 

Measurement Mode 

No. 

Date 

GMT 

A 

B 

Difference 

10 

7/23/71 

0700-0845 

2.8 

8.2 

8.3 

Pre-Calibration 


10.1 

14.0 

16.8 

Post-Calibration 




e 

- 

9-18 

Satellite Ranging 

11 

7/26/71 





No Data 

12 

7/27/71 

0100-0400 

- 

- 

10.45 

Drift Calibration Line 


4.1 

9.7 

10.3 

Pre-Calibration 




3.8 

8.1 

9.7 

Pre-Calibration 




3.7 

11.6 

12,4 

Post-Calibration 




- 

- 

B-13 

Satellite Ranging 

13 

7/28/71 

0200-0100 

- 

- 

10.25 

Drift Calibration Line 



3.7 

9.0 

10.0 

Pre-Calibration 




3.5 

8.4 

9.8 

Pre-Calibration 




3.5 

9.2 

10.0 

Post-Calibration 




- 

- 

0-14 

Satellite Ranging 

14 

7/28/71 

0330-0600 

3.7 

9.0 

9.8 

Pre-Calibration 



8.4 

8.9 

12.7 

Post-Calibration 






10-16 ' 

Satellite Ranging 




from the mean by more than + 2.66 a. . Since the number of points 
outside the mean + 2.66 a interval is assumed to be small, the 
distribution of these points is better described by a Poisson 
distribution. The mean (M) of the Poisson distribution is given 
by 


M = pN 

where 

p = 1 - q 

and N is the number of samples observed and q is the probability 
99.2%. 


The probability that k points fall outside the +2.66 a interval 
centered at the mean is given by 



-M 


Thus , 


Po 



0.5 


incorporates the 50% chance of any observation being outside the 
interval. Thus 

pN = -In 0.5 

= .69 

Since 

p = ,008 

N = 86 

The probability interval of 4.8a about the mean was chosen so that 
the number of observations would be collected in a relatively short 
interval. The difference between the smallest and largest range 
difference observed in 86 points was corrected for the relative 
drift between the rubidium standards, estimated at 600 nsec/hour. 
The observed corrected difference divided by 4.8 then yields a. 

The co-located test was considered in detail in order to determine 
the rms noise due to the receivers. The rms noise on receiver A 
was 0.2 ft during both calibration periods and 7.8 and 8.1 feet on 
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receiver B during pre- and post-calibration periods, respectively. 
The rms on the range difference was 7.8 and 8.2 feet during these 
calibration periods. Since 257 range readings on the satellite 
were available, o is calculated for three intervals of 86, 86 and 
85 data points. The average was 9.53 from the three calculated 
values of 8, 14,8 and 6.8. When was reduced by the average 
noise from the calibration periods (S^) , a remainder of 26.25 
resulted. Assuming equal noise gen erated by the receivers, the 
rms noise due to each receiver is /26 .^5^2 or 3.62 feet. This 
agrees closely with the predicted rms of 3.6 feet. 

4.8 IONOSPHERIC EFFECTS 

The collected data was analyzed in some detail in order to deter- 
mine if ionospheric effects were discernible. Since the relative 
ranging tests were not performed with this purpose in mind, 
definitive conclusions could not be reached. For example, it is 
not possible to extract information on the ionospheric-induced 
delays along the paths from ATS-5 to the receivers since such 
information is inbedded in the noise. Therefore, the only way in 
which ionospheric effects can be inferred is by considering the 
distribution of the LOP's as a function of time. 

Figures 4-11 through 4-15 show the west longitudes at one of the 
reference latitudes. The scatter on the points from test 2 and 5 
(0700-1930 GMT) is larger than for test 8 (0800-0845 GMT) which, 
in turn, is larger than for tests 12 and- 13 (0100-0400 GMT) . 

For tests 2 and 5, less than two or three successive LOP’s were to 
one side of the average LOP and, in general, any two successive 
LOP's were well separated, typically by 10 to 30 feet. This 
indicated the presence of considerable noise on the measurements. 
This noise could be due to the rubidium standards, the receivers 
or the ionosphere. During test 12, thirteen successive LOP's 
(23 minutes of observation time) were clustered to one side of 
the average LOP. These LOP's were all displaced by 30 to 50 feet 
from the mean LOP. Later, in the same test, 29 successive LOP’s 
taken during 30 minutes were all clustered to one side of the 
mean as shown in Figure 4-21. Test 13 yielded an even more striking 
result. The first 46 LOP's covering one hour and three minutes 
were to one side of their mean, while the remaining 34 LOP's 
representing one hour of data were to the opposite side. This 
effect may be due to inaccuracy in the slope of the calibration 
line (overestimating the drift between the rubidium standards in 
the first half of the test and underestimating it during the 
concluding half) or due to the motion of the ATS-5 relative to the 
predicted ephemeris which would introduce a bias-type error. 

The spread on the LOP's for tests 2 and 5 was due to considerable 
noise while most of the spread of the LOP's from tests 12 and 13 
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Fiaure 4-21. Test 12, West Longitude Crossing Points 
at Latitude 37.864253 Degrees 

ii 


resulted from a sequential motion of the LOP ' s during these tests. 
Figures 4—16 and 4—17 from tests 2 and 12 data, respectively, 
are illustrative. In these figures a single range difference was 
used to determine an LOP. Thus, no smoothing was done. Test 2 
data collected over 3 minutes 15 seconds indicates a considerable 
spread on the longitudes of approximately 200 feet. 

The rms exhibited here can be found by reducing the spread by the 
GDOP multiplier of 1.25 leaving 200/1.25 or 160 feet. Equating 
this value to approximately 5o or 66 (133 points) yields ala 
value of roughly 30. The spread for the test 12 data is 72 feet, 
for the first set and 86 feet for the second set. No dips are 
exhibited in test 12 which are exhibited by the test 2 data. 

The rms noise data presented in Table 4—7 was also considered. 
However, no statistically significant conclusions could be reached 
on the basis of that data. 
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SECTION 5 


EXPERIMENT INSTRUMENTATION 


5.1 GENERAL CONSIDERATIONS 

The test instrumentation assembled to meet the precise positioning 
requirements of the experiment included a modulator, two receivers 
and a data recording system. The design of the ALPHA II receiver 
system was based on an existing All L-band correlation receiver 
used in earlier ATS-5 ranging experiments. This experiment re- 
quired a receiver with similar precision and ranging capability. 

The ALPHA II receiver system has the following performance 
characteristics : (1) 10-MHz PRN code reception and demodulation 

capability, (2) ability to acquire and track signals from the 
spinning ATS-5 satellite, (3) computer-oriented format for re- 
cording or data-phone transmission, and (4) provision for precision 
range readout to observe small range fluctuations due to noise. 

The technical approach adopted in the ALPHA II receiver design 
was to demodulate the PRN code from the carrier without phase 
locking the carrier. The demodulating technique, explained in 
more detail later in this section, was capable of tolerating the 
small doppler anticipated due to the movement of the ATS-5 
satellite. Severe frequency uncertainties resulting from drift 
of the satellite master oscillator with temperature were accommo- 
dated by manual frequency adjustments. The digital range tracking 
loop, which is integral with the code generating and control logic 
circuitry, had natural memory characteristics that enable the 
signal to be rapidly re-acquired at the beginning of satellite 
signal burst. 

The system design characteristics which govern the ALPHA II 
receiver operation are summarized below. 

5.1.1 Modulation Design 

Coherent PSK modulation is used in the ALPHA II system. This 
implies a carrier system which transmits two signals differing 
only by a phase reversal of the carrier. The phase reversal 
represents the mark-space information. The carrier phase reversal 
produces sequential signals wich are the negative of each other. 
This means that they have a cross-correlation coefficient of -1, 
and thus are mutually uncorrelated. This modulation can be demon- 
strated to be the optimum binary system. It has been shown in 
many standard references that the coherent PSK system provides 
the lowest error probability for a given signal-to-noise ratio 
that one can expect from binary systems, including differential 
PSK, coherent frequency shift keying (CFSK) , and non-coherent FSK. 
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5.1.2 Detection Design 

The ALPHA II receiver system is a correlation design and a full 
description appears later in this section. Correlation techniques 
have long been used to extract a signal from noise or to reveal 
the presence of periodicity in an otherwise completely random 
signal. Thus, if a signal with the right properties is correlated 
with a relatively noise-free replica of itself, (i.e., autocorre- 
lation) , the random noise component of the signal will rapidly 
reduce to zero, and any periodic component will be revealed. 
Similarly, if two signals are cross-correlated, one with the 
other, an immediate result is the indication of mutual properties 
which the signals may share. 

Communication systems and radar systems apply correlation tech- 
niques to optimize performance in the presence of noise and a 
jamming environment. In the ALPHA II receiver. All uses a 
correlation technique in which the code automatically tracks by 
adjustment of the digital tracking loop. If the filter response 
of a receiver containing white noise is multiplied by a constant, 
both the signal and the noise are amplified linearly by the same 
factor, and the S/N ratio is unchanged. Therefore, to achieve 
maximum S/N ratio, the signal component must be maximized. The 
impulse response of the filter which maximizes the signal compon- 
ent and the S/N describes what is known as a matched filter. If 
the filter output is integrated over a long time period, then the 
S/N ratio can be very large. Matched filter operation can be 
obtained in this case with a correlator. The correlator multiplies 
the incoming signal plus noise by a signal which is identical to 
the anticipated signal and integrates the result. The result is 
sampled at a later time and a decision is made, concerning its 
value. The correlation principles adopted are described in the 
equipment description in this section. Appendix B describes a 
simple theoretical development of the correlation function in the 
receiver and demonstrates that the selected approach is truly an 
optimum system. 

5.1.3 Signal Design 

The code chosen for the ALPHA II receiver has a clock rate of 
10-MHz. The code generator contains a 14-stage shift register 
capable of being programmed for several maximum length codes. 

The code length necessary to provide unambiguous ranging of the 
satellite position is 16,383 bits. A shorter code of 2,047 bits 
is also available which has excellent multipath rejection and a 
relatively short search- time with a continuous signal. The 
choice of code is governed by the desired unambiguous range capa- 
bility, accuracy, the multipath performance, the acquisition 
time, system performance under doppler conditions and hardware 
complexity required for code generation. 
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It is established in Appendix C that maximum- length PN codes of 
511 bits or more are theoretically satisfactory to provide 
greater than 50 dB multipath rejection. Table C— 1 of Appendix C 
shows some of the various code lengths available from the 14-bit 
shift register code generator. This table provides an insight 
into the trade-off in search time, multipath rejection and range 
ambiguity. 

The ALPHA II system utilizes a proven correlation technique 
capable of tracking the ATS-5 L-band signal which is only present 
approximately 50 ms out of every 780 ms (one satellite 
revolution) . The utilization of digital loop to track the 
received PN code provides the memory capability between signal 
bursts. 

5 . 2 CODE CORRELATOR 

The receiver code correlator shown in Figure 5-1, compares the 
PRN code modulating the received signal with a locally generated 
duplicate of that PRN code. The modulation code is a pseudo- 
random, maximum-length code having transitions at a 10-raegabit 
per second rate . The spectrum occupancy of the code-modulated 
signal from the receiver RF front-end assembly, which feeds the 
correlator, is 20-MHz. The correlator collapses this wideband 
noise— like signal to a steady carrier when an exact phase match 
between the incoming code and the locally generated code is 
produced. 



locally Generated 
Codes to Feed the 
Decoders 


Figure 5-1. Code Correlator 
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The code correlator multiplies the received PN code modulated 
signal by a replica of the PN code and integrates the result. 

If a phase error exists between the receive code and the internal 
code, a d-c error signal is generated which causes a correction 
to be made in the phase of the internal code. When a "match" 
occurs, the correlation is complete and the phase difference can 
be measured. The code correlator operates in a sample-and-hold 
mode making phase corrections only if the decoded received signal 
exceeds a predetermined signal-to-noise ratio. The correlator 
has a predetection bandwidth of 1.6 kHz and a noise loop bandwidth 
of 50 Hz. Operation during the acquisition mode is a time search 
through all bit positions of the PN code which is terminated by 
an indication that the signal-to-noise ratio has exceeded a pre- 
determined threshold, and thus correlation has been accomplished. 
The characteristics of the code correlator are summarized in 
Table 5-1. 

The 70-MHz bi-phase PSK modulated signal is converted to 4.0-MHz 
and split into three channels. In the three correlator channels, 
the modulated waveform is multiplied with a replica of the PN 
code. The code replicas on each of two channels are displaced 
in time by plus and minus one-half bit interval from the pseudo- 
random code on the third channel. In the channel designated the 
in-phase channel, the modulation waveform is centered in time 
relative to the PN code in the other two channels. During the 


TABLE 5-1. CHARACTERISTICS OF THE CODE CORRELATOR 


Parameter 

Value 

Input 

70-MHz 

Input Bandwidth 

20 -MHz 

Input Level 

-20 dBm 

Sensitivity 

-30 dBm 

Modulation 

Bi-Phase, PSK 

Coarse Range Granularity 

100 feet 

Fine Range Quantization 

1.0 foot 

Code Length 

2047 or 16,383 

Clock Rage 

10 MB/S 

Acquisition Bandwidth 

1.6 kHz 

Control Loop Bandwidth 

50 Hz 

Multipath Rejection 

50 dB 

Duty Cycle of Tracking 

6% (ATS-5 burst) 
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acquisition mode, the time difference between the PN code replicas 
and the received modulation is not known. The function of the ^ 
correlator is to produce an adjustment of the PN code multiplying 
the in-phase channel so that it is in time coincidence with the 
received modulation waveform. This alignment is performed by 
stepping the local PN code in phase until correlation occurs at 
one of the three decoders. The output magnitude of the decoder 
increases at the instant of correlation and its presence is 
confirmed with a threshold detector. The amplitude increase can 
occur in any of the three 4.0-MHz channels depending on the time 
coincidence of the modulation waveform. 

The circuit behavior can best be explained by examining the auto- 
correlation function of each of the three decoders shown in Figure 
5-1. When the modulation waveform is matched in any of the three 
decoders the output amplitude of the decoder increases in level. 
When it exceeds a preset threshold, the correlator operating mode 
is switched from code acquisition to range tracking. 

Each decoder serves the role of comparing the incoming received 
code with a similar code locally generated in the digital sub- 
system. The 2,047 bit code correlator against itself over a 
full-code period has a single correlation peak and all other 
time side lobes are down approximately 66 dB below this peak. 
Maximum-length-code correlation characteristics are illus- 
trated in Figure 5-2. The triangle identified as A in the 



10 KCi’jbit 1 t t I 1-1^ 
Clock Pulses 


CORRtLATION PEAKS OF DECODERS A, B AND C 



Figure 5-2. Correlator Output Characteristics 
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left diagram is a correlation peak obtained by measuring the 
amplitude of the decoder output as the locally generated tracking 
code is moved through a ’’match” with the incoming code. A match 
is defined as the relative phase point at which all zeros and 
ones in the received code are aligned with corresponding zeros 
and ones in the tracking code. It is noted that the build-up of 
the correlation peak begins one bit length before the peak and 
ends one bit following the peak; thus, the correlation peak is 
two bits wide at the base. 

Decoders B and C in Figure 5-1 are identical in action to decoder 
A. Their peaks, however, are caused to occur one-half bit before 
(decoder B) and one-half bit later (decoder C) . This is accomp- 
lished by supplying the locally generated code one-half bit early 
and one-half bit late, respectively, to these decoders. 

The outputs from decoders B and C are shown in Figure 5-1 to be 
supplied to a differencing circuit. This circuit subtracts the 
decoder C output from the decoder B output. The resulting 
characteristic is shown in the right sketch of Figure 5-2. It 
is referred to as the Z curve. The mid-point of the Z curve, on 
the zero axis, is the point at which the two codes exactly match. 

It should be noted that this point coincides with the point at 
which the output from decoder A reaches its peak. 

The 4.0-MHz signals from the quadrature channel and the in-phase 
channel are filtered with crystal filters having a bandwidth of 
approximately 1.5 kHz. The in-phase and quadrature channels at 
the crystal filter outputs contain only the carrier; the modula- 
tion sidebands are suppressed by the filter selectivity. 

The voltage amplitude of the carrier at the output of the crystal 
filters contains the information concerning the time coincidence 
or correlation of the received modulation and the locally generated 
PN code. When the two range code signals are in time coincidence, 
the amplitude of the 4.0-MHz carrier in the in-phase channel will 
be maximized and the 4.0-MHz carrier in the quadrature channel 
will be minimized. 

The processing gain of the system is inversely proportional to 

filter bandwidth, and is computed by the ratxo of the incoming 
signal spectrum to the filter bandwidth. Thus, the 20-MHz spread- 
spectrum signal and the 1.6-kHz filter produces a processing gain 
of 41 dB. 

Other factors influencing the choice of the filter bandwidth are 
minimum search time and tolerance to frequency uncertainty. The 
search time is directly affected by the time required to build up 
(j 0 tectable output from the decoders. As the locally generated 
tracking code is stepped one bit at a time in its search mode, 
it must remain in each phase position for a time sufficient to 
integrate the signal in the filter. 
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The correlator configuration shown in Figure 5-1 perforins its 
decoding independent of the doppler components on the code- 
modulated carrier. Tolerance to frequency offset extends to a 
maximum value of +0.8 kHz, The correlator design enables the 
decoding action to take place despite the presence of small 
frequency changes on the carrier. The frequency effects are 
identical in both channels of the correlator and if the carrier 
falls within the bandpass filters, the decoding is independent of 
the carrier frequency. Code tracking was demonstrated by this 
receiver in the operation with the ATS-5 satellite for doppler 
produced by the satellite motion (approximately 35 feet per 
second). Severe frequency drifts due to the ATS-5 oscillator 
(up to approximately 200 Hz/min) were corrected for by tuning 
of the 2nd L.O. injection signal. 

At this point the code phase error information is in the form of 
amplitude. The ALPHA II correlation system converts this ampli- 
tude function into phase information, so that variations in the 
received signal amplitude will not influence the range measure- 
ment accuracy . The quadrature channel is amplified in a limiting 
amplifier and the in-phase .channel signal is amplified in a similar 
limiting amplifier. The two limiting amplifiers are phase-matched 
to witin 2® so that their phase response characteristics do not 
produce range errors. The limiting amplifier outputs are then 
multiplied in a phase detector which is capable of generating an 
error voltage proportional . to the phase error between the two 
inputs. The error voltage is integrated in a low-pass filter and 
amplified. This DC signal contains the sense and magnitude of 
the alignment error between the local PN code generator and the 
received range code y/aveform. 

The characteristic which describes the output of the balanced 
detector as the received code and locally generated code near the 
optimum tracking position is the Z curve shown in Figure 5-2. 

The 0 point midway on the Z curve is the balance position where a 
phase match exists between the two codes. As the phase match 
changes, a DC output is generated by the Z characteristic. The 
error signal from the balanced detector is band limited by a low- 
pass filter and sampled by an analog-to-digital converter. This 
circuit converts the DC error voltage into digital logic pulses 
to control the phase positioning of the tracking codes. 

5.3 DIGITAL TRACKING SUBSYSTEM 

5.3.1 General 

A correlation system requires that the receiver have exact knowledge 
of the code transmitted by the sender. This requirement is met by 
a reference code generator in the receiver identical to one in the 
modulator. Under normal circumstances where the receiver is oper- 
ating independent of the transmitter in a remote location, the 
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reference code generator in the receiver is driven by its own 
master clock. The block diagram of the ALPHA II code correlator 
digital subsystem is shown in Figure 5-3. The system has three 
basic functions to perform: (1) acquire and lock the internally 

generated PRN code to the received modulation, (2) measure the 
phase difference between the received code and the internal 
reference, and (3) display this measurement to the operator and 
make it available for external recording. 

5.3.2 Operation 

Two PRN code generators are used in the ALPHA II correlator. 

The first is a reference generator, while the second, the tracking 
code generator, is locked to the incoming signal. Each code 
generator is a 14-stage shift register with feedback taps con- 
figured to generate either 16,383- or 2,047-bit maximum length 
PN codes. The state (10000... 00) is decoded as a synchronization 
point, and a pulse is produced whenever this state is reached. 
Since the input clock rate is 10-MHz, the cycle time of the code 
generators is 204.7 or 1638.3 microseconds, depending upon the 
selected code length. The tracking code generator develops the 
three waveforms provided to the decoders. These three versions 
of the PRN code are offset from each other by one half of a code 
bit (50 nsec) . 

The phase selector provides a convenient method of shifting the 
tracking code generator to provide the equivalent of 40-MHz 
resolution. This technique, illustrated in Figure 5-4, produces 
four 10-MHz square waves displaced from each other by 25 nano- 
seconds (90®) . By selecting the waveform nearest the phase of 
the received signal as the clock for the tracking code generator, 
the correlation will be within +1/8 bit (+12.5 ns) and the error 
voltage remains in the linear region of the Z curve. 

The phase of the internal code can be easily shifted by changing 
from one clock stream to another, as the last two waveforms of 
Figure 5-4 indicate. If the negative edge of the waveform shifts 
the code generator and phase 2 is the active clock phase, there 
is a shift at the 0 and 100 nsec points. If the clock switches 
from phase 2 to phase 3 at 25 nsec, the "retard" waveform results, 
with shifts at the 0 and 125 nsec points. The subsequent shifts 
will occur at 100-nsec intervals provided no further clock changes 
are made. Thus, the effective phase of the code generator has 
been delayed by 25 nsec or 1/4 of a bit. Similarly, switching 
from phase 2 to phase 1 produces the "advance" waveform with 
shifts at the 0 and 75 nsec points and advances the code phase 
25 nsec. By successively stepping through the clocks (2, 3, 0,1, 2 
or 2, 1,0, 3, 2), the phase of the code can be changed to any 
desired position. 
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Figure 5-3* Simplified Block Diagram of Digital Subsystem 
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Figure 5-4. Clock Phase Timing Diagram 

If the internal code and the received code are within one bit 
(+100 nsec) of phase lock, an indication is given whenever the 
received signal strength is above a preset threshold. This 
receiver threshold should occur each time the satellite antenna 
illuminates the receiving sites, or approximately 50 milliseconds 
out of every 783 milliseconds. If the threshold signal is not 
received for two seconds , the control logic activates the search 
generator, changing the tracking code generator phase until a 
threshold (i.e., correlation) is achieved. 

When the receiver is in the acquisition or search mode, the code 
tracking generator is sequenced through the entire 2,047 phase 
positions relative to the reference generator. When the correla- 
j^j^j- 0 shold detectors provide a positive indication, the 
acquisition mode is stopped and correlation is within +1 bit. 

The range computing circuitry reverts to the tracking mode after 
the initial acquisition of the received code. In the tracking 
mode the 10-MHz clock is phase shifted (as previously discussed) 
in 90® increments to achieve the effective granularity of a 40- 
MHz clock system. The clock phase shifting is initiated when the 
correlator phase detector output exceeds the preset threshold 
criteria. 
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5.3.3 Search Routine 


The rotating ATS- 5 satellite imposes a unique search-sequence 
requirement on the range acquisition and tracking circuitry. 

The usable satellite illumination time at the Wallops Station 
is slightly more than 50 milliseconds of every 783 milliseconds 
and the correlation system integration time is approximately 1 
millisecond. These constraints established the search sequence 
shown in Figure 5-5, which achieves system range code phase lock- 
on with the ATS-5 satellite within 135 seconds (worst case) for a 
2,047-bit code length. If the signal is lost for more than 2 
seconds, as measured by the correlator threshold detector circuit, 
a search mode is initiated. During normal tracking mode, the 
inherent digital memory of the range tracking code generator will 
enable re-acquisition of phase lock within the first satellite 
illumination burst. 

The search pattern times in the receiver are shown in Figure 5-5. 
Beginning at an arbitrary zero reference point for the tracking 
code generator phase, 12.5-bit positions (fifty 25-nsec steps) 
are searched during a 51,2-msec period. These same positions 
are searched repeatedly for 821 msec, which is sufficient time 
(in excess of 780 msec) to assure that at least one signal burst 
has been received from the satellite. The search pattern then 
moves to the next 12.5-bit position for the same length of time, 



Figure 5-5. Code Search Pattern 
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progressing through all the positions in approximately 135 seconds. 
When correlation is achieved, a receiver threshold signal switches 
the control logic from the search mode to a tracking mode. The 
worst-case search time is determined by calculating the number of 
12.5-bit searches, that is 2,047/12.5 - 163.7, and multiplying by 
0.821 second. The maximum search time, then, is 135 seconds. 

The presence of a continuous signal (i.e., a non-spinning satellite) 
would reduce the acquisition time by a factor of approximately 16:1 
by allowing a simplified non-repeating search pattern to be 
implemented . 

5.3.4 Tracking Mode Routine 


In the tracking mode, the system operates along the Z curve, 
remaining within +l/8-bit of the center point (see Figure 5-6) . 
The fine range circuits accomplish this by performing analog-to- 
digital (A/D) conversions of the correlator DC error voltage and 
comparing these to preset bounds. If the positive bound (A/D = 
25) is exceeded, a signal is sent to the phase selectors to 
retard the code phase one quarter bit (25 nsec). Similarly, 
exceeding the negative threshold (0) initiates a 25 nsec phase 
advance in the code streams. 

A/D values measured between these two bounds are valid range 
readings between 0 and 25 nsec. The clock phase used to drive 
the tracking code generator provides an indication of the 25-nsec 
increments. Addition of 0, 25, 50 or 75 nanoseconds to the A/D 
value is accomplished in the fine range logic and produces the 
two least significant digits of the range reading. 



Figure 5-6 . " Z " Curve 
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The range counter determines the remaining range digits by 
measuring the phase difference between the reference code generator 
and the tracking code generator. The counter starts incrementing 
one count every 100 nsec when it receives the sync pulse from the 
reference generator and stops upon receipt of the sync from the 
tracking generator. 

5.3.5 Data Output and Disp lay 

As shown in Figure 5-11 the ALPHA II code correlator displays the 
range reading to the operator via a 7“segment LED numeric display. 
Storage elements are provided so that the most recent valid reading 
is displayed. The data output also is available as a serial bit 
stream containing the range reading and selected status indicators. 
The data format consists of four USASCII words using odd parity 
which allows the information to be transmitted via commonly 
available data sets and received by the recording device. The 
signals comply to EIA Standard RS232C levels for ease of interface. 
The four data words are transmitted when the receiver threshold 
signal drops if a valid reading had been made while the threshold 
was present. 

5.4 RECEIVER FRONT END 

A functional diagram of the receiver front end and IF section is 
presented in Figure 5~7. The RF preamplifier is a broadband 
transistorized unit with a measured noise figure of 4.4 dB and a 
bandwidth greater than 50-MKz. The gain at 1550 MKz is specified 
as 25 dB minimum. The down-converter is a * Schottky-barrier 
junction - diode mixer with a conversion loss of 8 dB and an 
optimum noise figure of 8 dB when driven with a local oscillator 
injection power level of 1 to 10 milliwatts. The total RF-to-IF 
conversion gain is approximately 110 dB. 

The injection chain which supplies the 10-milliwatt, 1480-MHz 
signal to the down-converter contains several transistor multi- 
plier stages, a step-recovery diode multiplier and a double 
balanced mixer. The injection chains require 5-MHz and 1-MHz 
input signals from the master frequency source. The block 
diagram also shows generation of the two injection signals 
required by the code correlator and the ranging computer 
subsystems . 

The 70-MHz IF input impedance is 50 ohms with a nominal sensitivity 
of approximately -84 dBm. The physical implementation of the 70- 
MHz IF output is such that other equipments may be interfaced. 

The IF section was designed with a controlled bandwidth of 25-MHz 
and a total gain of approximately 90 dB. A broadband AGC detector 
and amplifier were incorporated to provide gain control of the IF 
based upon a total detected noise level at the IF output. The 
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IF TEST POINT 



Figure 5-7. Functional Block Diagram of RF Subsystem 
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wide-bandwidth and stable-phase characteristics of the RF and IF 
sections are necessary in order to minimize the influence of 
environmental variations upon the signal delays through these 
stages. The design techniques associated with the RF and IF 
stages and injection chains are conventional but use the latest 
state-of-the-art devices. The multiplier stages are well filtered 
to remove all spurious products naturally generated in such a 
configuration . 

5 . 5 MODULATOR 

The ALPHA II modulator installed at the NASA Rosman STADAN Station 
is an all-solid-state signal generator used to phase-modulate the 
70-MHz IF carrier furnished to the NSAS Rosman transmitter system. 
The 70-MHz IF carrier is bi-phase modulated in a PSK manner by the 
PN range code, and is compatible with the STADAN interface stand- 
ards. The modulator is self-contained and requires only a 5-MHz 
master clock signal for operation. The stability of the clock 
source governs the stability of the 10-megabit pseudo-random code. 

A built-in frequency multiplier chain converts the 5-MHz signal to 
a stable 70-MHz carrier by means of a X 7 and a X 2 multiplication. 
The resultant carrier is then multiplied by the 10 -megabit PRN code 
in a phase modulator. The phase modulator shifts the carrier wave 
by 180 degrees at each transition of the PN code, thus producing 
the bi-phase PSK modulation of the carrier. The modulator is 
capable of generating either the 2,047 or 16,383 bit code. The 
normal experimental configuration of the modulator at the NASA 
Rosman Station consisted of feeding the output of a frequency 
synthesizer into the phase modulator at an IF signal frequency 
of 70.1-MHz. This was necessitated by the biased bandpass char- 
acteristics of the ATS-5 C-band to L-band cross-strapped trans- 
ponder. A similar modulator was also employed at the Wallops 
Station test configuration for pre- and post-calibration test 
purposes , 

5.6 RUBIDIUM FREQUENCY STANDARDS 

In general, four types of frequency standards are commonly used; 

(1) Atomic hydrogen maser 

(2) Cesium atomic-beam-controlled oscillator 

(3) Rubidium gas-cell-controlled oscillator 

(4) Quartz crystal oscillator 

The first two are referred to as primary frequency standards and 
the last two as secondary -frequency standards. The advantages and 
limitations of each standard type are summarized in Table 5-2. 

The requirements of the positioning experiment were basically 
portability and good short-term stability. Therefore, the rubidium 
gas-cell-controlled oscillator was recommended for the test. 
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TABLE 5~2. CHARACTERISTICS OF FREQUENCY STANDARDS 


Standard 

Construction 

Advantages 

Limitation 

Atomic 

Active Maser 

Primary 

Size and 

Hydrogen 

Maser 


Standard 

Long and Short 
Term Stability 

weight 

Cesium Atomic 

Atomic beam 

Primary 

Short-term 

Beam Resonator 

Controlled 

Oscillator 

interaction 
with fields 

Standard 

Long-term 

Stability 

stability 

Rubidium Gas 

Gas buffered 

Compact and 

Secondary 

Cell Resonator 

resonance 

lightweight 

Standard 

Controlled 

cell with 

Good short- 
term 

stability 

Requires 

calibration 

Oscillator 

optically 
pumped state 
selection 

Quartz Crystal 

Quartz 

Very compact, 

Long-term 

Oscillator 

crystal 

light, rugged 
and inexpensive 

stability 


Two identical Varian Associates Model R-20 Rubidium Frequency 
Standards were employed during the experiment as the frequency 
sources. The frequency standard consists primarily of a servo 
system in which a crystal oscillator is locked to an electron 
hyperfine transition frequency found in the rubidium atom. This 
provides the long-term stability reference, and the crystal 
oscillator provides the short-term stability and spectral purity. 

The Model R-20 Rubidium Frequency Standard has the published 
specifications summarized in Table 5-3, The instrument is 
designed for highly portable applications, employing a 24 to 30 
volt d-c battery-pack input for standby or field-use. During 
the experiment, the standards were operated from 24-volt, 200- 
AH batteries which were recharged following each test period. 

The battery for the remote system also supplied power for the 
fans and code tracking logic in the code correlator via a solid- 
state DC/AC inverter. Continuous power was required to preserve 
the code reference after calibration. 

5,7 DATA RECORDING SUBSYSTEM 

/ 

The serial data stream generated by each code correlator contains 
the range reading which was permanently recorded on magnetic tape 
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TABLE 5“3. SPECIFICATIONS FOR VARIAN MODEL R-2 
RUBIDIUM FREQUENCY STANDARD 


Parameter 

Value 

Output Frequencies 

5 MHz, 1 MHz and 100 kHz, 
simultaneously at front and rear panel 

Output Level 

1 V rms, into 50 ohms 

Frequency Stability 
(Standard Deviation) 

1 X 10"ii for a one-second averaging 
time 

—11 

5 X 10 over any one-year period 

Environmental 

Stability 

Above long- and short-term stability 
specifications will be maintained 
over the following conditions : 


• Temperature: 0° to 50 ®C 

o Humidity: 0 to 95% 

9 Input Voltage: +15% of nominal 
0 Load: Changed from open to short- 

circuit at any other terminal 
0 Magnetic Field: Less than 5 x 10 

frequency change for any 
orientation in the Earth’s field 

Radio Frequency 
Interference 

RFI shielded 

Warm-Up Time 

Instrument accuracy is +1 x 10"’^^ 
after less than a one-hour period 

Alarm Indicator 

Front panel continuity alarm light 
visually indicates that the output 
frequency is locked 

Input Requirements 

24 - 30 V dc nominal; less than 1 amp 
operating; 2,5 amps warm-up 

Packaging 

Dimensions are 4-7/8 inches by 7-5/8 
inches by 19-9/16 inches. Weight is 
20 pounds. 
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by the data recording subsystem. The tape recording subsystem 
consisted of a Kennedy Model 1510 incremental tape recorder 
driven by an All tape control unit, Model TC-701. The All tape 
control unit (TCU) , shown in the functional block diagram of 
Figure 5~8, is specifically designed to record the serial data 
from each code correlator. It accepts four consecutive ASCII- 
odd bytes (RS 232C or TTL compatible) at 110 baud, on each of 
two channels, cheeks the parity, formats the 28 bits of data, 
and writes the data on magnetic tape via the 7-track (6 data, 1 
parity) Kennedy incremental tape recorder. Successive blocks of 
four bytes on each TCU channel must be separated by a minimum of 
50 milliseconds. The delay between the arrival of data blocks 
on the two channels must be less than 50 milliseconds. The 
parity bits must arrive correctly on each channel for any infor- 
mation to be recorded. Inter-record gaps are inserted between 
information sets. 

As shown in Figure 5-8, the TCU receives four ASCII-odd bytes 
from the code correlator, strips the four control bits from each 
11-bit byte, stores the stripped data in a shift register, checks 
parity, verifies that all four bytes have arrived, and (for two- 
channel operation) verifies that the block of data on each channel 
has arrived within 50 milliseconds of the other channel. If these 
criteria are satisfied, the multiplexer formats the stored data 
and provides the proper write-clock and gap pulses for the 
magnetic tape transport. The data format generated by each code 
correlator is shown in Figure 5-9. 

The data format appearing on the magnetic tape can include a 
"time" representation. When the mode switch is placed in the 
RECORD position, an elapsed time counter is initiated. The 
maximum time indication before repeating is 59 minutes 59.9 
seconds and the resolution is 100 milliseconds. The time indica- 
tion is automatically recorded with each data file as three bytes. 
The data file written on 7-track magnetic tape is shown in Table 
5-4. The "C" or parity bit is automatically generated by the 
tape recorder and is therefore not shown in the table. 

During the experiment, the data link to the TCU from the fixed 
site receiver was simply a direct cable since the code correlator 
output and the TCU input are RS232C compatible. The data link 
for data transmission from the remote receiver site was via a 
pair of RS232C compatible acoustic couplers and a commercial 
telephone line. This is shown in simplified form in Figure 5-10. 

5.8 EQUIPMENT PHYSICAL DESCRIPTION 

The three main assemblies which comprise each receiver system 
are shown in the photographs of Figure 5—11. Each unit is con- 
tained in a standard 19-inch wide rack-mounted enclosure. 
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Figure 5-8. Functional Block Diagram of Tape Control Unit 
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Figure 5-9. Correlator Output Data Format 


TABLE 5-4. MAGNETIC TAPE DATA FORMAT 
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Figure 5-10. Experiment Data Transmission Link 

The RP front end is a 5-1/4 inch high drav/er having a depth of 
17 inches. The packaging consists principally of shielded RF 
subassemblies mounted internally on an aluminum shelf. RF signal 
input and output connectors are located in the rear, together 
with A-C power connection. This unit contains its own internal 
power supply. The receiver down-converter weighs about 25 pounds 
and occupies 1.3 cubic feet. 

The code correlator unit is 7 inches high by 21 inches deep. 

All input and output connnections are on the rear panel. The 
front panel includes test points, control switches and an LED 
readout display of the range reading in nanoseconds. The readout 
display was included to provide a "real-time" check of operation 
and reasonableness of data during testing. A "track" light 
on the panel is illuminated on each pass of the satellite beam 
to indicate signal presence in the receiver channel. Figure 5-12 
shows several of the RF subassemblies alongside the printed 
circuit cards which comprise the digital ranging circuits. 

The power supply furnished the various voltages required by the 
code correlator. Meters are included to monitor current and 
voltage. This unit is 7 inches high and 15 inches deep, occupying 
a total of 1.2 cubic feet. The weight is 35 pounds. 

Figure 5-13 shows the All modulator which was installed at the 
NASA Rosman Station, and is also a 19 -inch-wide , rack-mounted unit. 
The 3-1/2 inch by 19 inch Tape Control Unit is also shown in 
Figure 5-13. 
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APPENDIX A 


DATA SUMMARY 


In order to place the data reduction and summary method in proper 
perspective, the flow of information is traced from the time the 
two correlators make their measurements until the lines of 
position (lop’s) are determined u 

Assuming that at a given time the code • correlators make valid 
measurements on the respective ranges to the satellite, the 
range values are sent to the tape recorder subsystem and recorded. 
The decimal format of these binary numbers is shown in the last 
24 columns of Table A~1 (a segment of data from Test 12, July 27), 
divided in two groups of 12 columns each. This is the rawest form 
in which data is available. The computer, under control of the 
Raw Tape Decoder Program, decodes the raw data bits into time, 
status indicators and coarse and fine range values for both 
receivers. At this time record numbers are assigned, range values 
are converted to feet, and the range difference is determined. 

A magnetic tape of the converted data is then generated and 
submitted to the Lines of Position Computer Program. The first 
operation of this program is to align the recorded time with 
Greenwich Mean Time. Since the time counter of the incremental _ 
tape recorder overflows at 59 minutes 59.9 seconds, all succeeding 
recorded times must be increased by one hour until the next counter 
overflow occurs, after which two hours are added, etc. In addition, 
this program aligns the zero time, corresponding to tape recorder 
turn-on, with the actual time of the start of test. Thus, the 
true Greenwich Mean Times for the data presented in Table A-1 are 
achieved by adding a correction factor of one hour, 20 minutes to 
the record times indicated. 

The next step is the determination of the parabolic fits to the 
sets of satellite latitude, longitude and geocentric range. 

These parabolas for test 12 are: 

9 Latitude = —,723 + ,109t + ,019t^ deg. 

• Longitude = 105.216 - .048t + .OOlt^ deg. 

• Range = 22717.705 - 2.600t + 1.220t2 rani. 

The calibration data is fit to a straight line: 

L., = -79619.478 - 605.738t 

JL ^ 

The time t is in hours and decimal hours since 0000 GMT. 
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TABLE A~l. SAMPLE CODED AND DECODED 


time 

HH HH sS.Sj 

PECOPD 

N(jm8Fb 

sta 
1 2 

CAI 
1 2 

S/T 
t ? 

Th9 
1 2 

CHM 
1 2 

COA9SfcL rlNfl 
<NsEC) 

C0A9SE2 rilvf2 
t^sEC) 

0-58-49.40 

2C53 

1 

0 

0 

n 

1 

1 

1 

1 

1 

1 

561 

22 

1172 

40 

n-5B-5u . In 

2054 

1' 

u 

0 

n 

1 

1 

1 

1 

1 

1 

561 

5 

U72 

26 

0-58-50.90 

2C55 

1 


0 

0 

1 

1 

1 

1 

1 

1 

560 

91 

1172 

24 

0-58-51. Tf) 

2C56 

1 

•} 

0 

0 

1 

1 

1 

1 

1 

1 

56o 

70 

1172 

9 

n-56-52 . 5 n 

2C57 

1 

\) 

n 

h 

1 

1 

1 

1 

1 

1 

56q 

56 

ll7l 

9l 

n-50-54 .In 

2C58 

1 

It 

n 

n 

t 

1 

1 

1 

1 

1 

56CI 

21 

ll7l 

62 

0-50^54. 8n 

2C59 

1 

j 

n 

n 

1 

1 

1 

1 

1 

1 

56c 

15 

1171 

40 

0-50-56 . 4fi 


1 


0 

Fl 

1 

1 

1 

1 

1 

1 

559 

76 

ll7l 

4 

o-50i59. on 

?U51 

1 

1, 

0 

n 

1 

1 

1 

1 

1 

1 

559 

3^ 

1170 

77 

0-58-50. 7n 

?062 

1 


0 

n 

1 

1 

1 

1 

1 

t 

559 

22 

ll7o 

63 

0-50-59. 5n 

?C63 

1 

!l 

D 

0 

1 

1 

i 

1 

1 

1 

559 



57 

0-59- C'.3n 

?C64 

1 

.1 

n 

0 

1 

1 

1 

1 

1 

1 

558 

92 

1170 

41 

0-59- 1 . 9fi 

?C65 

1 


D 

n 

1 

1 

1 

1 

1 

1 

558 

66 

116 9 

02 

0-5?: ?,6o 

?C66 

1 

1. 

0 

0 

1 

1 

1 

1 

1 

1 

558 

46 

1169 

77 

n-59- 4 .?o 

2067 

1 

>J 

n 

n 

1 

1 

1 

1 

1 

1 

558 


1169 

38 

n-59- 5.00 

2C68 

1 


n 

n 

1 

1 

1 

1 

i 

1 

550 

1 ’ 

1169 

29 

n-59i 5.8n 

2069 

1 

.1 

0 

n 

1 

1 

1 

1 

1 

1 

557 

64 

1169 

10 

0-59- ft.5o 

2C70 

1 

fl 

n 

n 

1 

1 

1 

1 

1 

1 

557 

69 

1169 

18 

n- 59 ; 7.3n 

2071 

1 

l1 

0 

0 

1 

1 

1 

1 

1 

1 

557 

43 

1160 

9l 

0 - 59 ; 8.9o 

2D72 

1 

ij 

n 

n 

1 

1 

1 

1 

1 

1 

557 

16 

1168 

54 

n-59;n'.4o 

2073 

1 


n 

n 

1 

1 

1 

1 

1 

1 

556 

04 

1160 

2 

0-59-1 1 . 2n 

2074 

1 

u 

0 

n 

1 

1 

1 

1 

1 

1 

556 

60 

U66 

13 

n-59;i 2 .80 

2075 

1 

'J 

0 

R 

1 

t 

1 

1 

1 


556 

49 

1167 

05 

0-59;i 3 , 6n 

2076 

1 

:• 

n 

n 

1 

1 

1 

1 

1 

1 

556 

l7 • 

1167 

66 

n-59;i 4 . 3n 

2077 

1 

1 

n 

n 

1 

1 

1 

1 

1 

1 

556 

6 

1167 

28 

0-59;i5.9o 

2078 

1 

[• 

0 

n 

1 

1 

1 

1 

1 

1 

555 

89 

1167 

16 

n-59;i6. 7n 

2079 

1 

■J 

0 

n 

1 

1 

1 

1 

1 

1 

558 

53 

1167 

2 

n-59;i7 , 5o 

2000 

1 

0 

n 

0 

1 

1 

1 

1 

1 

1 

555 

4? 

1166 

94 

o-59;20 . isn 

2CP1 

1 

'1 

a 

n 

1 

1 

1 

1 

1 

1 

554 

77 

1166 

l9 

0 - 59 ;?! .40 

2082 


J, 

0 

n 

1 

1 

1 

1 

1 

1 

554 

67 

1166 

1 6 


2083 

1 

1 ; 

0 

ft 

1 

1 

1 

1 

1 

1 

554 

49 

1165 

9n 

n-59;?2.9n 

2084 

1 

:t 

0 

0 

1 

1 

1 

1 

1 

\ 

554 

40 

1165 

76 

n-59-24 . 5n 

20 85 

1 

• 1 

n 

fl 

t 

1 

1 

1 

1 

1 

554 

l4 

1165 

55 

0 -59-25 . 3n 

2086 

1 

1' 

0 

fl 

1 

1 

1 


1 

1 

553 

09 

1165 

37 

0-59-2? . 6o 

20 8 7 

1 

»i 

n 

n 

1 

1 

1 

1 

1 

1 

553 

36 

1164 

79 

f[-59;2fl . 4o 

2080 

1 


0 

n 

1 

1 

1 

1 

I 

1 

553 

1? 

1164 

59 

(1-59;29. 2o 

2u89 

1 

r 

0 

n 

1 

1 

1 

1 

1 

1 

553 

11 

1164 

59 

o-59;2V . 9n 

2U«0 

1 

_ 

(} 

0 

n 

1 

1 

1 

1 

1 

1 

552 

91 

1164 

44 


FROM TEST 12 (JULY 27, 1971) 




















It should be noted here that the range difference (AR) in Table 
A"1 equals M 2 - , i.e., the measurement by the B or second 

receiver, reduced by the measurement of the fixed, or second 
receiver. However, the Lines of Position Program calculated the 
lop's using - M 2 thus causing AR to change sign. 

The ephemeris and calibration line are evaluated at 2 hours, 18 
minutes, 49.4 seconds, corresponding to the first entry in Table 
A-1.' The relative position calculations (see Appendix B) are 
then made and the longitudes corresponding to the reference 
latitudes of 37.864253 and 37.864300 degrees are recorded. These 
longitudes are 75.5075923 and 75.5076569, as given in Table A-7 
for Test 12, (page 2 of 3). The other individual records of 
Table A“1 are also all used to determine the LOP ' s of Table A-7; 
but, no smoothing on the observed range differences was attempted 
in this case. The third longitude shown in the table is the 
running average of the longitudes corresponding to the first 
reference latitude. A similar segment of unsmoothed range data 
is given for Test 2 in Table A-3. 

In general, twenty consecutive range differences were fitted to 
parabolas. These parabolas were then evaluated at the center of 
the timespan. The LOP's were then determined as described above. 
The resulting LOP's were then plotted in Figures 4-6 through 
4-10 for data collected on July 12, 15, 21, 27 and 28, 
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TABLE A-2. EPHEMERIS AND LOP SUMMARY FOR TEST 2 
(JULY 12, 1971) (PAGE 1 OF 2) 


ATS-5 

EPHEMERIS 





TIME 

SAT RANGE 

LATITUDE 

LONGITUDE 



GMT 

NMI 

DEG 

DEG 



lG3000.f1 

22781.02 

0.150 

105.251 



1 70000.0 

22778.04 

0.038 

105.250 



173000.0 

22774.59 

-0.075 

105.268 



1P0000.0 

22770.74 

-0.187 

105.282 



IR3000.0 

22766.55 

-0.295 

105.294 



1 5)0000.0 

22762. 1 1 

-0.399 

105.303 



f>ARAPOUC 

FIT TO EPHEMERIS-LI NEAR TO CAL DATA 



CONSTANT 

LINEAR - QUADRATIC 

RMS 

LATITUDE 

5.1302525 -0.5724861 0 

.0042857 0. 

00105 

LOaGITUDE 

105. 0968293 0.2165545 -0 

.0052851 0. 

00028 

GECC RANGE 22676 .525 7P55 IB. -4257095 -0 

.7330221 0. 

03894 

CALIB LINE B5660. 2156161 -547.5374731 

14. 

02497 

TIME 

range DIF. 

CAL. VAL. 

LONGITUDES (DEG) 

GMT 

(FT) 

(FT) 

POINTl 

P0INT2 

CUMUL. AVQ. 

1R5645.8 

100259.69 

79286,54 

75.5075091 

75.5075732 

75.5075091 

IB5724.0 

;00262.80 

79280.75 

75.5075602 

75.5076244 

75.5075347 

1B5755.6 

100262.91 

79276.23 

75.5075854 

75.5076495 

75.5075516 

185855.4 

1 00199.66 

79270.18 

75,5072198 

75.5072840 

75.50746R6 

185904.5 

100244,25 

79265.45 

75.5075265 

75,5075904 

75.5074802 

185938.8 

100239.4.9 

79260.23 

75.5075247 

75.5075889 

75.5074876 

190021.7 

100224.50 

79253.71 

75.5074645 

75.5075287 

75.5074845 

190056,0 

1 00209.99 

79248.49 

75.5074027 

75.5074669 

75.5074741 

190130,3 

108250.67 

79243.27 

75.5075613 

75.5076255 

75.5074P3P 

190204.6 

100219.54 

79258.06 

75,5075195 

75.5075837 

75.5074874 

190258.8 

100164.17 

79232.86 

75,5071990 

75.5072632 

75.507461 1 

J90516.6 

100166.46 

79227.1 1 

75.5072447 

75.5073088 

75.5074431 

190404.6 

100186.19 

79219.81 

75,5074086 

75.5074728 

75.5074405 

I9l,'501 .ff 

1 00214.99 

79211.1 1 

75.5076375 

75.5077015 

75.5074545 

190^40.4 

100204 .87 

79205.23 

75,5076055 

75.5076697 

75.5074646 

1906;»8.5 

100191 .57 

79200.96 

75.5075449 

75,507609 1 

75.5074696 

190655.9 

1S0I59.1 1 

79196.34 

75.5075656 

75.5074298 

75.5074635 

190710.8 

100167.57 

79191.49 

75.5074440 

75.5075082 

75.5074624 

190754. & 

100176.09 

79187.93 

75.5075183 

75.5075825 

75.5074654 

1908 I0<.0 

100168.41 

79182.48 

75.5074995 

75.5075637 

75,5074671 

190!?»47.8 

100143.54 

79176.73 

75.5073741 

75.5074383 

75.5074626 

190917.8 

100155.57 

79172.17 

75.5074621 

75.5075263 

75.5074626 

190959.7 

1001 16.12 

79168.84 

75.5072443 

75.5073085 

75.5074531 

191021 .5 

100110.98 

79162.51 

75,5072463 

75.5075105 

75.5074445 

191101.5 

100124.59 

79156.40 

75.5073641 

75.5074283 

75.5074413 

191 149.8 

100117.59 

79149.05 

75.5073612 

75.5074254 

75.5074382 

191224.9 

100131.60 

79143.71 

75.5074785 

75.5075427 

75.5074397 

191300.7 

100135.06 

79138.27 

75.5075299 

75.5075941 

75.5074429 

191535.4 

100125.24 

79133.29 

75,5074951 

75.5073595 

75.5074447 

191415.5 

100081.40 

79126.89 

75.5072538 

75.5073180 

75.5074383 




TABLE A-2. EPHEMERIS AND LOP SUMMARY FOR TEST 2 
(JULY 12, 1971) (PAGE 2 OF 2) 



TIME 

RANGE PIE, 

CAL. VAL. 

GMT 

(FT ) 

(FT ) 

101443.2 

1001 14,60 

701P2.68 

101522.5 

1 00f’c3 .6 1 

701 16.70 

19161 1 .6 

100106.25 

70 1 0o . 23 

191643.6 

100106.73 

79104.37 

191721 .0 

100101.10 

70008 . 68 

101750.6 

100100.47 

700^4 , IP 

101821.0 

1000^3 .32 

7o0fq,55 

191002.3 

100082. 12 

70083.27 

191029.6 

100070 

79079. 12 

192006.2 

100081 . 19 

79073.55 

1 02044.4 

100095.34 

70067.74 

102118.7 

1 0008 1 .70 

70062,53 

102147.3 

100067.57 

70058. 15 

102217.1 

100074. 13 

70035.64 

192247.3 

100058.42 

70040.02 

102322.6 

1 00050.50 

70043 ,68 

102540, 1 

10005 1 .60 

70050.65 

1 02422.6 

1 0003 1 .07 

70054.36 

102431 .4 

100067.21 

70050. 18 

102321 .8 

100061 .02 

79023.33 

192530.2 

100052.08 

70021.23 

192612.5 

1 00045.68 

700 17.84 

102635.1 

100040.26 

79014,40 

102700.8 

1 00037.00 

790 10.49 

102728.8 

I 00014, 74 

70006.24 

102803.4 

1 00026.00 

700.00,67 

192842. 1 

1 00027.26 

78095.00 

102906.2 

100022. 10 

78991 .42 

192931.2 

100018.03 

78987.62 

102058 .4 

1 00016.26 

78983.48 

103021 . 1 

100008 , 1 0 

78980.03 


LONGI TUDES (DEG) 

POIMT! P0INT2 CUMUL. AVG. 

30755'^? 75.^P1739')<5 73.5074366 

307^06? 73.3073705 75.50743P7 

3073360 73.3076002 75.5074416 

30753:*^ 75.5075Q61 73.3074442 

,307331") 73.307616 1 73.5074472 

5074703 75.3075345 73.307447P 

,3074037 73.507339P 75.3074401 

,5075046 75.50756B9 73.3074505 

,3075427 73.3076070 75.307452P 

,3076633 75.5077277 73.5074379 

.5076060 75.3076702 75.5074613 

,3073402 73.3076031 73.5074633 

,3076066 73. 3076709 7‘j. 5074666 

,5075329 75.3073071 75.5074620 

.3073627 73.5076330 75.5074702 

,5073413 75.5076038 73,5074717 

>3073710 73 3076353 73.5074732 

,3076900 75.3077351 75.5074722 

.3076827 75.5077470 75,3074c23 

,5076443 75.5077085 73.5074p35 

.3076224 75.5076867 75.5074Ppl 

.5076070 75.5076713 75.5074004 

.3076134 73.3076777 75.5074027 

.50740.07 75.3075550 75.5074026 

.3073926 75.5076369 73.3074044 

.5076305 75,5076046 ?5.5074O6<7 

.5076178 73.5076821 75.50749«9 

.5076129 75,5076772 75.307300? 

.5076243 73.5076886 75.5073029 

.3073023 75.3076566 75.5073043 
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TABLE A-3. EPHEI^IERIS AND LOP SUMMARY FOR TEST 2 
(NO SMOOTHING) (PAGE 1 OF 4) 


ATF-5 FPHF.[«FPIf; 

TIME 

GMT 

TTprrirnvM' ■ 

ipp0rci.0 

17^3000.0'' 227SG.35 

1«)0000.0 2276:2.11 


DEG DEG 

0.150 105.231 

105.230 

-0.075 105.?6;p 

-0.1P7 105. 2P2 

- - >P5.7>05 7^ A 

-0.395 105.303 


SAT PAN OE 
NMI 

227P 1 .0? 

" 2277F.0V5- 


L ATTTirnr • — Lmrmvvr 


PAPAFOLIC FIT TO FPHFrFP I F-LI NFAP TO CAL DATA 

CONSTANT LINEAR ClIAOPATIC PMS 

nATITODE ' 3. 1502327 -pr.7724KFT 0;00'42P 5 7" 0.00T05 

LONGITUDE 1 05 . 0’^6P2<i3 0.2165345 -0.0052P51 0.0002P 

GEOC RANGE 22676.52^7255 IR. 4257095 -0.7330221 0.03P<i4 

CALTP'UTrF P9660.2I'i6I6T --5-4T.377-^31 14.T24Q7 


GMT 

(FT ) 

191340.0 

100120.51 

'T9T350.6 

100117. 7T ■ 

10135 1 ./ 

1 00007.22 

101350.0 

100113.77 

1017^-4 71 

IOO10P.26 ■■ 

191355.3 

100115,74 

191401.5 

1001 15,74 

"loT^ry.p 

T'00 105.0 1 ■ 

|01 /06,2 

100.058 .60 

101407.3 

1 00050.83 

191409.3 

••"100040 ,p/ 

191410.1 

1 00100.24 

101410,0 

100039.02 

10141 1.6- 

T00037.06 

101413.2 

100102.96 

101414.0 

100000,14 

19141 4. T“ 

lOOl 16.73 

191415.5 

1 00128.53 

101416.7 

1001 n .2 1 

1014L7.O 

1000pp. 20 

101 / IP .6 

1001 10.82 

lOj/10.4 

1 00 0P 2.30 

101421 .0 

100005.09 

191422.5 

100124.50 

101423.3 

100117.71 

101424.1 

100120.66 

101406.0 

1 00 107.87 

101407.0 

1 00120.66 

101422.4 

1001 19.62 

101422.2 

100116.73 


POINT 1 
75. 50753/- <5 
75.'5074’61^3 
75.5072705 
75.50743«>0 
73'. 5074050 


CALWATT 
(FT) 

75130.^2 
70T5F.6^ 

7<5 130.56 
70130.33 

■ -75 I30 .T5 
7512«>.56 75.5074534 
7512<J.02 75.50.74525 

— T5r2R .6T '75 .5075^25 
70102,5! 75.5071030 
7«>12P,14 75.50705// 

■ '75T27 .23T ■75'.307F4O2 
7<>127.7I 75.5074224 
70127.59 75.5069230 
TOT2r;4^T5V5F6^n 2“ 
70107.24 75.5073277 
70107.12 75.5073202 

■ '7OlTrT;01 T5.3IT74T5?; 

7oj0<i.29 75.5075506 
75.5074463 


“TONGTTUDES (DEG) 


79106.71 
ToI0-6'.^3 
70106.42 
70106.30 

■ 79106V05 
79125.23 

70125.71 
70 1 25'.'3? ■ 
70125.20 
7ol05.il 
70124.93 
79124.27 


73.'3072O'p6‘ 
75.5074416 
75.50.72607 
75.5073/45 
75.5075316 
75.507/229 
T5.“50T502 2 

75.5074000 

75,5075107 

'75.5075056 

73.5074273 


P0INT2 
75.507599 1 
75'.*30730 61' 
75.50733/7 
75 . 5075030 
75. 50 74 T3 2 
75.5075176 
75.5075227 
75,5074627 
75.507 1 670 
75.5071 IC6 
75.5071140 
75.5074007 
75,5070472 
75.'507035/ 

75.5074510 

75.507384/ 

'75.5075399 

75.5076140 
75.5075105 
■T5.5'0 7362? 
75.5075059 

75,5073260 

75.5074022 
75.5075952 
75.5075532 
■■75.50 75724 
75.5074034 
75.5075750 
75.50756O'8 
75.5075515 


CUMUL. AVG. 

75.5075340 

75 .5074OC' A 
75.5074204 
75.5074065 
75.5074030 
75.50740t' 1 
75.507432/ 
75.52742^2 
75.5075001 

75.5073503 
75.50.73303 
75.50733P4 
75.50.731 1 1 
75.50700 60 
75.5070035 
75.5072052 
75 .5^7'305p 

75.5073104 

75.5073261 

75.5073247' 

75.5073303 

75.5073070 

75;507‘32o0“ 

75.5073365 

75.5073406 

'75.5075 400 
75.5073510 
75.5073576 
'75.5073627 

75.507366P 
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TABLE A~3. EPHEMERIS AND LOP SUMMARY FOR TEST 2 
(NO SMOOTHING) (PAGE 2 OF 4 ) 


Tlt^F 

PANOF DIF. 

CAL. 

Gt^'T 

(FT ) 

(F 

1 oi /?c>^ ^ 

100136.4? 

701 

IOU31 . 1 

100111.?! 

701 

JO] /i3 1 

1 00131 . 47 

70 1 

lol 43?.7 

1 00 110.6? 

70] 

101435.? 

1 00045 .o 1 

70 1 

lol/:36.6 

1 0O1P5.5? 

70| 

101450.7 

100004.10 

70 1 

10144?.5 

1 0 0 1 0 o . c- 4 

701 

10144 1.:^ 

1 00 1 ! P.70 

70 1 

1 Q 1 44?. r 

1001P7. 54 

701 

1 O 1 4 4 5 . P 

1001 P?. 53 

70 1 

10 1445. 1 

1 0OO<^6.P3 

70 1 

I O 1 4 /; s , o 

1 000^5.05 

70 1 

101446,7 

100 1 1 o , rr 

701 

101447.5 

100 1 34 . 4,^ 

70 1 

101 440. r 

1 00 IP7.54 

701 

101440.0 

100110.60 

70 1 

101450.6 

100 IP?. 53 

701 

101451.4 

100111.?! 

•70 1 

1 O 1 45P .5 

100117.71 

701 

101455.1 

1 0 0 1 0 0 , O C! 

70 1 

101456.0 

1P0143.P? 

701 

101456.? 

1001 P5.5^^ 

70 1 

101457,6 

1 00 1 44 , ?7 

70 1 

1 1 50 1 . o 

100053.70 

.79 1 

lol5P3.r^ 

1 00 1 0P.O6 

701 

101504.6 

1 000r 1 ,.T9 

70 1 

101505.^ 

i 0r0F5.?5 

791 

1 01506. P, 

1 00 1 07. ?7 

70 1 

101507.0 

1 00000.(1? 

70 1 

1 01 507.7 

1 00073.45 

70 1 

101500.9 

10013 1 .4P 

701 

1015] 1 ,6 

1 00070 .50 

70 1 

10151?. 4 

1000??. 30 

70 1 

1 O 1 5 1 3 . P 

1001 I 1 .? 1 

70 1 

101514.0 

100121,64 

70 1 

101517.1 

1001 1 1 .P 1 

70 1 

1 o 1 5 1 P . 6 

1 00 1 60. c 4 

70 1 

101510.4 

1 00006.07 

70 1 

1O15P0.6 

100121.64 

701 

1015??. 5 

100 1 10. ?2 

70 I 

1015P7.P 

1 00 1 0O,p4 

70 1 

1015??. 0 

100020.10 

70 1 

1 ol5?o,6 

100112.70 

70 1 

JO] 530.3 

100031 .15 

70 1 

101531 .1 

1001 10.6P 

701 

101531.0 

100065.5? 

791 


LOrCI TUDES (DEG) 

T) rOlMTl POINT? CUMUL, AVG, 

75.5f7^;i?l 75.507R7K^ 75. 

?^.5? 75.5D7^5P? 75.5?75??^ 

75,5D75P?F 75.5P7(^/i7P 
?4.?P 75.5D75{^'»1 75.5P75755 
P.5.FG 75.507047P1 75.507111? 

75.5P754Q5 75.507G157 
75.5P7555F 75.5074170 
75.5074556 75.507517^ 

7S.507/i 7PP 75.5075570 
PP.Pfi 75.5PI75GR3 75.50763'PIE 
75.507575^^ *75.507637? 

??,3^ 75.50730P7 75.50737PP 
PP.P7 75.507505? 75.5073674 
PP.15 75.5P75Pf7 75. 5075040 
PP.OP 75. 5076145 75.5076705 
PI.FP 75.5P757P1 75.5076364 
PI. 67 75.5075170 75. 5075^13 
PI. 55 75.5075707 75.5076450 
PI. 45 75.5074750 75.507530? 

PI.P6 75.5075151 75.5075773 

PO. 07 75.5074000 75.507474? 

PP. 73 75.5P76771 75.5P77413 
PO. 61 75.5P7566P 75.5076305 
P0.4O 75.5076846 75.50774FP 

0.F3 75.507 MP? 75.5O710P4 
0,55 75.5074?o5 75.507405? 
o.4p 75.507PO5O 75.50755? 1 
0.P4 75.5073106 75.5075?3P5 
0.1? 75.50746P4 75.5075P67 

O. 06 75.5074074 75.5074716 
F.05 75,5070460 75.5073111 
r.7P 75.5076157 75.5076770 

P. 36 75.507P315 75.507?Q57 
?.P4 75.5075065 75.50737P7 
?.l? 75.5074051 75.5075573 
7.o<> 75,5075556 75.5076109 
7.5P 75.5074Q65 75. 50756^5 
7.po 75,5074051 75.5075403 
7.17 75.507300? 75.5074633 

6.00 75.50756 11 75.5tT76?55' 

6.70 75.5074^^45 75.50755P7 

5.00 75.50740?? 75.5075565 

5.P6 75, 5 060? 4 0 75.50TO2O-r 
5.6? 75.50.751PP 75.5075770 75.5074166 

5.51 75.506OOO] 75.5070634 75.5074110 

5.50 75.5075575 75 . 5 0T6PTT~~ 7^75 0 7 4 1 30 
5.P7 75.5P7?173 75.507p?16 75.5074104 


75.5073775 
75.5073^36 
7'5 , 5T1T3FT5" 
75.5073775 
75.5073P25 
75 .■5073p'1'5“ 
75.5073P34 
75.5073P57 
75'.'5{W5502 
75.5073947 
75.5073927 
75V50V5O06 
75.5073035 
75.5073004 
‘ 75;50r?402Z 
75.5074047 
75.50740P3 
757TP74P^ 
75.5074 I 17 
75.50741 17 
75'. 5074 I 6? 
75.5074196 
75.5074245 
75-. 507 4 1 90 
75.507410? 
75.5074170 
75,5074153 
75.5074161 
75.5074159 
'T5.5074TT2“ 
75,5074164 
75.5074135 

7T. '5 074 IT? 1- 

75.5074130 
75,5074152 
7575074164 

75.5074174 
75.5074172 

757507TT9'Z J" 

75.5074205 
75.5074213 
'75.5074153 
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TABLE 


TI 1 «E 

c,m 

0153?. 7 
01555. A 
01 53fi,fi 

01557.7 
0155? .5 

01530.5 

01550.7 
01541.2 
0 1 543 .P 

01545.1 
0154f^.7 
0154P.7 

01552.5 

01554.0 

01555 .5 
0155?..'^ 

01550.2 

01550.0 
;Olfif1p.7 
;oifi?3.«^ 

|Oi f,?/i 

101505. F 
OlfiPf;.? 
01 fi07.4 
oifipo.7 
01^11.5 

0 1 r. 1 5 . 

01514.0 

01515.1 

01 f.lf..5 
ol«^17.l 

oUl 7.0 
Ol<^|O.0 
01F21 .0 
0W^21 .P 
01fi22.f: 
Ol fi24.o 
01 f.2P .2 
oi62o.r. 
0 1 f.30.3 
0If;51 . 1 
01652.7 

01655.5 
0 1634.2 
Ol 655 ,P 

01656.6 
01657.4 


A-3. EPHEMERIS AND LOP SUMMARY FOR TEST 2 
(NO SMOOTHING) (PAGE 3 OF 4) 


PAMGE m F. 
(FT ) 

1 00 12P.6(^ 

- I001IP.GO 

1001 I 1 .PI 

1001 IP .60 

100105.04 

1 00007.05 
100076.40 
1 000QP .0j1 
100 I IP. 60 
1 00 I 0P .P6 
100104 .0?> 

100102.06 
1 P002P.20 
1 00054.76 

1000OO .0/r 
100105.01 
100112.60 
100047. PP 
1 00066.56 
100075.42 
I 00 152.46 

100125.6 1 
lOOPPO. lo 

100067.55 
1 00150.50 
1 00 1 OF .P6 
1 0O17O.F4 

100122.65 
100150.50 

100100.00 

100103.04 
17-2131 .40 
1 001 0o,p4 
102-132,46 
1001 17.71 

100045.01 
1P010F.P6 
1001 16.73 

100123.61 

120103.04 

1001 11. PI 

100126.56 

100104.02 
I00125.5P 

100122.66 
100102.o'6 
1001 10. P2 


C/^L. VAL. LOPCITLfPEF(DFn 

(FT) POINTl POirTP CUr^.llL. AV/G. 
70115.15 75.5075640 75.5076202 75.5074124 
70114.74 75.5075540 75.5076|o0 75.527414? 

70114.56 75.50.75104 75.5275767 75.5274154 

70114.30 75.5275567 75.5076200 75.5274172 
70114.27 75.5074644 75.50752P7 75.5074(7^ 
70114. 1 5 75. '5074217 75.5074^59 75,5074I7<’ 
701 14. OP 75.5072O1P 75.5073562 75.5274163 
70115.26 75. 5274005 75.5074037 75.5074165 

70115.40 75.5275616 75.507625^ 75.50741^1 
70115.26 75.5275000 75.50-75651 75.527.'MO] 

70115.02 75.5074774 75.5075416 75.5074ioa 
70112.72 75.5074667 75.507550P 75.5074225 
70112.14 75.506OQpo 75,50-70631 75.5274156 
70 111.77 75.52716^^2 75.527'^524 75.52741f^o 
70111.71 75.5274/11? 75.5275l'54 < 75.5074132 

70111.30 75.507403'! 75.5075572 75.5274141 
70111.12 75.5075745 75.50763P7 75.507415^' 

70111.01 75,5071200 75.507103? 75.5074127 
701 10. FP 75.5272473 75.5073116 75.50-74110 
70110.4? 75.5075057 75.50-73600 75.5r?4roo 

70110.30 75.5076657 75.5077200 75.5074125 
70110.1? 75,5076100 75.507675? 7 3 . 5 2-7 4 1 .'o 
70110.05 75,52.75044 75.50745F7 75,5074143 
7P10P.87 75.5072501 75.5P73233 75.50741?p 
70100,52 75.5076576 75,507721? 75.5074152 
70100.23 75.527522P 75.5075^70 75 . 527 4 1 6:’ 

70100.03 75,50 7 530 6 75,52 7 50 40 7 5.527.^117,3 
7O10F.F7 75.50761 15 75.507675^ 75.527410? 
7O10F.7O 75.507662.0 75.5077265 75.5274215 
7O10P.52 75.50.7477 1 75.50754 14 75.5274220 

70102.40 75.5074064 75.5075606 75.5074227 
7O10P.2F 75.50-76705 75.5077347 75.5074252 
7ol0F.ll 75.5275551 75.5275oo3 75,5074260 
70107, P0 75.5076705 75.5077435 75.5274?o3 
70107. 6« 75.5075P70 75.5076512 75.5074907 

70107.56 75.5071354 75. 5071006 75.5074271 
70107.21 75.507535P 75.5075PP1 75.50742c-0 
70 126.71 75.5075P61 75.5276504 75.5274?o/. 
70106,50 75.5076306 75.5076040 75.5074311 

70126.30 75.5075073 75.5075715 75.5074310 
701 06. 2.7' 75.5075575 75.507621P 75.50743?c^ 

70106.02 75.5076518 75.50-77 160 75.5074547 
70105.00 75.5075161 75.5075P04 75.5074353 
70105. 80 75.50'7646P 75.5077111 75.507/.371 
70105.55 75.5276172 75.5076^14 75. 50743^6 
70105.43 75.50750-63 75.507570(^ 75.5074301 

70105.31 75.5075565 75.5076207 75.5074401 
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TABLE A“3 


EPHEKERIS AND LOP SUMMARY FOR TEST 2 
(NO SMOOTHING) (PAGE 4 OF 4) 


TIFF 

PAMGE DIF. 

CAL. V'AL. 

GFT 

(FT ) 

(FT) 

1^155?. I 

IPPI P7.5/: 

70 1 P5 .^0 

10 )^ 50.7 

lf’(M 55.^1 

7«s 1 O/i.Qf. 

1^1 /ip. 5 

imP5.n j 

7^ 1 O/ 1 .C 4 

1^1 F4 1 .5 

] Pf^I 1 1 .P- 1 

1 OA . 77 

jolfi/p.r 

IPOlPP.Pr- 

7Q UM.G 1 


1 POMP. 33 

7Q10A.37 

|C: 1 ,yi 

1001 15.7/1 

7Q 1 OA, 0/ 

IoK>Z!5.o 

100113.77 

7 0 10/1,00 

101^5 I 

1001 p/i. 5-^ 

7-3 103. !<’ 

lQlfi53.7 

rooioo.oo 

7Q1 00.P3 


LONGiTnrtFf (HFr.) 

PoiMi roirT? cuMUL. a\t.. 
7‘s 75 , !^F77pr.r. 75 

75.5r-77 1 75.5:'77775 7 5,5r7/,/;/;^- 

7'’i.5r7S?c^ 1 75,5r'75^p/- 7 5 . 5P7 / 7 

75. 5F75'^.A;r 75.5:’7A5PP 75.5P74^5f^ 

75.5P7^o?o 75.5P755G/J 75.5P?i/;f;r 

75.5P77/*75 75.5P7c'|jr’ 75.577/^':'’. 

7‘S.5P'75'^35 75.577^575 75 . 5 77 / 

7 5.5P75PPJ 7 5. 57 7 6^/.. 5 7 5.577/^'^:’/ 
75,577r.5/c’ 75.5777 I'M 75 , 5/7 -M' ] ^ 
75.5P75P|= 75.5r75Gf.l 75,5r745P.5 


S 
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TABLE A~4. EPHEMERIS AND LOP SUMMARY FOR TEST 5 
(JULY 15, 1971) (PAGE 1 OF 3) 


ATS-5 EPHEMERIS 

TIME SAT RANGE LATITUDE LONGITUDE 
GMT NMI DEG DEG 

173000. H 22773.87 -0.118 105.284 

180000.0 22770.00 -0.228 105.2<)R 

185000.0 227G5.81 -0.334 105.505 

190000.0 22761.36 -0.455 105.318 

193000.0 22756.74 -0.528 105.323 


PARABOLIC FIT TO EPKEMERI S-LI NEAR TO CAL DATA 

CONSTANT LINEAR QUADRATIC RMS 
LATITUDE 7.2793723 -0.6176P58 -0.0111429 0.0005^ 
LONGITUDE 102.9904324 0.2310838 -0.0057158 0.00021 
GF.OC RANGE 22752.3074475 10.0394174 -0.5032261 0 
CALIB LINE -38074.8373629 -618.8420590 18.43029 


TIME RANGE DIF. CAL, VAL. LON GI TUDES (DEG) 

GMT (FT) - (FT) POI NT I POINT2CUMUL. AVG. 

170214.9 -27762.06 -48618.34 75.5076528 75.5076964 75.5076528 

170239.5 -27828.90 -48622.57 75.5072346 75.5072982 75.5074337 

170306.0 -27822.57 -48627.13 75.5072992 75.5073628 75.5073PP9 

170334.4 -27791.00 -48632.01 75.5075247,75.5075883 75.507422? 

170346.1 -27792.65 -48654.02 75.5075252 75.5075888 75.5074433 

170555.9 -27802.27 ,-48655.36 75.5074718 75.5075354 75.50744P1 

170426.5 -27842.61 -48640.93 75.5072479 75.5075115 75.5074195 

170500.6 -27814.96 -48646.83 75.5074542 75.5075178 75.507423? 

170527.5 -27826.65 -48651.45 75.5O74056 7^.5074692 75.507421P 

170552.0 -27838.94 -48655.66 75.5073511 75.5074147 75.5074147 

170613.1 -27863.97 -48659.29 75.5072131 75.5072767 75.5073964 

170630.6 -27841.00 -48662.30 75.5073741 75.507437? 75.5073945 

170719.0 -27859.21 -48670.62 75.5075046 75.5073682 75.5073876 

170754.8 -27851.35 -48676.77 75.5073876 75.5074512 75.5073876 

170824.9 -27887.76 -48681.94 75.5071863 75.5072500 75.5073742 

170843.9 -27821.46 -48685.21 75.5076218 75.5076854 75.5073897 

170901.5 -27892.82 -48688.24 75.5071887 75.5072523 75.507377P 

170926.4 -27864.75 -48692.52 75.5073888 75.5074524 75.50737^4 

170948.6 -27905.50 -48696.33 75.507I52P 75.5072164 75.5073666 

171025.3 -27906.90 -48702.64 75. 5071783 75.5072419 75.5073572 

171054.1 -27888. R”? -48707.59 75.5073182 75,5073818 75.5073553 

171128.4 -27890.24 -48713.49 75.5073423 75.5074059 75.5073547 

171212.1 -27917.56 -48721.0075.507211075.5872746 75.5075485 

171238.6 -27892.69 -48725.56 75.5073925 75.5074561 75.50735P3 

171315.6 -27913,16 -48731,9275.5072981 75.5073618 75.50734p2 

171550.3 -27855.54 -48737.88 75.5077062 75.5077699 75.5073620 

171432.0 -27889.97 -48745.05 75.5075157 75.5075794 75.5073677 

171449.2 -27915,41 -48748.01 75.5073715 75.5074352 75.507367? 

171521.5 -27094.50 -48755.56 75.507534^ 75.5075^84 75.5075756 

171542.2 -27905. /'7 -48757. !2 75.5074964 75.5075601 75.5075777 

171554.5 - •--.'cip/: ^ 75 -4875^.20 75,5075610 75.507424? 75.5075771 
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TABLE A“4. EPHEMERIS AND LOP SUMMARY FOR TEST 5 
(JULY 15, 1971) (PAGE 2 OF 3) 


TIME 

RANGE DIF. 

CAL. VAL. 

LONGITUDES(DEG) 

GMT 

(FT) 

(FT) 

POINTl 

P0INT2 

CUMUL. AVG. 

1 71624. J 

-27945.03 

-48764.35 

75.5072866 

75.5073502 

75.5073743 

171647.5 

-27927.21 

-48768.31 

75.5074078 

75.5074714 

75,5073753 

J 71 714.6 

-27943 . I 5 

-48773.00 

75.5073330 

75.5073967 

75.507374 1 

171757,9 

-27936.01 

-48780,44 

75.5074185 

75.5074821 

75.5073753 

1 7 1 8 2B . 6 

-27939.1 1 

-48785. 72 

75.5074277 

75.50749 1 4 

75.507376R 

171843.8 

-27980.30 

-48788.35 

75.5071824 

75.5072461 

75,5073715 

1719R8.8 

-27928.45 

-48792.63 

75.5075325 

75,5075962 

75.5073758 

171931,4 

-27973,54 

-48796.52 

75.5072709 

75.5073346 

75,5073731 

171945.0 

-27969. IS 

-48798.85 

75.5073100 

75.5073-736 

75.5073715 

1 /2U03, f 

-27951.78 

-48802,07 

75.5074370 

75.5075007 

75.5073731 

1 72030.7 

-28009.71 

-48806.71 

75.5070973 

75.5071610 

75.5073665 

1 72057,5 

-27955.28 

-488 1 1 .32 

75.5074654 

75.507529 1 

75.5073688 

1 17A22, 1 

-27953.09 

-48815.55 

75.5075022 

75.5075659 

75.5073719 

1 72159,1 

-27965.75 

-48821 .91 

75.5074572 

75.5075209 

75.507373? 

172225.3 

-27976.54 

-48826.41 

75.5074137 

75.5074774 

75,5073746 

1 72259,9 

-27945.25 

-48832.36 

75.5076433 

75.5077070 

75.5073803 

172319.0 

-27977,24 

-48835.64 

75.5074597 

75,5075234 

75.5073820 

172334.2 

-27971.21 

-48838.25 

75.50751 19 

75.5075756 

75.5073846 

172424.1 

-27954.23 

-48846. «3 

75.5076657 

75,5077294 

75,5073903 

172431.9 

-28015.73 

-48848.17 

75,5072982 

75.5073619 

75.5073885 

172458.4 

-28010.47 

-48852.75 

75.5073436 

75.5074073 

75.5073876 

1 72532.4 

-27987.95 

-48858.57 

75.5075174 

75.5075811 

75.5073901 

1 72607.0 

-28024.77 

-48864.52 

75.5075179 

75.5073816 

75.5073887 

1 72622.6 

-28006.07 

-48867.20 

75.5074504 

75.5075141 

75.5073898 

1 72638.2 

-28004.37 

-48869.88 

75.5074757 

75.5075394 

75.50759 1 4 

1 72700.0 

-27971 . 1 1 

-48873.63 

75.5077057 

75.5077694 

75.5073969 

172713.7 

-2804 1 .90 

-48875.99 

75.5072726 

75,5073363 

75.5073947 

172725.0 

-28050.05 

-48877.93 

75.5072319 

75.5072956 

75.5073920 

172744.5 

-28018.88 

-48881.28 

75.5074465 

75.5075105 

75.5073929 

172801 .6 

-27985.67 

-48884.22 

75.5076718 

75.5077555 

75,5073975 

172822.7 

-28059.27 

-48887.85 

75.5072280 

75.5072917 

75.5073947 

172835. Q 

-2806 1.06 

-48890.12 

75.5072291 

75.50729 28 

75.5073921 

1 72851.5 

-28000. 16 

-48892.80 

75.5076274 

75.50769 1 1 

75.5073958 

172912.2 

-28034.84 

-48896.36 

75.5074284 

75.5074921 

75.5073963 

172927.0 

-28012.09 

-48898.90 

75.5075856 

75.5076495 

75.5073991 

172953.9 

-28057.32 

-48905.52 

75.5073259 

75.5073897 

75.5073980 

173014.2 

-28007.87 

-48907.01 

75.5076565 

75.5077202 

75.5074010 

173052.9 

-28008.65 

-48910.23 

75.5076692 

75.5077329 

75.5074037 

173101.7 

-28042.49 

-48915.18 

75.5074831 

75,5075468 

75,5074068 

1 75123.9 

-28046, 1 4 

-48919.00 

75,5074809 

75.5075447 

75.5074079 

173141.5 

-28047.89 

-48922.02 

75.5074865 

75.5075502 

75.5074090 

173204. 1 

-27997.51 

-48925.91 

75.5078251 

75.5078888 

75,5074147 

173225.6 

-28045.87 

-48929.26 

75.5075387 

75.5076025 

75.5074163 

1 73238.4 

-28082. 16 

-48931 .80 

75.5073241 

75.5073878 

75.5074151 

173259.4 

-28023.75 

-48935.41 

75.507711 7 

75.5077755 

75.5074190 

173321 ,3 

-28056.09 

-48939.18 

75.5075286 

75.5075924 

75.5074204 

173548. 

-28053.87 

-48945.92 

75.5075686 

75.5076323 

75,5074223 
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EPHEMERIS AND LOP SUMMARY FOR TEST 5 
(JULY 15, 1971) (PAGE 3 OF 3) 


RANGE DIF. CAL. VAL. LONGI TUDES (DEG) 

(pj) (FT) POINTl P0INT2 CUMUL. AVG. 
-2807H A! -AP9A7.54 75.5074831 75.5075AGP 75.5074231 
-28107. 0Q -48<)4<5.5G 75.5072642, 75.5073279 75.5074211 
-28019.12 -48952.31 75.5078333 75.5078971 75.5074262 


L2 



TABLE A- 5, 


EPHEMERIS AND LOP SUMMARY FOR TEST 8 
(JULY 21, 1971) 


ATf -3 
TIMF 

FPPFMF.FI? 
OAT PAMCF 

LATI TUPF 

LOm TliPF 

GMT 

8 Ml 

PEG 

PFG 

70000.8 

02743 ,43 

0. 40J 

1 04 , 06*3 

73000.0 

007SO, \o. 

0.577 

1 04.® 66 

POOOO . 0 

2075 4 , o<^ 

0.63 1 

104.067 

P 3 00 0.0 

2275® .73 

0.7 1 4 

1 04.® 70 

oOOOO .0 

22764,3® 

0.764 

104.076 


PAPAPOLIC FIT TO F.PHF.t^^FPl P-Ll N’FAP TO CPL DATA 


LATI TUFF 
LOF'filTUr-F 
GFor PPFOF 
CAL-ir LI f'F 


CONSTAt'T 
“ I .P70PP5G 
lP-5.3(*5GGS/i 


LI MEAP 

-P*0P7P307 
I 0. 77G 1 3G I 


GPOIP .FP.77 1 PP “G7P ,?5^4‘:i5Afi 


nUADPATIC PM? 

“0\0PPri7|i5 0.0^0/13 
?!. 00371^4 0 . 00(^21 
-r.rwOOOFA 0,022G3 
73. 13237 


PAMOF DTF. CAL. VAL. LOrOI TIIDE? (PEG) 

(FT) , (FT) POIFTl P0INT2 CLIMUL. AVG, 

F51P/1.F3 ”g 23^A.13 75 .307/ i ! i3 73.307/747 73.307411/ 
FJc>f^0,F3 G2325.41 73,3073^^5.’ 73.3C7AG31 73.3073036 

r3162.'il 62373.33 73.3P74J31 73.30747P3 73.3074734 

F3173.21 62532.3® 73.3P764F7 73.30771!9 73.30731^7 

P3167.72 62342. ®3 75.3076662 73.30772^3 73.30734®? 

®3134.22 62333.42 75.3075173 75.3P73P0F 73.5075431 

P31^3.36 ■ 62321, 3'3 73.3C7341P 75.3076051 73.307342® 

®30®1.33 62301,22 73.3074644 73.3073276 73.5073331 

Pv''-0®7,70 624pp. PI 75.5273®®3 73.3076523 73i30733?3 

F5074.7P 62470,06 75.3073712 73.3076544 75.3075423 

P3021.32 62440,63 75.5074313 73.3P7494P 75.3073324 

P5022.I3 62424.63 75.3075433 73.30760F7 73.3075335 

Pr®2®.7F 62360. P2 73.3073®16 73.307434p 73.3073226 

P2®22.73 62343,6! 73.30744op 73.3073130 73,3073174 

F^®0®.63 6P334.44 75.5074422 73.3073034 75.5075124 

F2P73.2P 62314.38 75.3073471 73.5074J03 73.5075021 

®2FP3.60 62303.36 75.3074®P4 73.3073616 73.307301F 

®2®P4.14 62?®3.33 73.30753®! 75.3076213 75.5075030 

P2F37.63 622P4.F1 73.3073206 75.3073F37 73.3074®33 

F2P30.PO 62267.®! 73.50731F3 73.3073P14 73.5074Q64 

P2P17.P1 62236. P9 73.5073P32 73.3074464 73.5074®10 

F2P14.34 62257,®? 73,5074004 73.3073336 73.3074010 

F2766.36 62225.73 73.3072680 73.3073312 73,3074813 

P274P.75 62196.00 75.5075400 73.3074122 75.507475P 

F273F.M 6217P.33 73.3074063 73.3074607 75.3074730 

09743.61 62166.17 75.3075233 75.5075P63 75.3074730 

F270'J.74 62132.70 75.3073083 75.5074614 75.5074721 

F260B.21 62132.52 75.5074613 75.5075245 73.5074717 

P260l,22 62120.31 73.5074QQ0 73,5073621 75.5074727 

p9fi63.2P 62106.71 75.5074125 73.5074756 73.5074707 

•P2630.®! 62087,23 73,5074643 73,50732.76 73.3074703 

P264p,®4 62072. F3 73. 50754P6 75,3076117 73,307472® 



TABLE A-6. EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(JULY 27, 1971) (PAGE 1 OF 3) 


ATS-5 

EPHEMERIS 





TIME 

SAT RANGE 

LATITUDE 

LONGITUDE 



GMT 

NMI 

DEG 

DEG 



0.0 

22717.71 

-0.722 

105.216 



3000.0 

227 16.6B 

-0.664 

105.193 



10000.0 

22716.27 

-0,595 

105.169 



13000.0 

227 16.4B 

-0.515 

105.146 



20000.0 

22717.29 

-0.426 

105.123 



23000.0 

227 1R.70 

-0,330 

105.101 



PARABOLIC 

FIT TO EPHEMERIS-LINEAR TO CAL DATA 



CONSTANT 

LINEAR QUADRATIC 

RMS 

LATITUDE 

-0.7226429 0.1091571 0 

.0192857 0, 

0008 1 

LONGITUDE 

105.2162500 -0.0477786 0 

.0006429 0. 

00026 

GEOC RANGE 

22717.7054555 -2.6502561 1 

.2198813 0. 

00495 

CALIB LINE 

-79619.4784223 -605.7383593 

10. 

44717 

TIME 

RANGE DI F. 

CAL. VAL. 

LONGITUDES (DEG) 

GMT 

(FT) 

(FT) 

POINT 1 

P0INT2 

CUMUL. AVG. 

131 1S.0 

-59642.20 

-80611 .54 

75.5072061 

75.5072708 

75.5072061 

13212.2 

-59650.60 

-80621.00 

75.5072187 

75.5072835 

75,5072124 

13306.1 

-59666.42 

-80630.07 

75.5071818 

75,5072465 

75.5072022 

1340A.6 

-59675.39 

-80659.91 

75.5071936 

75.5072583 

75,5072000 

13505.6 

-5968 1.24 

-80650,18 

75.507228 1 

75.5072929 

75,5072057 

13603.4 

-59683.27 

-80659.90 

75.5072851 

75.5073478 

75.5072186 

15717.5 

-59700.56 

-80672.37 

75.5872607 

75.5075254 

75.5072246 

13R33.3 

-59714.00 

-80685. 12 

75.5072646 

75.5073293 

75.5072296 

141 1 1 .A 

-59742.44 

-80711 .79 

75.5072709 

75.5075356 

75.5072542 

14259.5 

-59766.58 

-80729.91 

75.5072460 

75N5073107 

75.5072354 

145 13. K 

-59792.66 

-80752.51 

75.5072378 

75.5073025 

75.5072356 

15243.6 

-59863.00 

-80828.20 

75.5073229 

75.5075876 

75.5072429 

15409.4 

-59859,64 

-80842.63 

75.5074456 

75.5075102 

75,5072584 

1 5547. g 

-59862.47 

-80859. 19 

75.5075440 

75.5076086 

75.5072788 

15716.1 

-59876.87 

-80874.05 

75.5075572 

75.5076219 

75.5072974 

1 5049,0 

-59893.44 

-80889.68 

75.5075622 

75.5076269 

75,5073140 

20031 .3 

-59918.06 

-80906.89 

75.5075275 

75.5075921 

75.5073265 

20157. 1 

-59931.67 

-80921 ,33 

75.5075429 

75.5076076 

75.5073385 

20339.4 

-59949.95 

-80938.54 

75.5075485 

75.5076151 

75.5073496 

20519.4 

-5998 1 .96 

-80955.37 

75.5874644 

75.5075290 

75.5073553 

20657.0 

-60009.51 

-80971 ,92 

75.5074067 

75.5074713 

75.507357? 

20R22.9 

-60023.82 

-80986,24 

75,5074171 

75.50748 18 

75.5073605 

20955.0 

-60021.60 

-8 1001.74 

75.5075406 

75.5076053 

75.5073685 

21 1 16.2 

-60011 .56 

-81015.40 

75.5077008 

75.5077654 

75.50735722 

21236.6 

-60029.94 

-81028.93 

75.5076800 

75.5077446 

75.5073941 

21402.5 

-60064.56 

-81045,38 

75.5075630 

75.5076276 

75.5074006 

21509,7 

-60077,07 

-81054.69 

75.5075638 

75.5076283 

75.5074066 

21616. g 

-60100.76 

-81065.98 

75,5074936 

75.5075582 

75,5074097 

21721 .6 

-60104.06 

-8 1076.88 

75.5075500 

75.5076145 

75.5074146 

21P4 1.2 

-60124.15 

-81090.28 

75,5075176 

75.50758 22 

75.5074180 
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TABLE A“6. EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(JULY 21 , 1971) (PAGE 2 OF 3) 


TIME 

RANGE DI F. 

CAL. VAL. 

LONGITUDES (DEG) 


GMT 

(FT ) 

(FT) 

POl NT 1 

POI NT2 CUMUL , AVG. 

21957.5 

-60134.45 

-81099.75 

75.5075)95 

75.5075841 

75.507421 3 

2204 U5 

-60147.31 

-8111 0.52 

75.5075144 

75,5075790 

75.5074242 

22143.2 

-60154.78 

-8 1 120.90 

75.5075408 

75.5076053 

75.5074277 

22250,3 

-60170.23 

-8 1 132. 19 

75.5075230 

75,50758 76 

75.507^305 

22356.7 

-60184. 12 

-8 1 1 43.36 

75.5075143 

75.50757?^ 

75.507432'i 

22510.9 

-60201 .22 

-8 1 155.85 

75.5074947 

75.5075592 

75.5074346 

22618.8 

-60217,97 

-8 1 167.27 

75.5074698 

75.5075343 

75.5074356 

22725.2 

-60236.79 

-81 178.45 

75.5074300 

75,5074945 

75.5074354 

22830.0 

-60265.85 

-81 189,35 

75.5073234 

75.50758 79 

75.5074326 

22^27.7 

-60280.79 

-P 1 199.06 

75.5072978 

75,5073624 

75.50742‘32 

23031 .8 

-60296.84 

-8 1209 .84 

75.507272'^ 

75.5075374 

75.5074254 

23138. 1 

-60308.71 

-3 1 221.00 

75.5072771 

75,5073417 

75.50742l« 

23252.3 

-60320.93 

-8 1233.49 

75.5072943 

75.5073588 

75,5074189 

23401 .8 

-60329.36 

-81 245, 18 

75.5073 185 

75.507585 1 

75.5074166 

23512.9 

-60360.65 

-81257.14 

75,5072056 

75.5072701 

75.5074 1 1<5 

23608.3 

-60366.57 

-8 1 266,46 

75.5072346 

75.5072991 

75.507408 1 

25712.3 

-60375.41 

-81277.23 

75,5072554 

75.5073 199 

75.5074040 

25807. 0 

-60389.99 

-81286.44 

75,5072287 

75.5072932 

75.507401 1 

23902.4 

-60394.53 

-81295.76 

75.5072665 

75.5073310 

75. 5073884 

24056.4 

-60418.24 

-81314.94 

75.5072532 

75.5075177 

75.5073955 

24202.8 

-60432.61 

-81326.1 1 

75.5072420 

75.5073065 

75,5073925 

24306.0 

-60442. 14 

-81336.75 

75.5072577 

75.5073222 

75.5073P8O 

24406.2 

-60452.36 

-8 1346.88 

75,5072654 

75.5073298 

75.5073875 

24505.5 

-60460.4 1 

-8 1356.85 

75,5072857 

75.5073502 

75.5073856 

24608.0 

-60474.80 

-81367.37 

75.5072698 

75.5073343 

75.5073835 

24 708. 1 

-60485.99 

-81377.48 

75.5072713 

75.5073358 

75.5073815 

24812.9 

-60495.07 

-81388.39 

75.5072919 

75.5073564 

75.5073800 

24^20.8 

-60510.58 

-8 1399.81 

75.5072755 

75.5073399 

75.5073782 

25043.6 

-60524.97 

-81413.74 

75.5072842 

75.5073486 

75.5073766 

25 144.5 

-60535.48 

-81425.99 

75.5072^1 1 

75.5073555 

75.5073751 

25246.9 

-60529.93 

-81434.49 

75.5074016 

75,5074660 

75.5073756 

2535 1 .8 

-60542.53 

-81445.41 

75.5074001 

75.5074645 

75.5073760 

25448.0 

-60552.39 

-8 1 454.87 

75.5074055 

75.5074699 

75.5073764 

25538.7 

-60559.96 

-8 1463.40 

75.5074 188 

75.5074853 

75.5073771 

25634 .<J 

-60567.81 

-8 1472.85 

75.5074570 

75.5075014 

75.5073780 

25726.5 

-60577.96 

-8 148 1 .53 

75.507435 1 

75.5074^^95 

75.50737OQ 

258 16.4 

-60584.61 

-81489.93 

75.5074533 

75.5075177 

75,5073800 

2591 1 .9 

-6059J>.25 

-8 1 499.27 

75.5074721 

75.5075565 

75.5073814 

30007.5 

-60599. 15 

-81508.59 

75.5074954 

75,5075598 

75.5073830 

30059.6 

-60609 .85 

-81517.39 

75,5074909 

75.5075555 

75.5075845 

30155.1 

-60627.69 

-81526.73 

75.5074451 

75.5075095 

75.5073854 

50246.6 

-60645,96 

-81535.40 

75.5074044 

75.5074688 

75.5073857 

30338 . 9 

-60648,54 

-81544.20 

75.5074388 

75.507503 1 

75.5073864 

30428 . 1 

-60664.71 

-81552.47 

75.5073959 

75.5074603 

75.5073865 

30552. 1 

-60685.13 

-81563.24 

75.5073441 

75.5074085 

75.5073860 

30650.7 

- 60695.99 

-8 1 575. 10 

75.507359 1 

75,5074234 

75.5073856 

30750.8 

-60708.32 

-8 1583.21 

75.5073412 

75.5074056 

75.5075850 
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TABLE A“6 . EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(JULY 27, 1971) (PAGE 3 OF 3) 


ME P/iNGE DIF. 


Gf^T 
imA I .P 

I 

31 15R.3 
5 123P.4 
31337. P 
3U4 1 .P 
3I53P.P 
31G35.P 
31 73P.3 
3 IP32.<J 
3 

3P.PI24 ,6 
32117.7 
3221P.pl 
323P3. I 
32354. G 
33'^!^. 1 
340! 1 .4 
34102. 1 
34 I 55.2 

34254.6 

34348.5 
34436. <5 

34525.3 

34614.5 

34709.9 

34759.9 

34849.8 

34«>39.P 

35026.7 

351 19.8 

35212.8 

35302.8 
3535^.8 

35451.4 

35544.4 

35640.7 

35737.7 

35831.5 

35924.6 

40038.8 

40146.7 
40304.0 


(FT) 

-60722.66 

-60735.41 

-60753.31 

-60761 .64 

-60773.54 

-60787.32 

-607<»7.80 

-60809.75 

-60820.34 

-60830.95 

-60841 . 17 

-60847.85 

-60863.06 

-60867.46 

-60875.04 

-60P86.64 

-61025.50 

-61032.57 

-61036.05 

-61056.32 

-61071 .03 

-61088.35 

-61092.65 

-6 1 102.51 

-61111.40 

-6! 122.65 

-61 133.79 

-61141.69 

-61151.71 

-61 155.87 

-61 171 .50 

-61 184.33 

-61195.10 

-61214.25 

-61219.11 

-6) 228.46 

-61241 .26 

-6124 1.07 

-61250,09 

-61244. <57 

-61265.04 

-61252.41 

-61250.87 


CAL. VAL. 

(FT) 

-81595. 18 

-81604.63 

-8 1624.86 

-81634.97 

-8 164.4.97 

-81655.74 

-81665.35 

-81674,92 

-8 1685.43 

-8 169^.62 

-81703.42 

-81713.41 

-81722.35 

-81731.15 

-81740.08 

-8 1748 .75 

-8 1904.31 

-8 1913.1 I 

-81921.64 

-81930.57 

-81940.57 

-8 194<5 .64 

-8 1«>57.78 

-81965.92 

-81974.20 

-8 1983.52 

-R 1 .<>4 


LONGITUDES (DEG) 


POINT 1 
75.5073366 
75.5073241 
75.5073567 
75.5073770 
75.5075758 
75.5073644 
75,5073673 
75.5073610 
75.5073700 
75.5075693 
75.5073683 
75.5073<^84 
75.5073667 
75.5074026 
75.5074194 
75 ,5074088 
75.5076607 
75.507680] 
75.5077203 
75.5076571 
75.5076368 
75.5075935 
75.5076256 
75.5076226 
75.5076268 
75.5076237 
7S 1 47 


POI NT2 
75.5074010 
75.5073884 
75.507421 1 
75.5074414 
75.5074382 
75.5074287 
75.5074317 
75.5074253 
75.5074343 
75.5074336 
75.5074326 
75.5074627 
75.5074310 
75.5074669 
75.5074837 
75.5074731 
75.507724<) 
75.5077443 
75.5077846 
75.5077213 
75,5077010 
75.5076575 
75.5076898 
75.5076«6P 
75.50769 10 
75.5076879 
75-50767P8 


CUMUL. AVG. 
75.5073844 
75.5073836 
75.5073833 
75 . 5073<732 
75,5073831 
75.5073p2<^ 
75.5073827 
75.5073^24 
75.5075825 
75.5075821 
75,5073<?20 
75 . 5075<722 
75.5073820 
75.5073022 
75.5073026 

75.5073020 

75.5073P5O 

75.5073090 

75.5073924 

75.5073051 

75.5073076 

75.5075096 

75.507401P 

75.5074040 

75.5074062 

75.5074003 

7*i-6074103 


-82000.35 
-82008.75 
-82016.64 
-82025.57 
-82034.49 
-82042.90 
-82052.40 
-8206 1 ♦ 18 
-82070.00 


75.5076261 

75.5076242 

75.5076556 

75,5076220 

75.5076060 

75.5075094 

75.5075485 

75.50758 12 

75.50758 74 

75.5075758 


75.5076002 

75,5076883 

75.5077190 

75.5076861 

75.5076702 

75.5076636 

75.5076125 

75.5076454 

75.5076515 

75.5076300 

75.5077116 

75.5077209 


75.5074124 
75.5074144 
75.5074166 
75.5074185 
75.5074202 
75.5074219 
75.5074230 
75.5074244 
75.5074250 
75.5074272 
75.5074291 
75.507431 1 
75.5074330 
75.5074362 
75.5074400 
75.5074446 


-82079.57 

-82089.16 75,5076475 
-82008.21 75.5076568 
-82107.14 75,5077550 75.5078191 
-82119.63 75.5077105 75,5077836 
-82131.05 75.5078837 75.5079478 
-82144.06 75.5070892 75.5080535 
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TABLE A- 7, 


EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(NO SMOOTHING) (PAGE 1 OF 3) 


fTr*‘=. 
TI r’F 
GMT 
r.G 
.1' 

I r» /■" 


f. puf r I f 

GAT RAf'GF 
f'Ml 
17 

Gr7l^> 

TR7 1 r. 

'^?1 I 7 


71 

r‘; 

T7 

/I” 

, 

. 7G 


LATi Tunr 

riG 

-(> , S^‘=i 

- G . 1 3 

-r. 


LO'T.l TITF 

nr.G 

1 

I (’ 3 

ir". 
irs 
1G3 


I '=5 

l.rp 

1 / f 

j 

iV I 


r’APAPOLI r. 

FIT TO FFlfFOrPl? 

-LKTAP TO CAL 

J 

1 

(1^ 


COOPT AM 

LI PFAP 

ptiAr^ATir 

LATITUDF 

-0 . 7 000/1 0‘3 

0. lOPl 57 1 

0.01'^0«'57 0. 000^^1 

LOrCl TLIDF 

105.01(^0500 

-0.0/!777«’0 

0 « 0 0 0 fi /t 0*^ 0.00003; 

GFOC PAf-GF 

007 l7.705-'i555 

“0.G500.5C 1 

1 .0 I'ipp 13 C .00/1^5 

CAL.1P U ‘T 

-7*5 5 l^./i7^.'i003 

-505 .73-^33<^3 

1 r , / . 1 7 1 7 

Tlf^F 

PAFGF DIF. .CAL. 

MA L , 

LOFGl TLIOFP (CFG) 


nr.T 

(FT ) 

.130.00.0- 

-53551 .50 

I ^ / 

-5355c. 50 

13003,0. 

-53555. ■p 

13003.3 

-53 5/' 3. PC 

13005.3 

-53553.50 

13007. 1 

-535/i3.P0 

.13007.3 

-5353A .30 

I 3003.''“ 

-53517.50 

1301 l.P 

-5 35/: 5. 70 

» 

c : 

-53553 . 50 

13010.5 

-5355 1 .50 

13013.3 

- 505-^0,70 

1 301/1. 1 

-53551 .50 

1 .301 A. 3 

-53 55 0./'0 

13017. 0 

-53550 .50 

1 301 ' ,P 

- 535/iP .70 

1301?. 

-53547. 7(^ 

13010,5 

_^q54/i. 70 

1 3 000. A 

-53553.50 

13001 . 1 

- 535/^7 . 70 

130;''! .o 

. -5Qr,/i/ , 7 0 

13000.7 

-535/17. 70 

13003.5 

-53557.50 

1300A .3 

-53550.50 

13005.0 

-5355P.C.O 

1 3005." 

- 5355.5 .50 

1 30.07. /• 

-53533 

13003.7 

-53535.33 

13030.5 

-535/1/: .70 

13031 .3 

-53544,70 


(FT) 

-rr'G 1 

GG 

-FGGGP'. G1 
-rrGGG. I ^ 
-?GGf’G.?r 

-7f’AGG 

^ O jT o 


.35 

.PG 


G3 

IP 

,3P 

, P ^ 
. oo 


-prGPG.^5 
-POGPl 
-PGGPl 

. -PPGPl 

-ppr.pi 

-PGGPP. I 1 

-PGIGP.2.5F 
-^?(^?.? • 3P 
“PPSP?. G3 

r,/»ro/*> 77 

^ r i> ♦ * * 

-PPGPP.^P 
-pPGp.3.r/‘ 
-?PG.r.5. 1 3 
-P PF P3 . 
-PPGP3 . 3f 
-ppGr.? 

PP 

-PPf-P-'i. PI 


73.3P7JQP I 
73 . Sp7rG73 
73.3G71 7PP 
73 . 3G7P3P.] 
73.3P71Q?P 
73.3G7P35P 
73.5073151 
75 . 30 7 ''■030 
7 5 . 30 70/:“^ 3 
75 . 50 7 1 r,F5 
75.3071501 
75.5070057 
75. 50 70.1/^ ft 

75.507 U70 
75.5071 
75 . 50 703 7 <3 
75 . 5070A^jP 
75.5070G^P 
75.5071 700 
75.5070/:75 
75.5070G75 
75 .5070^:'^ '' 
75.50710^3 

75.507 IP?«5 

75.507 1F3 1 
75 .50701 37 
75,3^73 1 1 0 
75 . 5073303 
75 .5070775 
75 .50707P5 


.50 7 050:' 

75 . 507 13r j 

/> 7 0 7 p. 

75. :)C7 0007 

.5073355 

75. 507 1301 

.50 73 1 5'’ 

75. 50700? 1 

.5070575 

75 

.5073005 

75.5070103 

. 5073773 

75.5070071 

,5074"3? 

75.5073513 

. 50 73 1 3 1 

75. .50705 1 5 

r.^fx 7 ^ 

75, 50 7''.'. 3? 

.5070143 

75.5070347 

.50 703 0 /* 

75.5070333 

.5070734 

75.5070305 

.507 0 ! 13 

7 5.50 7 00.G4 

0 7 1 

75.5070003 

.50 73.^03 

75.5070033 

.5073035 

75.507005 1 

.5073035 

75.5070073 

.5070357 

75.5070044 

. 5073 1 00 

75.5070055 

t f-i 7 7 T 0 

7 5. 50 ■^30.7 5 

.50 73 1 A 1 

75. ‘^07003^ 

.50 7 ''5 30 

75.507000.3 

.5070475 

75.507004C' 

.5070A7P 

75.5C70033 

C p T r* A 

7 5.507'''03.' 

. ^ sr?:^7r>r 

75 . 5:'7roro 

^ 7 7 

7 5.507030'.^ 

5 .50 734 03 

75,5070317 

>.50754.3'' 

75.5070330 




A-17 




TABLE A- 7 






% 




EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(NO SMOOTHING) (PAGE 2 OF 3) 




1 


TlrtJE 
6MT 
!3?3i>. ' 
\3T2,0.- 
13235.:' 
I iej 6 . j? 
132,".6.7 
li237.'i 

! 323£) . 1 


R/IWGE D1 r , 
(FT 1 

-5-' 647. 7:’ 
*S«'.667,/ir 

.s^(&7 9.5r 
-'^''653.^'' 
- 5-^72:. -r 


rfL. v-VL 
(FT) 


. r 'T P 
. c r. S 
— :*(“ T3 

.r ''S 

o ^ 


3 3 
.'I 7 


73. 

73. 

73. 

73. 

75. 

77. 

75. 


LorTiTur-^r (FIT) 

rcH’ir ciih’iL. Avn. 

3(‘7 1-'37 73 .3r7-’3;’S 7?.f’r7‘^.M7 


3T7^f;!- 7-.3;'7)''3r 

c f 7 7^ y o r 7 3 . 7 /. J, 7 r 

7 S. 3 r 7 ^5 3 7 

3r7i/.ii 73.5r7r:-3" 
3r7nir 73.3:'7i7 5r 
3 r 7 *■' /. 7 .3 7 3 . ‘ 7 3 I 

3 7 1 *' (■' 3 7 ‘' . 3 f ‘ 7 I f' ' 


7 7.~‘'7 

-7C ^ 7 - r 

7 c. C. !•• 7 I <j r, 
-? c, 7 .-■ .’ p’T 

73.3^7'’''*^'' 



1 

?:lP-5(’.^ 
31P51 .7 
J?IP5?.5 
?>1P54. 1 
2 I P 3 4 . p 
2 1 p 5 r. , 4 

2ir5F .r 

2l^5^^7 
2 If’ 39 .3 

2 if>ni ,9 

21022.3 
21024 .2 
21C25.2 
21«»03 .P 
2lf>23.3 

2102P.0 

21012.4 

?J o 1 1 , 2 
2 lo 12.F 
2 1 1 3 . 3 
2 J o I 4 . 3 
21015.0 
2 1 o i r. . 7 

21017.3 

21022.3 
21021 .4 


*301 

-32 1 13 .02 
-3212^.72 
-32134.32 
-32132.30 
-32133 .30 
-32 122. PO 
-32123.72 
-32 13^.32 

-321 33.52 

-30 13P.32 
-t;0144.4P 
-301 1 1.90 
-32 1 2.4 . 72 
-32 I 10.P2 
- 3 0 I 2 3 . 0 0 
-30 121 .P2 

-30144.40 

-30143,40 
-32135.30 
-32 1 15. '^0 
-32142.32 
-30 151.30 
-60 144,40 
-301 17. P2 
-30122. f'O 
-32144.42 
-30147.52 
-32157.50 
-30144.40 


-P 1091 

-9 \?^\ ,1° 

-r 1 0*= I 1 
-7 1 2'^2. 25 
-P 1092. IP 
-P 1092.45 
-F 1092.57 
_C 1 pop .P 4 

-f* 1205. 1 1 
-P 1 29s5.?2 
-p 1005,53 
-P 1 095,40 
-P 1095.76 
-f' 1 205. pp 
-P 1 204 . 1 5 
-P 1294. 2P 

-P 1 004.42 
_ p 1204,34 
-P 1094.37 

-p 1 O'^ 4 ,<34 

-F 1 203 .10 

- P I 2-9 5 . 5 5 
-p 1 P03 . 30 
-P 1005 .75 
-P 1P93.P5 

- f I 2.0 3 . 1 2 
-p 1203.23 
-P 1 293.50 
-P 1003.0 I 

-P 1097,04 


73 .3073025 
73.3273741 
73.5273224 

73.5074350 
75.5074^502 
75 . 507454P 
75.5075551 

75.5273523 
7 3 . 3 2 7 4 4 3f? 
73.3274305 
75.30744P5 
73.5074121 

75.50731OO 
75.50753f^3 
75.5273723 
73.3273312 
73.50753 IP 
75.3074 105 
75 .507423P 
75.5274Q0P 

75.5273174 

73 . 52744^? 
75 .32750P 1 
73 ,32742? 0 

75.5273075 

75.3273375 
75 . 52.745 17 
73.5074 1 45 

73.5074PO0 

75.3074375 


73.3273530 

73 .32'^'" 7 

? 5 . 52 733 40 
73 32732“^ 

75:5r72v." 
75.5275105 
75.52731*^3 
75.5073252 
73.50731 15 
75. 327324^' 
73.3073 1 5 1 
73.5074737 
75.3073P45 

73. 5:' 732.42 

73.3073572 
75.5077233 
73.3076234 
75.3074P41 
75.32.74Q I 4 
73.3073334 
73.3073C’ 1 o 
73 . 5075144 
75.5075727 
75.5074Q23 
75.5076610 
73.3076521 
75,3074065 

75 .30747?r 

75.3275446 
75.50750 10 


75.5073025 

73.5n"~7o 

73. 507*'" 3^ 
73.3073-27 
75. 3073''42 
75.5073 126 


73.507310? 

7 3.307321:’ 
73.3275120 

73.3073076 

73.3073225 

73.3074Q4<? 
73.3073244 
73. 307306Q 
73.32731 1 5 
73.3073206 
73.3073251 
75.5073175 
75.3073123 
75 . 30731 13 
75,30751 63 
75.5073155 
75.507515? 
75 ,5075007 
75.5075 152 
75.507 5 l‘^5 
75 .3073 12:’ 
75.30752-7 
75.5073077 


73.5073035 
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TABLE A-7. EPHEMERIS AND LOP SUMMARY FOR TEST 12 
(NO SMOOTHING) (PAGE 3 OF 3) 


TI !'T 

PfTF riF. 

PfL, VAL. 

LO^' 

PI TMOFf (PFC 

) 


(FT ) 

(FT ) 

POI r'T ] 

rOI f-’TP 

CUOUL. AUC, 

O 1 <> o n o 

-r.r i j r- . “-r 

1 P'^7 . 

7A . AP7/PPT 

7S , 9 7 A Apo 

75. 50750^9 

o j p p 

-^iP IT! . P P 

^ C j 

7A . A07AP0P 

7 r c i-> 7 c •- 7 

75 .5075P5T 

^ 1 . T’ 

-rp 1 1 p.^p 

IPo? . 

70 , A07A 1 Ap 

7 c, 79 7 f, 79 c 

75.50750PA 

n 

-^^P 1 />T , P 

-c’ J 907 . 7P 

7 A . A0 7/^/PT 

7A. AP7A 1 PP 

75.5075PA9 


“ <7 P 1 Ji =v 9 

-9 1 90 F ,09 

7A.507/;P?1 

7A. A!’7A/:AA 

75.507506 1 

1 0 94^’ , /l 

-r.p l/>p.^:p 

^r* 

7A, Ar?:^iAFT 

7A.AP7Ar.p9 

75. 50750^^9 

^ 1 Cj O ^ 

>FP |/!T. /;p 

-9 IPOp 

7 A . A0 7/iSTP 

75 . A07A1 77 

75. 5P75PT/’i 

9 1 O'"' 9 , 9 

- <.f 1 . S" 

-P ! PO'^ , / 7 

7 A . A 0 7/!9p7 

7A. AO 7^-9 7. A 

75. 5P75P IT 


-^P ITT.PP 

— ^ 1 P 0 9 A 

7A .5P7A 1 AO 

7 A, 50 7 Ar 1 A 

75 . 5P75P 1 7 

? 1 1 . 

-#^P 1 'S-': ,PP 

-r I P9p . 7A 

7A . AP7TP7T 

7A.5P7/:A IP 

75.5vP7/!i9pp 

1 

-(;p ip'.’.7p 

-9 1 poci ,p 1 

7A , AC’7AA97 

75 . AP7AT/.-T 

75.5P75PP5 

? 1 . f- 

- P 1 7 , T P 

-r 1 090 . 9A 

7A. AP7/:T/'7 

75 , AP7/;99p 

75 .507^1 'i'^0 

■< 

c 
: c 

_ rf 1 99 

_' • 1909,/; 9 

7A.AP7AAA) 

7 5. 50 7 A T9 7 

75 . 50 75005 

;■> 1 

-API//-. .' 1 T' 

1 999 , TT 

7A , AP 7/ A/o 

7 5,5C'7 5 [95 

75 , A 97 /:99 5 


- A, 9 ] /• A. , /i 9 

1999 . A? 

■'A . A 9 7 A A 7 f’ 

75 A 7 A 79 

75 . 5.97 /;99 p 


-^^PITT.r.P 

-r 1 p^^ p 

7A.5P7APAT 

75 , 5975999 

75 . 507^-^99 


-6P15 1 ,TP 

-0 1090,9/; 

7A . A07/. 1 /I 1 

75. A97A7P7 

75 .5074*^70 


-f?] 1 1 

-9 1 1 90.P A 

7A . A07AA/^ 

75.AP77P01 

75, 5975005 

p 1 Q 

- AP l/iS. /P 

n PP.T7 

7 A . AP7/A/.P 

75 . 507 A ]FP 

75 . 5P-7/i'^9 A 

P 1 ^ 1 . 7 

^ ^ * I r’ s 

-9 11 p'". A-A 

7A , 007/1999 

7 A 9 i'" 7 T, /: 

75 . 507^99f:, 

9 ] r; /; 9 , ‘=1 

-P.P 1 . 7P 

- 9 1 1 9 f' . A 9 

7A , AP7APp9 

7 A *' 9 7 9 / A S 

7 5 . A 0 7 A 9 1 9 

P 1 <=/-.^.P 

- Ap I /I 1 , /( P 

-5^ 1 IPP. 7 1 

7A. AP7/iFpT 

75 .507 A A AO 

75.507500P 

P [ 9 /: /; , 9 

- r P W: P . p P 

_ 0 1 1 99,99 

7A,AP7/iP90 

7A.AP7A5/:5 

75. 5P7500A 

PI 

- A P 1 p . P P 

1 IP 1 . AP 

7 A . AP7/0AP 

75.507/^7 1 

75 . 507/1990 

'-> 1 t>/;r" , ? 

-P'P 1 , 7p 

-91191.9/ 

7 A . A p 7 A9 /; 9 

75 .507AA99 

7 5. 50 7 5 00 A 

Pl^^P.T' 

- ^ P 1 r , 1 P 

-^ 1 1 PP. P^ 

7 A . A9 7Ti'iP7 

7 . A 9 7 /. 1 /; T 

75 , 507/!-9 79 

r ■ 1 9 J- , P 

-Ap ipp.7r 

-P 1 1 Op. AA 

7 A . AP7 AO -^'p- 

^ c s '*• 7 p 0 / 

75.A0 7/;99f; 

PI 

_ A, 9 1 p, 3 , fi 9 

-Pit PP,'^^ 

7A.AP7ATAA 

75!a07AP01 


P l . 1 

- n P 1 /• 7 . T P 

-F 1 1 9T . PP 

7A , A97/iAp7 

7A . A07AP7T 

7 5. 50 7/^<i9f; 


-r.plv^7,5p 

-f' 1 IP-T. A7 

7A.AP7/:A/jA 

75.A075PP1 

75.5074^-50 

''?rr ? . /■ 

- Ap 1 90 , Ap 

1 I OT .AO 

7A .AP7 A77A 

7A.AP7A/.P! 

75. A975PPT 

o p f ■ 9 ] 

-APIA) , A 9 

-^1 1 PT . 7 .': 

7 A.. A 0 7A9TP 

7A, AP7AAPT 

75 , 5075003 

PPPPP. ,P 

-rp l/ip./’p 

1 IPT.^’7 

7 A . A P 7 /f9 /: 

75 . A97AATP 

75,5075003 

PPPPP. 7 

-AP 1 af .TO 

-PI IPT.oo 

7A . A07/;A |P 

7A.AP7APA/1 

75, 507^"' 'i? 

PPPP.^.f^ 

- 6 P IT7.5P 

-p 1 1 P;4. IT 

75 . AP7AT ] ? 

75. 50750 AP 

75.507500? 

PPPP^i . 1 

-(^PJ5P.5P 

-P 1 IPv^. /IP 

7 3. A0 7/iA?P 

75,50751 AA 

75. 507^'2q5 

PPPPf> » 7 

- A, 9 1 /i J, . / 9 

-P 1 1 '.Vi . AA 

7 A , AP7/.9 7 7 

75 . AP 7 AApP 

7 5 . A P 7 A! 9 /; 

99pp7, /I 

*AP. j T-'i .ft 

-P I l{'/ . 7- 

7A,AP7AA/p 

7A , 5P7A 1 9" 

75.5075P0P 

?.?.???,?. 

-F,9 l/jP.^ P 

-? 1 1 ;v- .^P 

7A . AP7A 1 79 

7 5 . 50 759 p/; 

75 . 5075005 

ppppo , p 

-AP-ITo ,?0 

-P 1 lOA.PA 

7A, A07AP.A1 

75,5075PO7 

75.5P75P0P 

P.P 0 1 P . p 

-APIM ,TP 

-P 1 1 P-A ,TP 

7A . AO 7 / 1 APT 

75.5075 IA9 

75. 5075P0P. 

ppr I 1 

-f,r-l Af .TP 

-p 1 IP-5./<A. 

7A .A0747P1 

75.5P753A7 

75.507/lQ<5P 

P P J ^ C| 

- AP 1 ST.9P 

-PI IPA.7 1 

7 A . AP7/i/iTf 

75.5P75P7A 

75, 507/:C:99 

ppOK,.f 

- A .■* 1 /I . A 9 

-p nOA.TA 

?A . A07A 1 A 1 

75 .50 7 59 Of; 

75.5P7£Oop 

PPP 1 7 ,r 

-APIA 1 ./p 

-r 1 1 PA . Ap 

7A.5P7APTT 

75 . 50 759 70 

75.507/5«595 

pppir./ 

-AO 1 AF .PP 

^F 1 IP A. AT 

7A,AP7/!l 70 

75 .507/iPp^ 

7 5, 507 /i -59 5 




TABLE A“8. EPHEMERIS AND LOP SUMMARY FOR TEST 13 
(JULY 28, 1971) (PAGE 1 OF 3) 


ATS‘5 

EPHEMERIS 





TIME 

SAT RANGE 

LATITUDE 

LONGITUDE 


— 

GMT 

MM I 

DFG 

DEG 




2271G.GQ 

-0.502 

105.148 



5>RP5f?P ,Pl 

22717. 5P 

-0.413 

105.125 





2271«5.0G 

-0.316 

105.103 



3PI0P0.P 

22721 . I 0 

-0.214 

105.083 



33000.0 

227 23. G 7 

-0.108 

105.064 




PAPAPOLIC 

FIT TO EPKEMEPI S-LI NEAR TO CAL DATA 



CONSTANT 

LINEAR QUADRATIC 

RMS 

LATITUDE 

-0.7400286 0.14IG857 0 

.01 1 U29 0. 

00058 

LONGITUDE 

105.2260285 -0,0562857 0 

.0028571 0. 

00011 

GEOC RANGE 

22717.32000^^9 -2.I0400P4 1 

.1200017 0. 

008 49 

CALIP LIME 

397G5. 4700302 - 638 . 783«>655 

10. 

25137 

T1 ME 

RANGE DI F. 

CAL. VAL. 

LONGITUDES (DEG) 

GMT 

(FT) 

(FT) 

POINTl 

P0INT2 

CUMUL. AVG. ... 

P0020.7 

5<5425.75 

38484.2.3 

75.50721 18 

75.5072764 

75.5072118 

20205.4 

5Q40Q 1 

3P4r^5.65 

75.5072407 

75.5073053 

75.5072262 

2034^,2 

59382.72 

38447.23 

75.5071079 

75.5072625 

7 5.5072,168 _ 

20522. 1 

5Q365.30 

38450.75 

75.5071007 

75.5072553 

75.5072103 

20G4^ ,G 

5<>554.4» 

38415.22 

75.5072429 

75.5075075 

75.5072168 

208 3 G, 5 

59348.34 

38396.25 

75.5073367 

75.5074015 

75.5072368 

2 1 02P .2 

5^344 .0Q 

3^376.43 

75.5074485 

75.5075131 

75.5072670 

2I22R. 5 

5«>301 .04 

38355.00 

75.5075254 

75.5073900 

75.5072743 

2 1 4 1 G . 2. 

59283. 9G 

38335.98 

75.507351 1 

75.5074156 

75.5072828 

21G14. 1 

5^270.87 

38515,06 

75.5074148 

75.5074703 

75.5072960 

21P 1 6 . 7 

50231 .28 

38203.30 

75.5073165 

75.507381 1 

75.5072979 

22023.2 

59212.67 

58270.86 

75.5073561 

75.5074207 

75.5073027 

22235 .G 

50 102.72 

38247,72 

75.5073921 

75.5074567 

75.5073006 

22452.5 

501GG.54 

58223.07 

75.5073003 

75.5074659 

75,5073160 

22G38 .0 

50150. 17 

38204.35 

75.5074271 

75.5074917 

75.507J234 _ 

2280<^ .3 

50 1 19.63 

38 188. 15 

75.5073475 

75.5074121 

75.5073249 

22^23.5 

50102.89 

38 1 74.00 

75.5075340 

75.5073986 

75.5073255 

25047.0 

50089 . 19 

38160.17 

75.5073514 

75.5074160 

75.50732.69 , 

23205 . 1 

50074.43 

3P 146.31 

75.5073554 

75.5074190 

75.5073284 

23332. G 

50057.67 

38 130.79 

75.5073584 

75.5074229 

75.5073290 

23504 . 7 

50045.02 

381 14.44 

75.5073933 

75.5074578 

75.50.U329 

23G3R. 4 

50020.44 

38097.82 

75.5073546 

75.5074191 

75.5075339 

2375G.5 

50000.53 

38083.96 

75.5073851 

75.5074476 

75,5073361 

23^ 17.7 

59011 .12 

38069.55 

75.5074946 

75.5075591 

75,50734.27 

24022.5 

58078.07 

38058 ,05 

75.5073720 

75.5074365 

75,5073^3? 

24232. 1 

58956.47 

38055.06 

75.5073015 

75.5074560 

75.5073457 

2455G. 1 

58040.56 

38023.70 

75.5073708 

75.5074352 

75.5073466 _ 

24430.8 

58032.66 

38014.00 

75.5073892 

75.5074536 

75.507348 1 

24531 .7 

58919.24 

38003. 19 

75.5073803 

75.5074448 

75.5073492 

24G24 ,8 

58914.86 

37003.77 

75.5074191 

75.5074835 

75.507351.6 . 

24725.7 

58001 .24 

37982.96 

75.5074000 

75.5074735 

75.5073554 
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TABLE A-8. EPHEMERIS AND LOP SUMMARY FOR TEST 13 
(JULY 28 , 1971) (PAGE 2 OF 3) 


TIMK RANGE DIF. CAL. VAL. LONGI TUDFS (DEG) 

GMT (FT) (FT) POINT! POI NT P CUMUL . AVG. 

2AS52.P 5FP'"7.54 ,'S797l.0‘i 73.50740^0 75.5074^94 75.5073550 

24952.?. 5PP74.28 3796P.51 75.5073Qfifi 75.5074G11 75.5075563 

25055.1 5P866.5I 57Q4Q.71 75,5074257 75.50748?^? 75.50755P5 

25145.7 5FP51. 54 57056. R5 75.507419<^ 75,5074845 75.5075600 

25255.9 58856.50 57Q?4,57 75.5074125 75.5274770 75.5075615 

25405.^ 5pPl<i.|0 57912.50 75.5075P75 75.5074520 75.5075622 

25505.2 58810.65 57^01 *78 75.5074084 75.5074728 75.5075654 

25557.1 58810.18 57892.22 75.5074750 75.5075574 75.5075662 

25707.4 58795.15 37879.74 75.5074552 75.5075176 75.5075684 

25758.1 . 58775.49 57870.75 75.5074050 75.5074694 75.5075695 

25850.4 5877P.'^7 57861.47 75.5074414 75.5075058 75.5075710 

25^43.5 58755,55 57852.05 75.5074108 75.5074755 75.5075710 

50051.5 58742.41 57859. 'iP 75.507415 1 75.5074775 75.507572? 

50147.7 58757.01 57850.01 75.5074494 75.507513? 75.5075746 

30240.0 58725.50 57820.75 75.5074409 75.5075055 75.5075760 

50344.8 58722,44 57809.25 75.5075042 75.5075686 75.50757^7- 

50442.6 58717. 577?8.Q7 75.5075487 75.5^76151 75.5075025 

50549.7 58705.75 57787.07 75.5075428 75.5076072 75.5075o55 

5065P.2 586<^ 1 .51 57774.74 75.507551 5 75.5076157 75.507508? 

30757.0 58681.26 37764.48 75.5075605 75.5076247 75.5075922 

50905.7 58667.42 57752.29 75.5075590 75.5076234 75.5075954 

51002.7 58662.58 57742.18 75.5075988 75.5076632 75.5075995 

55937.5 58327. 7Q 57427,25 75.5077185 75.5077826 75.5074052 

54027.4 58510.54 57j518.40 75,5076709 75.5077551 75,5074 100 

5412Q.<i 582^5. PO 57407.31 75.5076464 75.5P77136 75.5074142 

54225.8 58284.85 57597.75 75.5076566 75.5077209 75.5074185 

54525.5 58271.77 57586.80 75.5076520 75.5077162 75.5074225 

54422.5 58259.84 37576.68 75.5076486 75.5077128 75.5074265 

54516.5 58249 . P8 57567.14 75.5076:)45 75.5077l?5 75.5074301 

54605.2 58245.25 57558.82 75.5076711 75.5077555 75.5074541 

54657.8 58255.54 37549.15 75.5076788 75.5O774'50 75.5074580 

54751,7 58224.57 57559.56 75.5076905 75,5077545 75.5074420 

34844.0 5K2II.54 57530.28 75.507674075.5077582 75.5074457 

54P37.1 58201.77 57520.86 75.5076795 75.5077457 75.5074495 

55027. P 58188.67 5751 1.P5 75,5076609 75.5077251 75.5074525 

55121.1 58178.45 57502.41 75.5076655 75,5077276 75.5074556 

35214, « 58169.75 37292.88 75.5076766 75.5077408 75.5074589 

35502.5 58157.75 57284.42 75.5076610 75.5077252 75.50746IP 

55555.6 58150.28 37274.99 75.5076812 75.5077454 75.507464P 

55446.5 58 159.90 37266.00 75.5076798 75.5077440 75.50746f?0 

55558.6 58127.37 57256.72 75.5076668 75.5077509 75.5074707 

55651.7 58120.69 57247.29 75i5076919 75.5077561 75.5074757 

55721.7 58106.89 57258.42 75.5076679 75.5077521 75.5074764 

55808.6 58101.55 5725O.10 75.5076937 75.5077579 75.5074795 

5‘^P54.6 560PO.96 37221 .94 75.5076850 75.5077492 75.5074p20 

55942.2 58088.24 57213.49 75.5077285 75.5077925 75.5074852 

4f'P2P,l 58074. 57205.17 75.5077056 75.5077678 75.5074^80 
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TABLE A“8. EPHEMERIS AND LOP SUMMARY FOR TEST 13 
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APPENDIX B 


RELATIVE POSITION DETERMINATION 


The equations derived in this appendix permit the determination of 
Lines of Position (LOP*s) which locate the movable remote receiver 
with respect to the fixed receiver at a known position. 

The receivers are located in areas for which Type 1 surveys are 
available. This permits determination of the distance from the 
fixed receiver to the satellite with minimal error (qoverned only 
by the uncertainty in satellite position) . 


LINE OF POSITION 

A Line of Position is defined to be a straight line containing 
the receiver's position. It is obtained by assuming a reference 
latitude near the true latitude. On the latitude circle, a point 
of longitudinal intersection is located. The same procedure may 
be repeated in order to establish a second point, following which 
the two points are connected by a straight line. Alternately, a 
line perpendicular to the line of sight to the satellite may be 
drawn through one of these points. In either case, a line is 
obtained which contains the remote receiver. A range reading from 
the fixed receiver and a range reading from the remote receiver are 
required for the determination of an LGP. 


THE GEOMETRY OF THE LOP 


The position of the remote receiver is defined by three quantities, 

R, and X, where R is the radius of the Earth, and and X represent 
the latitude and longitude, respectively. With one satellite avail- 
able, only one of these unknowns can be determined. It is mathe- 
matically convenient to select latitude and Earth radius as the 
known quantities and then calculate longitude. Latitude and radius 
are determined from the international geoid model. 

Since position on the Earth's surface is normally given in geodetic 
coordinates, the position of the fixed receiver as well as the 
geodetic latitude and Earth radius of the remote receiver are trans- 
lated to the geocentric coordinates shown in Figure B-1. 

If D is the distance between either receiver and the satellite, it 
follows from spherical trigonometry that 

+ Rg^ - 2RRg cos Y (B-1) 

where 


cos Y = cos ip cos cos (Xg - X) + sin sin (B-2) 


B-1 




Figure B-1. Geocentric Coordinate System 



In the sequel the subscript S refers to the satellite while sub- 
scripts A and B pertain to the fixed and movable receivers, respective 

The angle y defined in equation (B-2) is located at the center of the 
Earth and between the rays to the satellite and receiver under consid- 
eration. 

In this framework, the distances and Dg are measured ambiguously 
as and M^. The range readings Mn and are contaminated by 
receiver noise, atmospheric noise, delays in the equipment and 
atmosphere as well as the relative drift in phase, and frequency be- 
tween the two rubidium frequency standards used. In the absence of 
these error sources, the range differences Dg - and Mg - both 
measure the distance between the two receivers . In the presence of 
these error sources we still write 


Dg - Da = Mg - (B-3) 

where the equivalence is understood to be not mathematically exact. 
Since the errors induced by the rubidium standards are correctable, 
equation (B-3) may be rewritten as 


Dg = + Mg - + At + B (B-4) 

where the linear function describes the correction factor to the 
relative phase drift. The relative frequency drift, which would 
introduce a correction term of the form Ct^ , is negligible. 


The determination of the longitude of the remote site proceeds 
as follows. Since the latitude, longitude and geocentric range of 
the fixed site are well known and the position of the satellite 
is available from NASA ephemeris predictions, cos y can be calculated. 
Hence, Da can be calculated from equation (B-1) and the distance Dg 
from the remote site to the satellite follows from equation (B-4) . 
Since the geocentric range of the remote site is assumed known. 


cos yg 





(B-5) 


equation (B-2) is used to determine cos (X_ - X^) 
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cos (Xg - \g) 


(B-6) 


COS Yta ” sin 

^ D b 

COS COS 


Then, 

sin (Xg - Xg) 


±[1 - cos (X, 


- VI 


^2 


(B-7) 


and 

, sin (X„ - X„) 

X„ = X ± tan"-^ ^ ^ 

cos (Xg ~ Xg) 


(B-8) 


If west longitude is defined positive and increasing westward, 
the upper (+) signs are used in equations (B-7) and (B-8) . 


Equation (B-8) shows that Xg is completely determined in terms of 
and Rg, By changing and correspondingly, a new longitude 
can be calculated. The resulting point when connected to the first 
point completely defines the LOP. 

The LOP can also be specified by one of the points and the azimuth 
heading, which is given by 

Az = 3 + TT (B-9) 

where 

3 = sin 


sin 


(Xg - Xg) cos 


sin Y 


B 


(B-10) 


West longitude is taken to be positive in equation (B-10) . In 
either case, the line of position is completely determined. 

The accuracy with which Xg is determined is a function of the 
accuracy with which the satellite position is known. It is implicit 
in the structure of the solution given here that any error in sat^^ellite- 
position appears as an error in Xg. Estimates of the effect of in- 
accuracies in satellite ephemeris are treated elsewhere. 

The listing of the Lines of Position Program is given in Table B-1. 
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TABLE B-1, LINES OF POSITION PROGRAM (PAGE 1 OF 3) 


nir'FNsi or-' st ( ipd ,stt( ir> c ipd ,stlmc( ifi) , stgpa ck^) al( 2 ) 

CSLT (2) ,S NLT (2) ,S LT (4 ) ,S LG ( ^ ) ,?R/i (4 ) .CALPAR (^) ,C ALTr ( IPIPPI) ,CALDR (1000) 
FA(.^),FAN'(4),FAfi(^),AL0Pl( IPD.ALOPPCjOOl.AVGfJOP) DP( 10P) T^FS (100) 
CLVLdOO) 

* cnPVEPFION rOf'STANTG 

* FIXF.n SITF nFM'oTFP PY TFP!*iIMAL F IN VARIAPLE NAMES 

* PFFFPF.rCF LATITIIDFS OFMOTED PY TERMIN(5>L 1 OP 2 
DRAD :.017A532P25 

FTMM 1 1 5Apf, 

A L A T F =PP A D *3 7 . Q3 OQ 

ALFGFrpPAP*73.A73'^52 
CALL CONV(ALATF,ALTF,PF) 

RF=PF+<53.fil/FTNM 
ALAT 1 -PPAP*37.PG4253 
AL A T? = PP A 0*3 7.8 64300 
CALL COF\/(ALAT 1 , ALT 1 ,P { 1 ) ) 

CALL CGFV(ALAT2,ALT2,P (2) ) 

CSLTF-CfiS (ALTF) 

CFLT ( 1 )=COS (ALT 1) 

CSLT (2)rC0S (ALT2) 

SNLTFrSI r (ALTF) 

SNLT ( 1 )rSI N(ALT 1 ) 

STLT (2>-SI V(ALT2) 

* ACCESS FPHFMFRIS AND CALJPPATION DATA 

* TITLE STOPFP AS FAN 

* CALIPPATIOM DATA STORED AS FA 

lOr VPJTF (1,($ TEST NllMPEP AND DATE:^.?)) 

PFAP (1,(AA6)) FAN 

VPITE (I , (S OPEN FILE:S,Z )) 

^FAP (1,(3A6)> FA 
OPEN (3,1 NPIIT .FA ) 

1 PFAD (3,(15)) NSAT 

2 nr 3 I = I .NSAT 

PFAP (3,(4F20.0)) fT(I ) .STLATU ).STLNC(I ) STP 
STGPA (I ) = STP/FTNM 

3 CALL TCONV(ST (I ) ,S TT (I ) ) 

* FIT FPHFMEPIS DATA TO PARABOLAS 
CALL LSTSO (NSAT .2.STT .STLAT.SLT) 

CALL LSTSP (NSAT .2.STT .STLNG ,SLC) 

CALL LSTSO (NSAT.2.STT .STGPA ,SPA ) 

* DETERMINATION CF THE RMS OF THE EPHEMEPIS POINTS ACOUT THE INDIVIDUAL 

* PAPAPOLAS 
KrJ 

CALL PPEAC (NSAT.STT.STLAT.SLT.ALAT.K ) 

CALL ’^PEAC (NSAT, STT, STLNG, SLG.ALNG.v ) 

CALL PPEAC (NSAT ,STT , STGPA ,SPA .AGFA ,K > 
z PEA.P (3, (J 5) ) NCAL 
DO 5 I r 1 ,NCAL 

READ (3,(2F20.0)> CALT ,C ALDR (I ) 

5 CALL TCONV(CALT ,CALTM(I ) ) 

* FIT CALIPPATIOM DATA TO A STRAIGHT LI NF 
CALL LSTSO (NCAL , 1 .CALTM .CALDP ,C ALPAP ) 

CLOSF (3) 

* DETERMINATION OF THE PMS OF T)^E CAL POINTS ABOUT THE LINE 
Kr2 
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TABLE B“l, LINES OF POSITION PROGRAM (PAGE 2 OF 3) 


CALL PPFAC <NCAL .CALTM ,CALDP ,CALPAR ,ACAL ,K ) 
e JJrp! 

CNTP=P. 

’k ACCP?? FATFLLITE MF Af^UP FI^FNT DATA FI LF ^'AMF FAf« 

VPITF (1,($ CPFr-! DATA FILE:$,Z)) 

RFAn (l,(AAfi)D RAM 
OPFN (3,1 fJPUT ,FAM ) 

7 JJrJJ+1 

^FAP (3,(PFPP.P)) TMFF (JJ ) ,DP (JJ ) 

IF (TnF? (.1^1 ) .FQ .0. ) GO TO Pp 
CALL TCON'V(T[*iF? (JJ ) .TrFAS ) 

* INTFPPOLATF ATS-5 FPHEMFPIS AT TIME=TMFAS 
CLATD=SLT ( J ) + (FLT (PJ+SLT (3)*TMEAS )*T(*1FA? 

SLNGD:$LG< 1 ) + (SLG ( P ) +S L G ( 3 )^Tf*’EA,S ) ^Tf'EA? 

r PAMG:?PA < 1 5 + (PPA (P)+SRA (3 )*Tf«^FAF ) *Tf**FAS 

rLATPrPLATP^PPAD 

CL^’GP=SLMGP*PPAP 

CSSATiCOG (GLATP ) 

GNGATrSI NCSLATP ) 

C?FS:CSLTF*CSSAT*COS (SLNGR- ALNGF)+SMLTF*SNSAT 
+ EVALUATE CALIPPATIOP LI PE AT TI f^F =Tf*'FAS 
CLVL (.lU )=CALPAP ( 1 )+CALPAP {P)*TMEA? 

RPS=RE/GPANC 

PFS = GPA.NG*SOPT(PPG* + P+I . - P ,*RRG*C? FS ) 

PPGrPFG- (PP (JJ )-CLVL (tU ) )/FTMtn+^5/FTN[« 

* DETERMINE LONGITUDES FOR THE GIVEN LATITUDES 
no P 1 1 1 , P 

f^PrP (I )/51p A NG 
DQ= DPS /SPANG 

COS J - (POjk*P+l .-no**P)/{p,*RO ) 

rOSDL: (COS 1-SNLT (I )=fcSNSAT) /(CSLT (I )*CSSAT) 

n NDLrSQPT ( 1 .-C0SDL**2) 

P AL (I) = (SLNGP-ATAN (SI NDL/COSDL ) )/DPAD 
ALOP 1 (JJ ):AL ( 1 ) 

ALOPP(JJ )r AL(2) 

CNTR=CMTR+AL( ] ) 

AVG(JJ)rCNTP AM 
GO TO 7 

+ GO TO 7 FOP MOPE MEASUPEMENT DATA 
PO CLOSE (3) 

* OUTPUT 
IPGrl 

WPITE (1,P1) FAN, IPG 

PI FOP MA T ( " /20X , 4 A G , 20x , iPA GE S , I 2 / ) 

WPITE (1,P2) 

P2 FORMAT( lOX, TATS-5 FPHFMFRI S S/l OX , T>TI ME $,3x , SSAT RANGE LATI 
TUOE LONGITUDE 4/ IPX.S GMT S , fix , XNMI I , PX , iDEGS ,PX , $DEG S) 

WPITE (1,01) (ST (I ) ,STGPA (I ) ,STLAT(I ) .STLNGd ) ,I = 1 ,NSAT) 

01 FOPMAT(fiX,Fo. 1 ,4X,F8 .2,F1 1 ,3. FI 2.3) 

WPITF (1,02) (SLT (I ) ,Ir 1 ,3) ,ALAT , (SLGU) ,I = 1 ,3) ,ALNG, (SRACI ) .1 =1. 3 ) AGPA 
,CALPAP ( I ) .CALPAR (2) ,ACAL 

02 F0RMAT(//5y,T PAPAPOLIC FIT TO FPHFMER I S-LI NFAP TO CAL DATA*^/ 

2 3X , $C0NSTANTS.7X , SL I NFAP 3. , 4X , SOUA PR AT I C S , 5X , $PMS $/ 
€X,$LATITUPFS,4X,3F13.7 ,FP,5/6X, SL0NGITUDES,3y ,3F1 3,7 ,FP.5/ 

6X,5iGFOC PANGFS2X,3FI3.7,FP.5/GX,K:ALIB LI NE S, IX ,F 1 4 . 7 ,F 1 3. 7 1 3v FP.5) 
WRITE (1,03) 
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TABLE B~l. LINES OF POSITION PROGRAM (PAGE 3 OF 3) 


FOPMATC //inx, $nME RANGE PIE. CAL. VAL . S,oy , $LONGI TUDES (DEG) $/l Ply , 
% GMTS.5X, J(FT)S, 7X,5.(FT)$,?X, SPPINTl S,5X,$P0INT?.$, 1X,$CUMUL. AVG.S) 
LlfJiNCAT+l^ 

DO ^7 Y~l 1 

WPITE (1.95) TFF? (K ) ,PP (K ) .CL^'L (K ) ,AL0P1 (K ) ,ALOP2(K) ,AVC(K) 

05 FORMAT (5X.FP. 1 ,F I 2.? ,F 1 1 .2 ,2F U . 7 ,F 1 2 . 7 ) 

LI N = LI N+1 

IF (LIN-50) o7,t)7,Ofi 

06 IPGrJPG+l 

VPITE (l.PD FAN, I PC 

U'PITE (1.95) 

LIM:A 

9 7 CONTINUE 

WRITE (!,('')) . . . 

CO TO 100 
DTOP 
E 

CUPPOUTINE CONV/(AA ,PP ,P ) 

* THI9 SUPPOUTI NF CALCULATES GEOCENTRIC LATITUDE AND EART 'S RADIUS 
FTNMr6076. 1 154R6 

PAr20oj>574 1 .00/FTNM 
PP =20«5550 I .00/FTNM 
E2- 1 .-(RP /PA )**2 
A rF2/(P.-E2) 

T UA*?IN(2.*AA ) 

T2-A*A*SI N(4.*AA )/2. 

T 5zA*A*A*SI N (fi.*AA )/5 . 

PP =AA-T 1+T2-T5 
P rpp /SORT( 1 ,-F2*C0S (EF )**2> 

RETURN 
E MD 

SUFPOUTIME TC0NV(T1 ME ,TMD) . .. 

* THIS SUPPOHTINF CONVERTS INPUT TIME TO DECIMAL TIME 

11 rTlM.F/10000. 

RES-TIME-I 1*10020. 

12 =RES/1P!0. 

RET=RES-I 2*100. 

TMD = I 1+(I 2+RET/GO. )/60. . . 

RETUPN 

ENP 

CUPROUTINE PPFAC (NPNTS ,TI ME .OFfVAL ,FI T ,CUT ,L) 

* THE PMS APOLIT THE MOVING MEAN IS CALCULATED IN THIS SUPROUTINE 
DI MENS I ON TI ME ( ^00) .OFFVAL (^00) ,FI T(5) 

DIF:0. 

no 2 T^I.MPNTS 

VAL:0. 

no 1 U- 1 .L 

1 VAL;VAL+FI T(J)*(TIME (I )**(J- 1 ) ) 

2 DI FrPIF+(VAL-OPSVAL (I >)**?, 

OUT=SOPT (DI F/NPNTS ) 

RFTIIPN 

END 
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APPENDIX C 


SIGNAL DESIGN 


The ALPHA II receiver uses correlation between the received 
signal and a locally generated signal replica in order to reject 
the noise on the signal and enhance the accuracy of the range 
determination. The correlator has excellent noise characteristics, 
as developed in Appendix D. The signal design is therefore 
dedicated to the rejection of multipath while retaining the 
inherent characteristics of the correlator-receiver. 

For these reasons pseudo-random codes were selected as signal 
candidates. The most important properties of these codes* are: 

(1) The niimber of positive bits in one length of the periodic 
sequence exceeds the nur.iber of negative bits by exactly 
one; 

(2) A code sequence multiplied by a shifted replica of itself 
is the same code sequence shifted to a further location; 

(3) Pseudo-random (PRN) codes are conveniently generated 

by a linear shift register. If the shift register is N 
stages long, codes of length 2^-1 can be generated. 

Since the essential difference between the various pseudo-random 
codes is the number of bits L (= 2^-1) before the code repeats, 
time sidelobe rejection can be obtained by selecting a code of 
sufficient length. 

The autocorrelation function of some function f(t) is defined by 


r(t) = Lim 


f(t) f (t-x) dt 

2T -T 


Relation (C-1) can be written as 


(C-1) 


R(T) 


1 

L 


L 

^ f(t) f (t-T) dt 


(C-2) 


*For these and the subsequent discussion reference, e.g., Berkowitz 
Modern Radar, Analysis, Evaluation and System Design, Wiley, 1965. 



for maximum length codes, since these codes are periodic. The 
value of the autocorrelation function when t equals zero, and the 
two sequences are aligned bit by bit is thus 

R(0) = i f^(t) dt 



(C-3) 


The amplitude f ^ (t) equals unity for all times t since it is the 
bit-by-bit products of +1 and +1 or -1 and -1 amplitudes. If t 
is not equal to zero 


R(t) 


1 

L 


L 

J f(t) f (t-T) dt 


= i jl. f(t+T) dt 
li u 


(C-4) 


by the second property. Property one states that there is one 
more positive bit than there are negative bits . Hence , the time 
sidelobes described by R(t) are the residual of this one positive 
bit over the total number of bits in the code, or 


R(t) = 



(C-5) 


The autocorrelation function is as follows: 



The value 1/L is assumed at T^, the length of a single bit, and 
will remain at this level through the (2^-2) -th bit. If the code 
is infinitely long, the autocorrelation function will begin to* 
repeat itself assuming the value unity at the (2 -l)-th bit. If 
the code is one code length long, R(t) is zero beyond the (2 -2)-th 
bit. 
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and is shown 


The power spectrum, S (w) , has the form 
schematically in the figure 




It may be observed that ^he autocorrelation function and power 
spectrum are identical to those of a pulse of width 

The codes shown in Table C-1 are maximum length pseudo-random 
codes generated from a 14-bit shift register. Of course only the 
first 9 stages are used in feedback connections in order to 
generate the code of 511 bits in length. 


TABLE C“l. CHARACTERISTICS OF MAXIMUM LENGTH 
PSEUDO- RANDOM CODES 


Code Length 

ATS -5 

Search Time 

*Multipath 

Rejection 

Comments 

511 

35 secs 

-54 dB 


1023 

65 secs 

-60 dB 


2047** 

135 secs 

-66 dB 

Fundamental ambiguous 
search code 

16,383** 

17.4 min 

-70 dB 

Code for unambiguous 
motion of satellite 


*Multipath rejection obtainable on the basis of the codes 
structure only and in a 50 Hz bandwidth 

**2047 or 16,383 are available in ALPHA II system 
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APPENDIX D 


CORRELATOR ANALYSIS 


A matched filter is a device which treats a signal immersed in 
noise in such a way as to maximize the output signal with respect 
to a fixed output noise level. 


A correlation system acts as a matched filter. If the modulation 
signal for navigation is matched and the modulations due to path 
length dynamics and for communications are not matched because 
they are not known a priori, then these modulations remain at the 
output of the matched filter and the effective bandwidth should 
be made to satisfy the relation 


H (iw) I 


S(w) 

S{w) + N(w) 


(D-1) 


where |H(iw) | is the magnitude of the filter characteristic, S (w) 
is the power spectral density of the signal including modulation 
due to path dynamics and the communication signal, and N (w) is the 
noise power spectral density. The expression gives minimum error 
in reproduction of the residual modulation with noise present. 


In the All ALPHA II receiver, the response approaches the opti- 
mization criteria and, hence, has optimum signal-to-noise behavior. 
The correlator functional block diagram is shown in Figure D-1. 


I-Channel 



Figure D-1. Correlator Functional Block Diagram 
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The basic correlation system is a matched filter for the incoming 
signal. The transmitted bi-phase PSK signal is given by 

Vq = B f (t) cos 03 q t (D-2) 

where B is the signal magnitude at the source, f(t) is the 
modulation and has values of +1 only, and Wq is the radian 
carrier frequency. 

This signal is delayed by the propagation time, reduced in 
amplitude during transmission and is translated in the receiver 
to the intermediate frequency 

The signal applied to the correlator is 

= A f(t - t) cos t - u)q t) (D-3) 

This signal is contaminated by noise, n(t), which may be written 
with in-phase and quadrature components as 

n(t) - cos t - co^ t) + ng sin (wi t - x) (D-4) 

The n^ und components are independent variables of time and may 
be assumed to be normally distributed with spectral density 


Without change in rotation, the signal plus noise will also be 
denoted by 

The delay, x = x (t) , contains the doppler terms. Let 


X (t) 


R(t) 

C 


(D-5) 


where R(t) is the distance the wave travels between transmitter 
and receiver and C is the speed of light. Then, 


U) 


(Oq x(t) 


R(t) = 


R(t) 


Expanding , 






2tt 


2tt 


1 

2 

+ 



(D-6) 

(D-7) 

(D-8) 
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The terra R(0) is the initial range, ^ R(0)t is the doppler 

term whose corresponding frequency is -5— and ^ R*( 0 )t 2 is the 
acceleration terra. o o 

The signal is multiplied in the correlator I~channel balanced 
raodulated by the locally generated modulation function, F 2 = f(t-T). 
The balanced modulator perforins the function described by: 

^ (t-T) (D-9) 

or 

= A f(t”X)f(t-T) cos t - t) + n^ f (t-T) x 

cos t - 00 ^ t) + ng f (t-T) sin t - t) 

(D-10) 

Similarly, the Q-channel balanced modulator multiplies with 
V 3 , where 

V 3 = f (t-T + TQ/2)-f (t-T - Tq/2) (D-11) 

The output becomes 

V 5 = (f(t-T + Tq/ 2) -f( t-T -To/2)] (D-12) 

or, 

V 5 = Af(t-T) [f(t-T + Tq/ 2) -f(t-T - To/2)] X 

cos (wi t - ( 1)0 t - ^o t) + 

ng sin (w^ t - Wq t)1 x [f(t-T + To/2) - 

f(t-T - To/2)] (D-13) 

The output voltages (V 4 , V^) from the balanced modulators are 
filtered by the bandpass filters, whose outputs (Vg, V 7 ) respectively 
become approximately 

Vg = i^f (t-Yj f (t-T) cos (Wj^ t - Wq t) + (D-14) 

•Where the noise terms are 

~ ^ ^ ^2R ^^1 ^ ” ^o^^ (D-15) 

= n^ f (t-T) , (D-16) 

q 2 j^ = Ug f(t-T) (D-17) 
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In the difference channel, the signal plus noise assumes the 
form 


= A f(t-T)[f(t-T + Tq/ 2) - f(t-T - To/2)!| cos 
(“l t T) + 

where 

"a " “^lA t - Uq + 'J2A ^ " “o 

while 

= n^ [f(t-T + T^/2) -f(t-T T^/2) J 
= ng [£(t-T V T^/1) -i (t-T - T^/m 


(D-18) 

(D-19) 

(D-20) 

(D-21) 


Averages are indicated by superposed bars. Autocorrelation 
functions R (t-T) and (t-T) are 

R (t-T) = f (t-T) f (t-Ty 


r^(t-T) = f (t-T) [f (t-T + T^/i) -f (t-T - To/2yT (D~23) 


The noise and code are .independent. The spectrum of the code is 
flat while f2 equals unity. The following spectral densities 
result in the single sided IF bandwidth of Bx/a 




(D-24) 



(D-25) 


The expressions for Vg and V 7 may be rewritten as 

= A R (t - T) cos (w^ t - u)qT) + n^^ (D-26) 

V? = A (t-T) cos (Wj^ t - x) + (D-27) 

Note that the cosine terms contain the received phase information 
including doppler. 


The output available at point V. can be used to measure doppler 
when lock-on occurs (x = T) . The error signal is Vg and is given 
by 

Vg = Vg . V 7 = a2 R (x-T) R^ (x-T) cos2 t - 0 )^ x) + 

noise terms (D-28) 
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After filtering in the low pass filter the error voltage resulting 
from the first term of Vg becomes 

V 9 = i R (t-T) • (t-T) (D-29) 

This is used to control the phase of the tracking code generator. 
Near lock-on, we have (x T) 

Vg = I (t-T) (D-30) 

for the error signal output and AR(0) for the doppler signal 
output. 

The lock-on point is chosen such that the autocorrelation function 
nulls at T = T and is linear. With pseudo-random code modulation, 
two decoders are required in the Q-channel. 

The receiver noise results from the noise terms indicated in the 
expression for V 3 . These noise terms are 

^10 2 ^IR ■'* 2 ^ (t-T) + 


^"^lA ^IR ^2A ^2R^/^ (D-31) 

The bar now indicates averaging in the low pass filter bandwidth 
B 2 * In order to determine the error induced by the noise, the 
noise power must be determined. Since the q parameters are normal, 
zero mean and each is assumed independent of the other three, 
equivalently to the determination of power, the variances are to 
be determined. 

The quantities q, ^ and q . . are both flat over the IF bandwidth 
hence over tne bandwidth B2 they contribute 2 Nq B2 and 
4 Nq B 2 f respectively. 

The power associated with q^^ q^^j^ can be determined by convoluting 

the power spectral densities of q^^ and Since both are flat 

over the same bandwidth B^^, the result is a triangle with base- 

2 

width and magnitude at zero frequency of 8 B]^. Thus, the 

variance of q~^ qj^j ^/2 (and because of symmetry also q^^ q^^/2) is 
8 Nq^ B^B 2 / 4 . Those two terms together thus contribute 4 Nq^ ®1®2* 


D-5 



The variance of V 


becomes thus 


10 


II 

( 1 

T 2 

0 iq^^) 

2 

+ [|r(t-T)] o2(q^^) 

4 . 

^2 1 ^lA ^IR 

^ '^2A 

<J2R 


^ L 

2 

2 

J 

r ,2 

II 

(T-T)J 

2N^B« + 
O 2 

[1 R(t-T)J 4N^B2 

+ 



(D-32) 


The noise voltage will cause a timing error At as can be 

seen from the figure (Z-curve) shown below* 



Since near lock-on (t=^T) , R(t-T). « 1 . while R^(x-T) is equal to 
unity when the error correcting signal has to correct for a 50 
nsec shift 


At 


50 


A ^/2 


Hence, 

At 

The variance is then 

o ^ 

a ^ = (100)^ (D-33) 

A 


AV, 


= 100 
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Substituting the expression of into this result and taking 

the square root 


= 100 


(T-T) 

2a2 


N B . 


4N 


^ 1^2 


1/2 


(D-34) 


J 


near lock-on. Introducing the bitwidth W, which is 100 nsec, i.e., 
twice the 50-nsec shift corresponding to maximum signal error 
voltage, then 


W 

2 


— r 2 ^ ^ 

Ti 2 a « 2 


16N 


^ 1^2 


1/2 


(D-35) 


The signal power, S, is A^/2, The noise power N depends upon the 
way bandwidth is defined. In terms of the audio bandwidth 


N = 

O 2 



(D-36) 


Under the assumption of a sufficiently large signal-to-noise ratio 
in the IF bandwidth (considered subsequently) , the third term 
under the square roof^raay be neglected. 

The magnitude of the first term would be 50% the magnitude of the 
second term if the maximum error voltage, a 2/2, were permitted, 
corresponding to a displacement of 50 nsec. 


I (t-T) = I at t-T 


50 nsec 


(D-37) 


In actual operations, a four phase clock is used. The effect of 
this clock is that range measurements are only made when the 
signal plus noise voltage is within +12.5 nsec of the nulling 
point, O, along the slanted line. Hence, 


N 

S 


(t-T) 



A- 

16 S 


(D-38) 
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The maximum value assumed by this term is 1/3 2th that of the 
second term and is thus neglected. Hence, 



(D-39) 


W 



if this second term is negligible. 


This is so when 
2 


16 N 


T 




<< 


O 2 
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which results from equation (D-35) - Simplifying, 


4N B 
o 


1 



« 1 


or 



>> 2 


This is the signal-to-noise in the IF bandwidth- The measured 

S/N in the 1600 Hz IF bandwidth at NASA Wallops was 6 to 10 dB or 
4 to 10 numerically. The indicated assumption is thus justified. 

If instead of the audio bandwidth, a double-sided bandwidth is 
used in tracing the S/N from IF to the filter in which the 
measurements are made, equation (D-39) assumes the form 



(D-41) 



APPENDIX E 

EFFECTS OF BAND LIMITING ON CORRELATION RECEIVER PERFORMANCE 


A signal of the type 

T (t) = A f(t) cos t (E-1) 

is transmitted through the satellite transponder which band limits 
the signal. The consequence of this limiting is a modified signal 
defined by 

(t) = A g(t) cos (jO^ t (E-*2) 

where f (t) is the original code and g(t) is the modified code 
due to bandwidth restraints. The received signal, after delay of 
T seconds, is 

R(t) = g(t - t) cos ’(t “ t) (E-3) 

The correlation receiver provides an output amplitude of 

|s| = f (t-T) g (t-T) * (E-4) 

where the bar over the product of code and bandwidth-modified 
code implies averaging over some time interval (integration time) . 


The frequency transform of the product, f (t-T) g (t-x) , is the 
convolution of the two signals in the frequency domain. Setting 
T = T and transforming f (t-T) = F(w) and transforming g(t-T) = G (co) , 

00 

F(V) * G(w-V) dV (E-5) 

G(cd) = H(oj) F(w) (E-6) 

where H (w) is the band characteristic function of the satellite 
transponder. Hence, 

P(iw) =_J[,“f(V) * F(o)-V) H(oi-V) dV '(E-7) 
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The average signal amplitude output, equation E-4, is the D-C 
{o)=0) spectrum of equation E-7. 

|S| = P(0) = F(V) * F(-V) H(-V) dV (E-8) 

F(w) is symmetrical for a maximum length PN code and 

|s| = P(0) = |F(V)!^ H(V) dV (E-9) 

s. 

The largest value |s| can take is for H (V) = 1, then, 

I^Lax = |F(V)|^dV (E-10) 

If H (V) = 1, when V ^ 0 and H (V) = 0 otherwise, then 

|S| = |f(V)|^ dV = 1/2 |s| (E-11) 

' ' 0 'max 

A loss of half the spectriam results in a signal amplitude reduction 
to half its value and the signal output power is reduced by 6 dB. 
This is the case for the ATS-5 C-band to L-band cross-strap mode 
and the 10 megabit ALPHA II PN code. 
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APPENDIX F 


AN ALTERNATE METHOD OF LOP COMPUTATION 


An alternate method of computing relative position is described 
here . 

*+- 

The quantity Ar is the magnitude of the vector AR which is 
directed from the fixed receiver to the satellite. The geometry 
involved may be understood by referring to Figure F-1. The 
distances and Rg are radii of concentric spheres and S£. 

If the elevation and azimuth of the satellite are known as 
functions of time, the line of sight (LOS) from the fixed re- 
ceiver A to the satellite can be established. The distance R 
is measured along the LOS and the surface through the end point 
C perpendicular , to the LOS is established. To a first order 
approximation, this surface is the sphere Sb since the radii of 
curvatures are extremely large, and thus, virtually planar. The 
line CD may be drawn in the "plane” Sb locating D on the surface 
of the Earth. The perpendicular to AD on the surface of the 
Earth is then the LOP containing B. The origin of some x-y 
coordinate system may be placed at the fixed site and the observed 
lop’s described in this coordinate system. 

One could then introduce the known location of the fixed receiver 
on the surface of the Earth which implies a transformation of 
coordinates from the fixed site to the center of the Earth. If 
the azimuth and elevation of the satellite at the fixed site are 
known AR can be determined in this Earth-centered coordinate 
system. The terminus of AR can then be expressed in longitude 
and latitude. 

Knowledge of the satellite position and the location of the fixed 
site is equivalent to knowledge of the azimuth, elevation and 
geocentric coordinates of the fixed site. Using spherical trigo- 
nometry, equations can be determined relating these two sets of 
parameters. In either case the LOP ' s are determined relative to 
the fixed site. In both cases satellite position error is present 
either directly or indirectly. 

Since the purpose of the relative ranging tests at Wallops was to 
determine range and thereby to accurately measure position, the 
relative separation between the two sites had to be known precisely 
as a check on the LOP calculations. Therefore, sites were selected 
for which a first order survey was available. In addition, the 
orientation of the range difference vector in three dimensional 
space needed to be known with high precision. This information 
can only be obtained from measured elevation and azimuth angles 
or equivalently from satellite ephemeris and fixed site coordinates. 
This information permitted determination of the LOP's for the 
remote site. 
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Figure F-1, Local Line of Position Geometry 



